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Prefaţă 



Prefaţă 

Cartea a fost scrisă pentru aceia care sunt pasionaţi de acest domeniu 
îngust dar foarte vast al dinamicii rotorilor: modelarea dinamică a rotorilor. 

In speţă, cartea se adresează studenţilor masteranzi, doctoranzilor şi 
specialiştilor în domeniu. 

Această carte oferă cititorului atât istoricul dezvoltării acestei ramuri a 
ştiinţei, cât şi unele noutăţi - la momentul scrierii cărţii - în ceea ce priveşte 
punerea problemelor şi rezolvarea acestora. 

Cartea continuă preocupările mele începute în teza de doctorat având 
aceeaşi temă şi sunt recunoscătoare profesorului Mircea Radeş, conducătorul 
meu de doctorat, pentru ajutorul şi îndrumările oferite, pline de profesionalism şi 
de talent. 



Dana-Codruţa Vişan 
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Simboluri şi notaţii utilizate în lucrare 



2 

A s = %A aria efectivă de forfecare a secţiunii arborelui, [m ] 

E modulul de elasticitate longitudinal, [N/m 2 ] 

EI modulul de rigiditate la încovoiere 

G modulul de elasticitate transversal, [N/m ] 

GA S modulul de rigiditate efectiv la forfecare 

i , j nodurile care definesc capetele elementului de arbore 

I = I y = I z momentul de inerţie axial al secţiunii transversale, [m 4 ] 

I P =2I momentul de inerţie polar al secţiunii transversale, [m 4 ] 

J P , J T momentul de inerţie masic polar, respectiv transversal, 

[kg -ni 2 ] 

EI 

K = 12 constanta de forfecare 

GA S L 

L lungimea elementului de arbore, [m] 

r raza interioară a elementului de arbore, [m] 

R raza exterioară a elementului de arbore, [m] 

u deplasarea laterală în lungul axei Ox 

v deplasarea în lungul axei Oy 

w deplasarea în lungul axei Oz 

y xy deformaţia de forfecare (lunecarea specifică) în planul x-y 

y xz deformaţia de forfecare (lunecarea specifică) în planul x-z 
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fi masa pe unitatea de lungime, fi=pA 

^ masa de rotaţie pe unitatea de lungime, ji = pi 

v coeficientul lui Poisson 

£ = — coordonata adimensională 
* L 

p densitatea materialului arborelui, [kg/m 3 ] 

(p rotirea în jurul axei Oy 

/ factorul de forfecare 

v|/ rotirea în jurul axei Oz 

co pulsaţia proprie de precesie 

Q viteza unghiulară a rotorului, [rad/sec] 

K xy rotirea specifică în planul x-y 

K xz rotirea specifică în planul x-z 

[c b ] matricea de rigiditate pentru lagărele cu alunecare 

jf y s j subvectorul forţelor nodale pentru elementul de arbore 

(şhaft) în planul x-y 

|f z s j subvectorul forţelor nodale pentru elementul de arbore 
în planul x-z 

jf s } vectorul forţelor nodale pentru elementul de arbore 

|f y d j subvectorul forţelor de dezechilibru din planul X-Y pentru 
disc 
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jf z d | subvectorul forţelor de dezechilibru din planul X-Z pentru 

disc 

jf d } vectorul forţelor de dezechilibru pentru disc 

[g s ] submatricea giroscopică pentru elementul de arbore în 

planul x-y sau x-z 

[g d ] matricea giroscopică pentru un disc axial simetric 

[g s ] matricea giroscopică globală pentru elementul de arbore 

[k s ] submatricea de rigiditate pentru elementul de arbore 

în planul x-y sau x-z 

[k b ] submatricea coeficienţilor de influenţă pentru lagărele cu 

alunecare 

[k s ] matricea de rigiditate globală pentru elementul de arbore 

[m s ] submatricea de masă pentru elementul de arbore în 

planul x-y sau x-z [m s ] = [mţ ]+ [m^ ] 

[mţ ], [m| ] submatricea de masă de translaţie, respectiv de rotaţie 

pentru elementul de arbore în planul x-y sau x-z 

[m s ] matricea de masă globală pentru elementul de arbore 

[m d ] matricea de masă pentru un disc axial simetric în planul 

X-Y sau X-Z 

|_Nj funcţiile de formă pentru translaţii 

|_N J funcţiile de formă pentru rotiri 
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|_nJ alt set de funcţii de formă [nJ= [n_|- — ^ 

{q} vectorul coordonatelor generalizate 

{q} vectorul forţelor generalizate 

|u y s j subvectorul deplasărilor nodale pentru elementul 

de arbore în planul x-y 

|u z s | subvectorul deplasărilor nodale pentru elementul 

de arbore în planul x-z 

|u s ] vectorul deplasărilor nodale pentru elementul de arbore 

|u y d j subvectorul deplasare din planul X-Y pentru disc 

|u z d j subvectorul deplasare din planul X-Z pentru disc 

|_%J funcţiile de formă pentru deformaţiile de forfecare la 

elementul CON 12. 
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1. Introducere 



1.1. Definirea problemei 

Tema lucrării de faţă este modelarea dinamică a sistemelor rotor- 
lagăre, temă izvorâtă din necesitatea practică de proiectare şi analiză dinamică a 
maşinilor cu rotor şi din nevoia de cunoaştere şi interpretare a răspunsului 
dinamic al acestora. 

In modelarea dinamică a rotorului, lagărelor şi a structurii de 
suport, se ţine cont că acestea lucrează ca un sistem, răspunzând împreună 
solicitărilor dinamice, intercondiţionându-se reciproc. Rotorul este o parte 
integrantă a unui sistem dinamic, comportarea lui fiind determinată de poziţia şi 
rigiditatea lagărelor, a etanşărilor, a piedestalurilor şi a fundaţiei, precum şi de 
amortizarea acestora. 

Modelul de calcul necesar studiului dinamic al sistemelor rotor- 
lagăre se poate realiza prin două metode: 

1 . Metoda matricelor de transfer; 

2. Metoda elementelor finite. 

In lucrarea de faţă se foloseşte metoda elementelor finite, pentru 
posibilitatea mai mare pe care o are de modelare a tuturor elementelor, efectelor 
şi interacţiunilor dintre părţile componente ale maşinilor cu rotor, spre deosebire 
de posibilităţile limitate de modelare pe care o are metoda matricelor de transfer, 
în aceleaşi situaţii. 

Corespunzător metodei de calcul alese, pentru modelarea dinamică 
a rotorilor există o mare diversitate de elemente finite cu care se poate realiza 
modelul: elemente brick, elemente shell, elemente de tip bară cu două noduri la 
capete etc. 

In lucrarea de faţă, arborele este modelat cu elemente finite de tip 
bară cu două noduri la capete, de tip C 1 . 

In Catedra de Rezistenţa materialelor din Universitatea Politehnica 
Bucureşti a fost elaborat un program de calcul de dinamica rotorilor, în mediul 
MATLAB. Faţă de programul existent, lucrarea de faţă aduce îmbunătăţiri, 
completând elementul de arbore de tip Bernoulli-Euler cu care lucrează acesta, 
cu elementul de arbore de tip Timoshenko. De asemenea pentru modelarea 
arborilor cu porţiuni tronconice se introduc două elemente finite conice de tip 
C 1 . 
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1.2. Scopul lucrării 

Scopul lucrării a fost elaborarea unui model îmbunătăţit pentru 
calculul răspunsului dinamic al sistemelor rotor-lagăre. S-a urmărit modelarea 
arborelui rotorului cu elemente finite de grindă Timoshenko, de tip C 1 , precum 
şi, atunci când este cazul, cu elemente conice cu două noduri, de tip C 1 . 

In acest sens, după o largă consultare bibliografică (aproximativ 
200 de titluri), au fost alese pentru studiu elementele CON8 şi CON 12 care au 
fost implementate în programul DRACULA (Dinamica Rotarilor Asistată de 
CalCULAtor) în mediul MATLAB şi au fost făcute simulări numerice pentru 
validarea programului. 

1.3. Organizarea lucrării 

Lucrarea este structurată pe şase capitole, primul fiind de 
introducere, în care s-a definit problema modelării sistemelor rotor-lagăre, s-a 
prezentat scopul şi organizarea lucrării. 

In capitolul al doilea s-a descris stadiul actual al cercetărilor în 
domeniul modelării dinamice a sistemelor rotor-lagăre şi problemele derivate 
din acesta, cu peste 200 de referiri bibliografice. 

începând cu al treilea capitol s-au prezentat mai întâi elementele 
finite cilindrice de bară cu două noduri la capete, adică elementul de bară de tip 
Bernoulli-Euler şi elementul de bară de tip Timoshenko. Pentru fiecare element 
în parte s-a arătat modul de obţinere a ecuaţiilor de mişcare, a funcţiilor de 
formă în două plane perpendiculare între ele, apoi modul de obţinere a 
matricelor de masă, de rigiditate şi a vectorului forţelor de dezechilibru masic 
pentru elementul de tip Bernoulli-Euler şi pentru elementul de tip Timoshenko, 
în plus şi matricea giroscopică. Apoi s-au prezentat două tipuri de elemente 
finite tronconice de bară cu două noduri, numite CON 12 şi CON 8, primul 
având şase grade de libertate pe nod, iar al doilea patru grade de libertate pe nod. 
Ambele se bazează pe teoria de bară Timoshenko, considerând în calcule şi 
deformaţia de forfecare, dar în moduri diferite de la un element la celălalt. S-au 
determinat matricele elementelor care au dimensiuni 8x8 pentru elementul CON 
8 şi 12x12 pentru elementul CON 12, pentru care în vederea reducerii ordinului 
matricelor la 8x8, se face o condensare statică. S-a prezentat şi ecuaţia 
matriceală de mişcare pentru fiecare element conic. Tot în capitolul trei s-au 
tratat şi discurile rigide şi lagărele. S-au prezentat matricele prin care acestea se 
introduc în modelul de calcul. 

In capitolul patru s-a descris răspunsul dinamic al rotarilor modelaţi 
cu elemente de bară de tip Timoshenko cilindrice şi conice. Au fost apoi 
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realizate simulări numerice cu ajutorul programului mai sus amintit şi s-au făcut 
observaţii în urma efectelor puse în evidenţă. 

In capitolul cinci s-au prezentat rezultatele experimentale obţinute 
în urma măsurărilor realizate în laboratorul de Dinamica rotorilor din Catedra de 
Rezistenţa materialelor din Universitatea Politehnica Bucureşti. Aceste rezultate 
au fost comparate cu rezultatele obţinute pentru acelaşi rotor în urma simulării 
numerice, în vederea validării programului considerat. 

In capitolul şase s-au prezentat concluziile desprinse în urma 
studiilor efectuate şi contribuţiile originale ale autoarei. 

In încheiere s-au adăugat patru Anexe cu diverse calcule ajutătoare. 
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2. Stadiul actual al cercetărilor în domeniul modelării 
dinamice a sistemelor rotor-lagăre 

Dinamica sistemelor rotor-lagăre s-a conturat ca o disciplină de 
sine stătătoare şi un domeniu de cercetare interdisciplinară în momentul în care 
s-a recunoscut importanţa efectului lagărelor şi etanşărilor asupra răspunsului 
dinamic al rotorilor. 

Interesul crescând privind studiul vibraţiilor maşinilor cu rotor se 
datoreşte faptului că peste 80% din problemele ce apar la aceste maşini se 
datoresc vibraţiilor. Efortul continuu de a micşora greutatea maşinii pe unitatea 
de putere dezvoltată, precum şi costul extrem de ridicat al întreruperilor 
accidentale din funcţionare, la maşinile ce lucrează în instalaţii cu flux 
tehnologic continuu, subliniază importanţa prevenirii vibraţiilor la maşinile cu 
rotor. Pe măsură ce noile soluţii constructive s-au apropiat de limitele de 
rezistenţă ale materialelor, problemele de vibraţii s-au înmulţit, impunând 
dezvoltarea cercetării şi metodelor de proiectare în două domenii de prim interes 
practic: dinamica sistemelor rotor-lagăre şi vibraţiile ansamblului disc-palete. 

Creşterea spectaculoasă a puterii maşinilor nu a fost însoţită de o 
creştere corespunzătoare a dimensiunilor acestora datorită cerinţelor privind 
creşterea puterii dezvoltate pe un kilogram de material. Semnificativă pentru 
progresul tehnic în acest domeniu este dezvoltarea turbinelor cu abur în Europa. 
De la prima turbină cu acţiune construită în 1883 de inginerul suedez Gustav de 
Laval şi prima turbină cu reacţiune cu mai multe trepte, construită în 1884 de 
englezul Charles Parsons şi până la turbinele pentru centrale nucleare de astăzi, 
evoluţia este spectaculoasă. 

Dacă la maşinile mici, de exemplu motoarele electrice, cu turaţii 
relativ mici şi lagăre cu rulmenţi, echilibrarea şi calculul dinamic al rotorului nu 
ridică probleme deosebite, în schimb, maşinile mari, cu rotori lungi şi flexibili, 
cu lagăre de alunecare, etanşări, piedestaluri şi carcase relativ elastice, cu turaţii 
mari, au impus dezvoltarea şi perfecţionarea continuă a metodelor de calcul 
dinamic şi de măsurare a vibraţiilor. 

Primul studiu al turaţiilor critice ale unui arbore elastic de secţiune 
constantă a fost făcut în 1869 de Rankine [142], care a introdus termenul de 
turaţie critică, crezând că este vorba de o instabilitate, rotorul neputând 
funcţiona deasupra turaţiei critice. Deşi prima rezolvare corectă a fost dată de 
Foppl [53], confuzia a persistat până la publicarea în 1919 a articolului lui 
Jeffcott [87] asupra precesiei datorită dezechilibrului masic. 
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In 1 894 Dunkerley [46] a publicat rezultatele studiilor sale privind 
turaţiile critice ale arborilor cu mai multe discuri şi a propus metoda şi formula 
care îi poartă numele. In 1916 Stodola [164] a publicat un studiu despre 
influenţa lagărelor asupra precesiei rotorului elastic. 

Influenţa amortizării interne asupra mişcării rotorilor a fost 
considerată pentru prima oară de Newkirk [124] în 1924 în studiul vibraţiilor 
compresoarelor de la furnalele înalte. In 1925, Newkirk şi Taylor [123] 
semnalează fenomenul de instabilitate a rotorilor în lagăre cu ulei (fenomenul 
„oii whirl") explicat abia în 1952 de Poritsky [136]. 

In 1933 Smith [157] a publicat un articol de sinteză asupra 
principalelor probleme ale Dinamicii rotorilor. Intre 1932-1935 Robertson [144] 
a studiat în detaliu principalele probleme de instabilitate a mişcării rotorilor şi 
precesia acestora în regim tranzitoriu, la trecerea prin turaţiile critice. 

In 1945 Prohl [137] a propus metoda care îi poartă numele pentru 
calculul turaţiilor critice la rotorii cu mase concentrate. Lund [111, 107] a 
dezvoltat analiza cu matrice de transfer a sistemului rotor-lagăre. Ruhl [151] şi 
Nordmann [127] au utilizat în tezele de doctorat metoda elementelor finite 
pentru analiza sistemelor rotor-lagăre, dar primele articole au fost publicate de 
Ruhl şi Booker [150] în 1972 şi Gash [57] în 1973. Reducerea modelului cu 
elemente finite a fost studiată începând cu 1980 de Rouch şi Kao [148] şi Jăcker 
[86], ultimul introducând şi efectul fundaţiei asupra răspunsului rotorului. 

In 1984, Hashish şi Sankar publică un articol [20], în care un sistem 
rotor-lagăr este modelat utilizând metoda elementelor finite. Modelul matematic 
ia în considerare momentul giroscopic, inerţia la rotaţie, deformaţia de forfecare, 
amortizarea internă vâscoasă, amortizarea histeretică, rigiditatea arborelui 
amortizarea din lagăre şi elasticitatea suportului lagărelor. Se foloseşte 
elementul de grindă Timoshenko şi se prezintă modul de construire a matricelor 
de masă, de rigiditate, de amortizare şi giroscopică, de ordin 4x4, necesare 
pentru calcularea răspunsului sistemului rotor-lagăr. 

O problemă importantă care se pune în tehnica modelării arborelui 
din sistemele rotor-lagăre şi care se tratează în lucrarea de faţă pe larg este 
alegerea tipului de element de arbore utilizat în calcule, respectiv elemente de 
arbore de tip Timoshenko care ţin cont pe lângă încovoiere şi de inerţia la rotaţie 
a elementului şi de deformaţiile de forfecare ale acestuia, sau elemente de arbore 
Bernoulli-Euler care nu ţin cont decât de încovoierea elementului. De asemenea, 
se pune problema eficienţei modelării secţiunilor tronconice de arbore, prin 
elemente finite cilindrice în trepte sau prin elemente finite conice, acesta fiind şi 
obiectul lucrării. 

In ceea ce priveşte vibraţiile grinzilor Timoshenko, literatura 
tehnică este foarte bogată. Există publicaţii referitoare la influenţa secţiunii 
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transversale neuniforme a grinzilor Timoshenko asupra răspunsului dinamic al 
acestora, la influenţa maselor rigide, la influenţa deformaţiilor de forfecare, la 
controversatul „spectru de frecvenţe de ordinul doi", la stabilitatea dinamică a 
acestor grinzi ş.a. 

Cele două articole clasice publicate de Timoshenko în 1921 şi 1922 
[172, 173], cu privire la efectul deformaţiilor de forfecare asupra vibraţiilor 
transversale ale grinzilor, au avut o influenţă covârşitoare în dezvoltarea 
dinamicii structurilor în secolul al douăzecilea. 

Trei decenii mai târziu, după ce Timoshenko a obţinut faimoasa sa 
ecuaţie, Mindlin a obţinut teoria cu privire la plăcile în vibraţie [119], care 
constituiau o altă piatră de hotar în teoria dinamică a elementelor structurale. 

La aplicarea teoriei Timoshenko la grinzile care vibrează există 5 
aspecte importante. Ele pot fi extinse şi asupra elementelor de arbore de tip 
Timoshenko din alcătuirea unui rotor. Acestea sunt: 

1 . Considerentele teoretice de bază; 

2. Predicţia „spectrului de frecvenţe de ordinul doi"; 

3. Consideraţii asupra elementelor structurale cu secţiune transversală 
neuniformă; prezenţa maselor concentrate sau a rezemărilor elastice; 

4. Prezenţa fundaţiilor elastice; cazuri de stabilitate dinamică; efecte 
neliniare; 

5. Rezolvarea problemelor utilizând metoda elementelor finite. 

Cele două articole clasice ale lui Timoshenko [172, 173] au pregătit 
terenul pentru investigaţii analitice în legătură cu: 

a) determinarea pulsaţiilor proprii [97, 8, 175, 40, 68, 118, 78, 80- 
82, 33, 50, 21, 27, 167, 24, 1, 89, 159, 165, 96, 116]; 

b) analiza mişcării forţate [77, 19, 47]; 

c) studii privind factorul de forfecare [34, 88]. 

După cum se ştie, unul din fundamentele teoriei Timoshenko este 
determinarea factorului de forfecare. Acesta, înmulţit cu aria secţiunii 
transversale, permite calculul unei arii efective de forfecare, care se utilizează în 
calcule în locul ariei reale. 

Introducerea lui a fost necesară datorită aproximaţiei din teoria lui 
Timoshenko, în care se neglijează deplanarea secţiunii transversale, şi se 
calculează o tensiune tangenţială medie, constantă pe înălţimea barei, necesară 
menţinerii secţiunii plane. Un studiu exhaustiv, citat în multe lucrări actuale, 
este articolul lui Cowper [34]. In general, se acceptă o valoare statică a 
factorului de forfecare. în cazul elementelor finite de arbore tronconice, se 
adoptă o valoare medie, de obicei calculată la mijlocul elementului. In alte 
lucrări [2 1 6] se arată că factorii de forfecare ar trebui să fie funcţie de frecvenţa 
vibraţiilor. 
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Un număr considerabil de lucrări au fost prezentate pe tema 
intitulată „Predicţia spectrului de frecvenţe de ordinul doi pentru grinzile 
Timoshenko" [42, 3, 15, 102, 160]. Stephen [160] a publicat un excelent studiu 
examinând spectrul de ordinul doi din punctul de vedere al propagării undelor. 
El a realizat această analiză utilizând teoria Pochhammer-Chree (în contextul 
teoriei elasticităţii) în cazul unor bare subţiri, de secţiune circulară şi de lungime 
infinită şi a ajuns la următoarea concluzie: „Teoria exactă Pochhammer-Chree 
cu privire la propagarea unei unde într-o grindă de secţiune circulară, prezice 
posibilitatea unei multitudini de moduri de vibraţie de acelaşi ordin. Teoria de 
bară Timoshenko prezice existenţa a câte două moduri de vibraţie de acelaşi 
ordin. Spectrul de frecvenţe de ordinul doi trebuie să fie considerat, pe de o 
parte, o consecinţă inevitabilă a comportării structurilor, iar pe de altă parte 
rezultatul unei teorii inginereşti remarcabile". 

Articolul lui Lee [101] constituie una din primele cercetări privind 
grinzile Timoshenko de secţiune transversală neuniformă. Alte studii interesante 
pe aceeaşi temă sunt [45, 84]. 

Cazul rezemărilor elastice ale capetelor barelor este tratat în [4] iar 
influenţa secţiunilor transversale neuniforme asupra acestora, influenţa maselor 
rigide etc. sunt luate în considerare în lucrările [115, 72, 146, 147]. 

Grinzile Timoshenko rezemate pe fundaţii elastice sunt tratate în 
[36, 187, 189, 2, 114, 18, 203, 204, 35, 55, 103, 176]. 

Cazul grinzilor vâscoelastice este prezentat în [56, 100, 130], iar 
răspunsul plastic în [191]. 

Au fost publicate numeroase studii privind stabilitatea şi analiza 
dinamică a stabilităţii barelor Timoshenko [122, 20, 92, 93, 2, 128, 152, 155, 18, 
131, 156, 176]. 

Vibraţiile barelor şi cadrelor Timoshenko sunt prezentate în [22, 23, 
28,188,41,79,200,51]. 

Cazul grinzilor Timoshenko pretensionate este tratat de Sun [166], 
cazul vibraţiilor neliniare este prezentat de Rao, Raju şi Raju [143], iar situaţiile 
termoelastice sunt considerate de Massalas şi Kalpakidis [113]. 

Grinzile Timoshenko pentru structurile de poduri sunt studiate în 
190, 76], iar barele curbe în [145, 85]. 

In literatură sunt prezentate o mulţime de elemente finite de bară 
Timoshenko [11, 17, 90, 25, 154, 6, 37, 125, 169, 168, 170, 171, 38, 71, 174, 95, 
112, 153, 39, 44, 5, 29]. Thomas, Wilson şi Wilson [169] au împărţit aceste 
elemente în două clase: elemente simple şi elemente complexe. Un element 
simplu este un element unidimensional, solicitat la încovoiere într-un singur 
plan şi care are în total 4 grade de libertate, câte două în fiecare dintre cele două 
noduri; un element complex este un element care are mai mult de patru grade de 
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libertate, adică are mai mult de două grade de libertate în fiecare din cele două 
noduri, sau are mai mult de două noduri. 

Articolul lui Akella şi Craggs [5] are aplicaţii la barele în rotaţie. 
Studiul lui Davis, Henshell şi Warburton [37] extinde aplicarea elementelor 
Timoshenko la studiul vibraţiilor. 

Dezvoltarea metodei elementelor finite pentru analiză în dinamica 
rotorilor este subiectul a foarte multe publicaţii recente [10, 30, 49, 58, 61, 67, 
70, 71, 75, 83, 90, 103, 117, 120,121, 129, 133-135, 139-141, 148, 149-151, 
153, 161, 162, 168-171, 174, 177, 179-181, 186, 205]. Acestea au inclus în 
studiul efectuat cu diferitele tipuri de elemente finite efectele inerţiei la rotaţie, 
efectul momentului giroscopic, al încărcării axiale, al amortizării interne şi al 
deformaţiei de forfecare. Totuşi, în majoritatea cazurilor studiate, geometria 
elementului finit utilizat a fost limitată la elementul cilindric. 

In ultimii ani, mulţi cercetători au extins activitatea de analiză în 
dinamica rotorilor utilizând metoda elementelor finite, după cum urmează: 

Nelson şi McVaugh [120], au prezentat elementul de grindă 
cilindric de tip Rayleigh care ţine cont de încovoiere, de efectul inerţiei la rotaţie 
şi la translaţie, de momentele giroscopice şi de încărcarea axială, dar neglijează 
forfecarea. Acest element a fost îmbunătăţit de Zorzi şi Nelson [205], care au 
introdus efectele amortizării vâscoase interne şi a amortizării histeretice. Mai 
târziu, Nelson [121] a adăugat la acest tip de element deformaţia de forfecare, ce 
a condus la elementul de grindă tip Timoshenko care a fost apoi perfecţionat de 

rr rr 

Ozgiiven şi Ozkan [129] luând în calcul şi amortizarea internă după modelul lui 
Zorzi şi Nelson. Toate aceste elemente erau cilindrice, având aria secţiunii 
transversale şi momentele de inerţie constante în lungul lor. Ecuaţia de mişcare a 
fost determinată utilizând 8 grade de libertate pe: câte două translaţii şi două 
rotaţii la capetele elementului, în plane perpendiculare între ele. 

Maşinile moderne au rotori cu geometrii simplificate în modelele 
utilizate pentru calculul lor. De obicei rotorul se modelează cu o serie de 
elemente cilindrice în trepte. Dar, pentru arborii conici, erorile introduse de 
această aproximare pot fi foarte mari. Ca urmare, Rouch şi Kao [149] au 
elaborat un element finit de grindă de tip Timoshenko utilizabil la modelarea 
rotorilor în trepte. Acest element se bazează pe concluziile lui Thomas, Wilson 
şi Wilson [169] potrivit cărora cea mai bună metodă pentru luarea în calcul a 
deformaţiei de forfecare este de a adăuga două grade de libertate adiţionale la 
fiecare capăt de element, obţinând astfel 12 grade de libertate pe element. 

Rouch şi Kao, pornind de la aceste observaţii, au inclus efectele 
giroscopice şi au exprimat aria şi momentele de inerţie ale secţiunii elementului 
ca polinoame de gradul 2, respectiv 4 în funcţie de raza secţiunii transversale. 
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Modelarea porţiunilor tronconice ale unui rotor prin elemente 
cilindrice în trepte nu este satisfăcătoare pentru că ea afectează atât simplitatea 
cât şi exactitatea modelului rotorului. Este preferabilă modelarea unui rotor 
conic, direct cu elemente conice, atât din punctul de vedere al comodităţii, cât şi 
din punct de vedere al preciziei calculelor. 

In 1985, Greenhill, Bickford şi Nelson [70] completează lucrările 
existente prezentând un element conic care ia în considerare toate efectele de 
mai sus. 

Un alt element finit conic utilizat în analiza liniară în dinamica 
rotorilor a fost elaborat de Genta, Gugliotta şi Brussino în 1984 [60]. In 
continuare, Genta şi Gugliota prezintă în 1988 [58] acest tip de element conic în 
detaliu, afirmând că rotorii în general şi turbinele cu gaz în particular, includ 
adesea părţi care nu sunt uşor de modelat cu elemente cilindrice, decât dacă se 
foloseşte o discretizare foarte fină, ceea ce duce la un număr foarte mare de 
grade de libertate pentru întreg modelul, deci la calcule foarte laborioase. 

Utilizând în calcule modelele prezentate în [58, 60, 70] se pot 
obţine următoarele facilităţi: analiză statică, calculul turaţiilor critice, 
determinarea diagramei Campbell, calculul stabilităţii şi răspunsul la 
dezechilibru (amortizat şi neamortizat). 

Elementul conic descris de Greenhill în [70], bazat pe teoria de bară 
Timoshenko prezintă, în fiecare plan de încovoiere, în plus faţă de gradele de 
libertate obişnuite - cele două deplasări şi două rotiri- un grad de libertate 
suplimentar numit deformaţie de forfecare. 

In ceea ce priveşte comportarea discurilor montate pe arborele unui 
rotor, încă din 1948 Green [69] a studiat efectul giroscopic al unui disc rigid 
asupra precesiei unui rotor flexibil în consolă, iar în 1957 Downham [43] a 
confirmat experimental existenţa precesiei inverse. 

In 1916 Stodola [164] a publicat un studiu asupra influenţei 
lagărelor asupra precesiei rotorului elastic. 

In 1925, Newkirk şi Taylor [123] semnalează fenomenul de 
instabilitate a rotorilor în lagăre cu ulei (oii whirl) explicat abia în 1952 de 
Poritsky [136]. 

între 1955 şi 1965, Hagg şi Sankey [74], Sternlicht [163], Lund 
[105] şi alţii au dezvoltat teoria lagărelor hidrodinamice, Yamamoto [202] a 
studiat lagărele de rostogolire. 

Intre 1963-1967, Lund [110] şi Glienicke [66] au publicat valori ale 
coeficienţilor dinamici elastici şi de amortizare pentru mai multe tipuri de 
lagăre. 

Problema lagărelor este prezentată şi în lucrările [12, 32, 52, 65, 94, 
98, 108, 109, 132, 192]. 
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In prezent, arborii confecţionaţi din materiale compozite sunt 
folosiţi în special ca arbori de transmisie în domeniul subcritic [195]. 

Greutatea unui sistem de arbori de transmisie, determinată de 
lungimea sa şi de puterea transmisă, poate fi optimizată folosind materialele 
compozite şi utilizând rotorul în domeniul supracritic, adică la o turaţie mai 
mare decât prima turaţie critică. Cea mai mare turaţie, este, pentru un sistem 
particular cea care transmite cel mai mic cuplu în scopul furnizării puterii. 
Există mai multe publicaţii unde au fost utilizate tehnici de optimizare prin 
metode numerice pentru a duce la o proiectare foarte bună din punct de vedere al 
greutăţii arborelui. 

Wettergren [194] a prezentat comportarea dinamică a rotorilor 
confecţionaţi din materiale compozite arătând că este diferită de cea a rotorilor 
confecţionaţi din materiale convenţionale datorită imperfecţiunilor de fabricaţie 
a materialelor compozite. 

O altă diferenţă importantă între rotorii confecţionaţi din materiale 
compozite şi cei confecţionaţi din oţel este amortizarea pe care o introduce 
materialul. Aceasta poate fi până la de două ori mai mare pentru arborii din 
materiale compozite, în special atunci când unghiul dintre fibrele materialului 
rotorului este mai mic de 45°. Pentru arborii cu o turaţie mică, acesta este un 
mare avantaj. Există multe cazuri când torsiunea arborelui a fost evitată folosind 
materiale compozite pentru realizarea acestuia în loc de oţel. Din nefericire, 
amortizarea materialului este un dezavantaj pentru rotorii utilizaţi în domeniul 
supracritic. 

Wettergren şi Olson [192] au studiat vibraţiile laterale ale unui 
rotor orizontal, elastic, rezemat pe lagăre elastice. 

Wettergren şi Olson au observat că instabilitatea rotorului apare, în 
general, în apropierea rezonanţei de dezechilibru. Această instabilitate poate fi 
evitată realizând diverse combinaţii între amortizarea materialului, amortizarea 
din lagăre şi dispunerea în diferite poziţii ale discului între lagăre. Rezonanţa, 
care se datoreşte greutăţii, care apare în apropierea unei turaţii egale cu jumătate 
din prima turaţie critică, poate fi redusă simţitor cu ajutorul amortizării interne 
vâscoase. Amortizarea internă este, totuşi, rareori vâscoasă. Pentru un rotor cu 
amortizare histeretică, în ceea ce priveşte amortizarea materialului, analiza 
devine foarte complicată. Wettergren [196] a elaborat o metodă pentru analiza 
amortizării histeretice în cazul în care orbita mişcării nu este circulară. 

Wettergren [193] a arătat că se iau în calcul în analiza arborilor din 
materiale compozite, rezemaţi pe lagăre hidrodinamice, efectele combinate ale 
peliculei de lubrifiant din lagăre şi ale amortizării materialului. Cea mai mare 
turaţie permisă pentru un rotor scade din cauza amortizării interne. 
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Surse de amortizare internă în materialele compozite sunt 
fenomenele microplastice sau vâscoelastice asociate matricei materialului şi 
lunecările relative dintre interfeţele matricei şi armăturii. Amortizarea internă a 
compozitelor poate fi influenţată de proprietăţile matricei şi ale armăturii, de 
fracţia de volum a fibrei, de mărimea golurilor, de orientarea fibrelor armăturii 
în raport cu direcţiile de solicitare şi de tratarea suprafeţelor fibrelor. Dacă 
turaţia este mai mare decât aproximativ 1000 rad/sec, factorul de pierdere începe 
să crească cu creşterea turaţiei. 

Willway [198] ajunge la concluzia că în răşinile epoxidice există o 
variaţie considerabilă în disiparea energiei. Este evident că proprietăţile 
dinamice ale răşinilor vâscoelastice depind de turaţia rotorului şi de temperatura 
de lucru a acestuia. 

Wettergren [195] tratează amortizarea ca fiind proporţională cu 
modulul de elasticitate şi prezintă o metodă de utilizare a amortizării 
materialului la rotorii din compozite şi analizează şi alte diferenţe ce există între 
rotorii confecţionaţi din materiale obişnuite şi cei confecţionaţi din compozite în 
sensul amortizării materialului. Wettergren studiază amortizarea materialului 
pentru o răşină epoxy cu fibre de carbon având un modul de elasticitate ridicat, 
pentru diferite unghiuri de dispunere a fibrelor în raport cu direcţia tensiunilor. 
Au fost luate în calcul situaţiile în care fibrele sunt în lungul tensiunilor, 
perpendiculare pe acestea şi tangenţiale. S-a ajuns la concluzia că amortizarea, 
deci coeficientul de amortizare are un minim când unghiul de dispunere al 
fibrelor este cuprins între 10° şi 45°, care din fericire, pentru majoritatea 
aplicaţiilor este intervalul care asigură obţinerea celei mai mici greutăţi posibile 
a rotorului. 

Mai mult, materialele compozite pot avea amortizarea mai mare 
decât a oţelului sau a aluminiului, o dată până la de două ori. 

Studiul a fost făcut utilizând elemente finite de grindă specifice 
laminatelor compozite cu pereţi subţiri, având câte 10 grade de libertate pe nod, 
element dezvoltat de Wu şi Sun [199]. 

Problema comportării materialelor compozite în dinamică este 
tratată şi în [7, 64, 178]. 

Studii cuprinzătoare privind dinamica rotorilor sunt prezentate şi în 
[9,16,59,63,91,99,106,126,182-185,201]. 

Instrumentul matematic utilizat în calculul dinamic al rotorilor este 
prezentat în [13, 14, 26, 31, 54, 62, 138, 197]. 
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3. Modele cu elemente finite ale sistemelor rotor-lagăre 



Proiectarea unor sisteme rotor-lagăre de performanţe ridicate 
necesită modele analitice şi experimentale care să simuleze cât mai bine 
comportarea sistemelor reale. Totodată se urmăreşte şi obţinerea de modele 
eficiente din punct de vedere al timpului şi costului calculelor, pentru că nu 
trebuie uitat că modelele sunt o reflectare a stadiului dezvoltării tehnologiilor de 
calcul. 

La proiectarea mecanică a rotorului, lagărelor şi a structurii de 
suport, se ţine cont că acestea lucrează ca un sistem, răspunzând împreună 
solicitărilor dinamice, intercondiţionându-se reciproc. Rotorul este o parte 
integrantă a unui sistem dinamic, cunoaşterea lui fiind determinată de poziţia şi 
rigiditatea lagărelor, a etanşărilor, a piedestalurilor şi a fundaţiei, precum şi de 
amortizarea acestora. 

Rotorul este elementul principal al unei maşini rotative, funcţia lui 
fiind generarea sau transmiterea puterii. El constă dintr-un arbore pe care se 
montează discuri cu palete sau roţi cu canale radiale, iar în cazul maşinilor 
electrice - înfăşurările bobinelor. Rotorul nu este niciodată perfect rigid. Totuşi, 
în practică, se numesc rotori rigizi cei care lucrează la 1/3 din prima turaţie 
critică de încovoiere. Rotorii elastici sunt cei care lucrează aproape de, sau 
deasupra primei turaţii critice de încovoiere, astfel că forţele centrifuge datorate 
dezechilibrului remanent influenţează deformaţiile rotorului. 

La cele mai multe maşini, rotorii au arbori cu secţiune transversală 
axial-simetrică. Dacă pe anumite porţiuni secţiunea transversală este 
nesimetrică, atunci rigiditatea la încovoiere a rotorului faţă de o axă transversală 
fixă variază în timpul rotirii şi apar mişcări de precesie nesincronă şi 
instabilităţi. Arborele rotorilor poate fi modelat ca bara de tip Bernoulli-Euler 
sau Timoshenko, ţinând cont sau nu de efectul inerţiei la rotaţie faţă de axa 
transversală, de efectul forfecării, incluzând în plus efectul cuplurilor 
giroscopice. Discurile, considerate de obicei rigide, intervin prin parametrii 
concentraţi: masa şi momentele de inerţie masice, polar şi transversal. Lagărele 
se aleg în funcţie de încărcarea pe reazem şi de turaţie, ţinând cont de solicitarea 
dinamică, de spaţiul disponibil, de pierderile de putere, de simplitatea soluţiei 
constructive şi de cerinţele de durabilitate şi fiabilitate. La început, lagărele au 
fost considerate rigide. Ulterior s-a ţinut cont şi de elasticitatea radială a 
lagărelor, apoi de amortizarea acestora. La lagărele cu rulmenţi şi la lagărele cu 
aer, amortizarea poate fi neglijată. La lagărele de alunecare, proprietăţile elastice 
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şi de amortizare depind de turaţie şi de încărcare. La lagărele cu rulmenţi, 
constanta elastică este independentă de turaţia maşinii şi de încărcare. In general 
se consideră numai rigiditatea radială a lagărelor, cea unghiulară fiind relativ 
mică. 

La lagărele de alunecare, în cazul regimului dinamic staţionar, forţa 
datorită presiunii totale egalează sarcina statică pe lagăr. Dacă fusul în mişcare 
are o mişcare de precesie, în filmul de lubrifiant apar tensiuni suplimentare, care 
produc forţe dinamice suplimentare asupra fusului. Forţa dinamică depinde de 
amplitudinea deplasării şi de viteza mişcării centrului fusului, însă, spre 
deosebire de forţele din alte lagăre, nu are aceeaşi direcţie cu mişcarea pe care o 
produce, fiind defazată în timp şi în spaţiu. 

Descompunând forţa dinamică în două componente în lungul unor 
axe de coordonate fixe Oy şi Oz legate de lagăr, notând cu y, respectiv cu z, 
componentele deplasării centrului fusului pe cele două direcţii, se stabilesc 
relaţiile: 

Fz=-K^y-C^y-K K z-C zz z. 

Aceste relaţii sunt exacte numai pentru deplasări mici, însă în 
practică s-au dovedit valabile pentru deplasări de trei ori mai mari decât jocul în 
lagăr. Coeficienţii elastici K yy , K yz , K zy , K zz şi coeficienţii de amortizare C yy , C yz , 
C zy , C zz se calculează cu ajutorul relaţiilor stabilite în teoria lubrificaţiei şi sunt 
caracteristici unui anumit lagăr. Ei depind de configuraţia lagărului, de 
proprietăţile lubrifiantului, fiind funcţie de poziţia centrului fusului, variind deci 
cu turaţia maşinii. Neegalitatea termenilor de cuplaj K yz ^K zy este cauza unui 
anumit tip de precesie instabilă a rotorilor. 

Etanşările inelare cu fluid sau cu gaz se consideră de obicei a fi 
izotrope, deci coeficienţii diagonali ai matricelor de rigiditate şi de amortizare 
sunt egali, iar coeficienţii de cuplaj sunt egali şi de sens contrar. Cele două 
componente ale forţei cu care etanşările lucrează asupra rotorului se exprimă sub 
forma: 
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Termenul inerţial este neglijabil la etanşările cu gaz, la care 
rigidităţile directe pot fi foarte mici, chiar negative. In cazul etanşărilor relativ 
lungi se introduc şi coeficienţi dinamici unghiulari, deoarece asupra rotorului 
acţionează şi cupluri, iar forţele produc şi rotiri, iar cuplurile produc şi deplasări 
liniare. 

Piedestalurile deformabile intervin în răspunsul dinamic al 
maşinilor, în special în cazul ventilatoarelor, suflantelor, pompelor centrifuge şi 
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al turbomaşinilor cu carcase flexibile şi lagăre în consolă. Modelul de calcul 
include de obicei rigiditatea şi amortizarea suporţilor lagărelor, iar în unele 
cazuri şi masa echivalentă a acestora. 

Fundaţia, placa de bază şi terenul de rezemare sunt incluse mai rar 
în modelul de calcul, influenţa lor asupra răspunsului rotorului fiind în general 
mică. Totuşi, în unele cazuri, mai ales la ventilatoare mari, pe piedestaluri de 
beton, se ţine cont de elasticitatea terenului de fundaţie. 

De obicei, la lagăre şi piedestaluri, chiar în cazul caracteristicilor 
independente de turaţie, rigiditatea în plan orizontal este mai mică decât cea în 
plan vertical. Această anizotropie duce la o dedublare a turaţiilor critice. Uneori 
însă, datorită amortizării puternice, separarea celor două turaţii critice din cauza 
ortotropiei lagărelor nu apare distinct în răspunsul la dezechilibru al rotorilor. De 
asemenea, este posibil ca unele moduri de precesie, în special cele retrograde, să 
fie amortizate supracritic, deci să nu apară în diagramele frecvenţelor proprii şi 
nici în cele de răspuns la dezechilibru. 

Metodele de calcul şi modul de interpretare a rezultatelor analizei 
dinamice a rotorilor au avut o evoluţie remarcabilă. In urmă cu 40 de ani, 
calculele se făceau de preferinţă grafic, considerând reazeme rigide în dreptul 
lagărelor şi tratând separat fiecare rotor dintr-o linie de arbori. Analiza se limita 
la determinarea turaţiilor critice neamortizate şi proiectarea rotorului astfel încât 
între domeniul turaţiilor de lucru şi orice turaţie critică să existe oarecare 
separare. 

Ulterior, reazemele rigide au fost înlocuite cu elemente deformabile 
având constanta elastică egală cu cea a filmului de lubrifiant din lagăre, 
adăugându-se apoi şi efectul elasticităţii piedestalurilor. S-a luat în considerare 
efectul cuplării mai multor rotori dintr-o linie de arbori şi s-a apelat la 
calculatoare numerice de mare capacitate pentru analiza rotorilor cu mai multe 
discuri. 

In unele cazuri au apărut discrepanţe mari între calcule şi încercări, 
care au fost micşorate introducând efectul amortizării şi determinând turaţiile 
critice amortizate. Calculând răspunsul rotorului la dezechilibru, s-a constatat că 
turaţiile la care orbitele precesiei sincrone au raze vectoare de valori maxime, 
numite turaţii critice de dezechilibru, diferă de turaţiile critice neamortizate, 
fiind mai aproape de cele amortizate. Astfel s-a observat că nu toate turaţiile 
potenţial critice sunt într-adevăr critice, amortizarea puternică aplatizând 
diagrama răspunsului la dezechilibru iar trecerea prin turaţia critică facându-se 
fără o creştere a amplitudinii mişcării de precesie. 

In domeniul compresoarelor şi turbinelor pentru instalaţii 
tehnologice şi pentru centrale electrice, tendinţa actuală de creştere a gabaritului 
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maşinilor şi a turaţiilor de lucru a dus la o nouă generaţie de maşini la care, 
inevitabil, una sau două turaţii critice sunt situate în domeniul turaţiilor de lucru. 

Pe măsură ce maşinile devin mai mari, elasticitatea filmului de 
lubrifiant din lagăre şi a suporţilor acestora joacă un rol tot mai important în 
comparaţie cu rigiditatea rotorului, ducând la scăderea turaţiilor critice şi 
interferarea lor cu domeniul turaţiilor de lucru. Deoarece caracteristicile 
dinamice ale lagărelor nu sunt cunoscute exact, turaţiile critice nu pot fi 
determinate precis, deci criteriul tradiţional de proiectare urmărind evitarea 
funcţionării la, sau aproape de, o turaţie critică nu mai poate oferi siguranţa 
necesară. 

In practică s-a constatat că se poate funcţiona perfect de sigur şi 
fiabil la turaţii critice bine amortizate dacă nivelele vibraţiilor nu depăşesc 
limitele admisibile şi dacă rotorul nu are o sensibilitate pronunţată la 
dezechilibre masice. Rezultă că răspunsul rotorului la dezechilibru poate furniza 
informaţiile cele mai utile privind viabilitatea unei soluţii constructive. 
Efectuând acest calcul pentru diferite distribuţii ale dezechilibrului, judicios 
alese, astfel încât să se accentueze deformaţia rotorului la diferite turaţii critice 
de dezechilibru, se poate stabili „cât de critică" este fiecare dintre aceste turaţii 
şi ce măsuri trebuie luate pentru ca amplitudinile vibraţiilor să rămână în limite 
normale, chiar în prezenţa dezechilibrelor ce apar în timpul funcţionării normale 
(eroziuni, depuneri, ruperi de componente, deformări termice etc). 

3.1. Elemente finite de arbore în rotaţie 

La modelarea cu elemente finite a sistemului rotor-lagăre, arborele 
este discretizat utilizând un număr mai mic de elemente decât numărul de 
secţiuni care ar fi necesar pentru rezolvarea problemei cu metoda matricelor de 
transfer. Acesta este unul dintre avantajele metodei elementelor finite faţă de 
metoda matricelor de transfer. Modelarea cu elemente finite conduce la scrierea 
sistemului de ecuaţii dinamice generale. De asemenea, în modelarea cu elemente 
finite, pot fi uşor incorporate alte efecte, cum ar fi elasticitatea suportului 
lagărelor, amortizarea din lagăre, deformaţiile de forfecare, inerţia la rotaţie, 
efectul încărcării axiale, efectul amortizării interne vâscoase şi amortizarea 
histeretică. 

La baza discretizării în elemente finite stau câteva reguli care 
folosesc la modelarea unui rotor de secţiune circulară, caracterizat prin elemente 
nedeformabile, fără lungime (masice) şi elemente deformabile (tronsoane de 
arbore): 

• Rotorul se discretizează, de obicei, prin elemente de tip grindă; 
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• Fiecare element se leagă prin nodurile de la capete de elementele vecine; 

• Elementele nedeformabile (masice) se concentrează în unul sau mai multe 
noduri corespunzătoare poziţiei elementului respectiv în raport cu unul 
din capetele rotorului, masele mai importante netrebuind să fie deplasate; 

• în zonele cu variaţie mare de diametru, lungimea elementelor trebuie să 
fie relativ mică (discretizare fină); 

• Pentru un arbore cu variaţii mici de diametru în lungul său, discretizarea 
este bine să fie cu atât mai fină cu cât lagărul este mai rigid, adică 
lungimea elementelor să fie cât mai mică în vecinătatea lagărelor rigide şi 
mai mare în vecinătatea lagărelor elastice; 

• Discretizarea trebuie să fie omogenă. Cu alte cuvinte, dacă situaţia reală o 
permite, este bine, pe cât posibil, ca matricele de inerţie şi de rigiditate ale 
elementelor vecine să aibă elemente corespunzătoare de valori apropiate 
(de exemplu, un rotor de diametru constant este bine să fie discretizat în 
elemente de aceeaşi lungime); 

• Concentrarea unui disc în unul sau mai multe noduri se face în funcţie de 
modul de asamblare al discului pe arbore şi de grosimea discului raportată 
la diametrul arborelui, astfel: 

- dacă discul este montat pe arbore prin intermediul unei pene sau 
prin caneluri (blocat cu inele de siguranţă), se poate considera că el 
nu modifică rigiditatea arborelui, efectele lui fiind concentrate în 
nodul corespunzător de pe arbore; 

- dacă discul este fixat pe arbore prin fretaj la cald, atunci se poate 
considera că discul a mărit rigiditatea arborelui, raza acestuia pe 
lungimea h (h=grosimea discului în zona de contact disc-arbore) 
fiind mărită cu h/2, efectul discului concentrându-se în nodul de la 
mijlocul grosimii acestuia; 

- dacă grosimea h a discului este mare, efectele discului pot fi 
concentrate în trei noduri distribuite pe grosimea acestuia. 

Elementele finite de arbore cele mai des utilizate sunt elementele de bară 
cilindrice. In lucrarea de faţă s-au tratat îndeosebi elementele de bară tronconice. 
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3.1.1. Elementul cilindric de bară cu două noduri 

Se consideră un segment de arbore de secţiune constantă, cu două 
noduri la capete, modelat printr-un element de bară Timoshenko (Fig.3.1). 




Fig.3.1 

Pentru acest element se folosesc notaţiile: 

• L - lungimea elementului de arbore; 

• EI - modulul de rigiditate la încovoiere; 

• GA S - modulul de rigiditate la forfecare efectiv; 

• As - aria de forfecare redusă; 

• % - factorul de forfecare; 

• \i - masa pe unitatea de lungime, |u = p A; 

• (i - masa de rotaţie pe unitatea de lungime, p. = pi ; 

• p - densitatea materialului arborelui; 

• E - modulul de elasticitate longitudinal; 
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• G - modulul de elasticitate transversal; 

• I - momentul de inerţie axial al secţiunii transversale. 

Se definesc doi subvectori ai deplasărilor nodale de ordin 4x1, ce cuprind câte 
două deplasări liniare şi câte două rotiri în planele X-Y şi X-Z: 





Vi 




Wj 




Vi 

v j 
Ti. 




w j ' 
-<Pi. 



(3.3) 



Se definesc, de asemenea, vectorii corespunzători forţelor nodale: 





T 

y 












> 











(3.4) 



De remarcat semnul minus al rotirilor şi momentelor faţă de axa y. 

Regula de semn folosită aici este aceea că forţele interioare sunt 
pozitive dacă acţionează în sensul pozitiv al axelor de coordonate într-o secţiune 
a barei având normala (exterioară) pozitivă, respectiv, forţele interioare sunt 
pozitive dacă acţionează în sensul negativ al axelor de coordonate într-o secţiune 
a barei având normala (exterioară) negativă. 
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3.1.1.1. Elementul Bernoulli-Euler 

Elementul cilindric de arbore Bernoulli-Euler este un element finit 
cu două noduri la capete, având fiecare câte patru grade de libertate: două 
translaţii şi două rotaţii. In teoria de bară Bernoulli-Euler se neglijează 
deformaţiile datorită forfecării şi inerţia la rotaţie. Se ţine cont de încovoierea 
barei. 

Ecuaţia de mişcare în planul x-v 

Fie un element de bară de lungime dx, situat în planul x-y, izolat 
dintr-un arbore (Fig.3.2). 

y 4 











i 


Mz+ Mzdx 

'-^ 




T 


1 


}i v dx 

| dx 


Ty- 


■Ty' dx 



Fig.3.2 
a. Ecuaţiile de echilibru 



(£M i(2) =0) M 2 -T y .dx-n**.Ş-(M,+M,'dx) = 0 



feX - o) 

unde 



T y -|T y +T y dx ] + |uVdx=0 



V = pA ; v = 

™ ' ^ M z 

dx 

y ax 



dt 2 



(3.5) 



(3.6) 



(3.7) 



(3.8) 
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iar după simplificări rezultă: 



M z +T y =0 
-T v '+^iv = 0 



(3.9) 



b. Ecuaţia deformatei elastice 



(3.10) 



unde reprezintă inversul razei de curbură a fibrei medii deformate, adică 

rotirea specifică, sau, cu alte cuvinte, unghiul cu care se rotesc, una faţă de alta, 
în jurul axei z, două secţiuni situate la o distanţă egală cu unitatea, adică: 



(3.11) 



c. Relaţii cinematice 

Ţinând cont de figura 3.3 se poate scrie: 




Fig.3.3 



(3.12) 
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sau 



0 = \|/-v' 



(3.13) 



d. Ecuaţia de mişcare 

Introducând (3.13) în (3.10) se obţine: 

M z =m z y=my (3.14) 

deci 

M z =EI z v'' (3.15) 
Derivând prima ecuaţie din relaţia (3.9) în raport cu x se obţine: 

M z "+T y '=0 (3.16) 
dar tot din relaţia (3.9) rezultă: 

T y '=nv (3.17) 

adică 

M z "+nv = 0 (3.18) 
Derivând (3.15) de două ori în raport cu x se obţine: 

M z "=EI z v IV (3.19) 

Introducând (3.19) în (3.18) se obţine ecuaţia de mişcare pentru elementul 
Bernoulli-Euler în planul x-y: 

EI z v IV + ( iV = 0 (3.20) 

Un calcul asemănător pentru stabilirea ecuaţiei de mişcare se face şi 
în planul x-z, de unde rezultă următoarele: 
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Ecuaţia de mişcare în planul x-z 



Fie un element de bară de lungime dx, situat în planul x-z, izolat 
dintr-un arbore (Fig.3.4). 



My Tz 



. >j w dx 



Tz+ Tz' dx 




iz 



My+My* dx 



Fig.3.4 
a. Ecuaţiile de echilibru 



(X M v(2)=<>) M y + T Z dx+^wdxy-^M y +M v 'dxJ = 0 



(£F Z = O) T z + ^iwdx - (t z + T z 'dx) = 0 



unde 



A .. d 2 w 
^ = pA; w = — 



(3.21) 



(3.22) 



M. = y 



ax 



(3.23) 



' dT 

J _ z 



ax 



iar după simplificări rezultă: 



M y -T z =0 
T z - ua v - 0 



(3.24) 



(3.25) 
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b. Ecuaţia deformatei elastice 



M y =EI y K*z 



(3.26) 



unde reprezintă rotirea specifică, adică, unghiul cu care se rotesc, una faţă 

de alta, în jurul axei y, două secţiuni situate la o distanţă egală cu unitatea, 
adică: 



Nxz = — = <P 
dx 



(3.27) 



c. Relaţii cinematice 



Ţinând cont de figura 3.5 se poate scrie: 



= -w 



(3.28) 



sau 



Xz=<P' 

0 = m + w' 



(3.29) 



d. Ecuaţia de mişcare 



Introducând (3.29) în (3.26) se obţine: 



M y = EI > (p' = -EI % w" 



(3.30) 



deci 



M y =-e yW " 



(3.31) 
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Derivând prima ecuaţie din relaţia (3.25) în raport cu x se obţine: 



M y "-T z '=0 



(3.32) 




Fig.3.5 



dar tot din (3.25) 



(3.33) 



adică 



M - jiw = 0 



(3.34) 



Derivând (3.31) de două ori în raport cu x se obţine: 



M y = -EI y w 



IV 



(3.35) 



Introducând (3.35) în (3.34) se obţine ecuaţia de mişcare pentru elementul 
Bernoulli-Euler în planul x-z: 



EI w IV +|iw = 0 



(3.36) 
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Funcţii de formă pentru elementul Bernoulli-Euler în planul x-y 

Pentru aproximarea deplasărilor şi rotirilor elementului, se 

x 

introduce coordonata adimensională £ de forma \ = — şi se folosesc funcţii 
polinomiale de gradul trei de forma: 



v(0 = A^ 3 +A 2 £ 2 +A^ + A 



(3.37) 



Dar din (3.12) rezultă: 



, dv dv d£ 1 dv 

\i/ = v= — = = 

dx dx Ld£ 



Deci 



Dar 



1 dv ti ^ 

L dc, 



— = 3A 3 £ 2 + 2A 2 £ + A t (3.39) 
dx 



Rezultă că: 

^(3A 3 ^ 2 +2A 2 ^ + A 1 )=B 3 ^ 3 +B 2 ^ 2 +B^ + B 0 (3.40) 
Ordonând în membrul drept şi în membrul stâng după puterile lui \ , rezultă: 

^1^ + ^1Ş + ^. = B 3 ^ + B 2 § 2 + B 1 § + B 0 (3.41) 
iar prin identificare se obţine: 
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B 3 =0 



3 A 



2A ; 

L 

Ai 



= B 



= B 1 



= B, 



B 3 =0 

B 2 =-A 3 
2 L 3 

B,=^A 2 



astfel, pentru că B 3 = 0 , funcţia \|/(%) este de gradul doi: 

V(0 = B 2^ 2 +B^ + B 0 = 

= — A 3 £ 2 + — A 2 E + — A,. 
L L L 



(3.42) 



(3.43) 



Deplasarea liniară v şi rotirea \\f se aproximează în funcţie de 
deplasările nodale, cu ajutorul funcţiilor de formă statice N, în felul următor: 



v(0 = N 1 fe)v i +N 2 (0yi +N 3 fe)v j +N 4 (0i|/j =|N_ 
yfe) = N 1 ft)v i +N a fe)qr.+N3fe)v J + N 4 fe)|rj 

unde s-a folosit notaţia (3.3) şi în plus 



u. 



u. 



LnJ=Ln 1 
[n]=[n 1 



N 2 N 3 N 
N, N, N 



4. 



4. 



(3.44) 



(3.45) 



Calculul funcţiilor de formă 

Funcţiile de formă au următoarele proprietăţi importante: 

1 . Suma funcţiilor de formă are valoarea egală cu 1 în orice nod; 

2. Corespunzător coordonatei n, funcţia de formă N n =l, iar celelalte 
funcţii de formă au valori nule. 

Pentru \ { = v(x = o)= vfe = 0) -> N, (^) = 1;N, (^)= 1; 

Pentru Vl = V|/(x = 0) = = o) -+N 2 (^=1;N 2 O0 = 1; 
Pentru Vj = v(x = L)= v(^ = l) ->N 3 ($)=1;N 3 ($)=1; 

Pentru ^ = \|/(x = L) = \|/(^ = l) -> N 4 fe) = l;N 4 fe)= L 

Calculul funcţiilor de formă este prezentat pe larg în Anexa 1. 
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Funcţii de formă pentru elementul Bernoulli-Euler în planul x-z 

Pentru aproximarea deplasărilor şi rotirilor elementului de bară în 
acest plan, se folosesc tot funcţii polinomiale de gradul trei de forma: 



w(0=A3^ + A 2 ^ 2 + A^ + A 0 
l-cpfe) = B 3^+B 2 ^ 2 +B^ + B ( 



(3.46) 



Din relaţia (3.28) se poate scrie: 



dw 
dx 



= -cp 



şi folosind în continuare coordonata adimensională \ este valabilă relaţia: 



1 dw 
L d£ 



(3.47) 



Dar 



dw i 

— = 3A 3 ^ 2 +2A 2 4 + A 1 
dx 



(3.48) 



Rezultă că: 



^(3A 3 ^ 2 +2A 2 ^ + A 1 )=B 3 ^ 3 +B 2 ^ 2 +B^ + B 0 (3.49) 



Ordonând şi în membrul drept şi în membrul stâng după puterile lui \ , rezultă: 



3A^ 2 + A^ = B ^3 +B g +B ^ + Bo 



(3.50) 



iar prin identificare se obţine: 
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< 



B 3 =0 
2 L 





(3.51) 



astfel, pentru că B 3 = 0 , şi funcţia [ -cp(%) ] este de gradul doi: 



-9fe) = B 2^ 2 +B^ + B 0 = 




(3.52) 



Deplasarea liniară w şi rotirea cp se aproximează în funcţie de 
deplasările nodale, cu ajutorul funcţiilor de formă statice N, în felul următor: 



Folosind proprietăţile funcţiilor de formă, se calculează expresiile 
acestora şi în planul x-z. 

Calculul funcţiilor de formă este prezentat pe larg în Anexa 1, în 
urma căruia s-au obţinut: 

Funcţii de formă pentru deplasări liniare, valabile atât în planul x-y 
cât şi în planul x-z sunt: 



< 



w 



K0=N 1 feK -N 2 (^)cp i +N 3 (^)w j -N 4 (^)(p j =|Njk S } 

«pfeMifeV. - N,©*, +N 3 fek -N.fefo =LnJK s } 



(3.53) 



Calculul funcţiilor de formă 




'^(0=1-3^+2^ 



(3.54) 



N 4 fe)=L(-$ 2 +$ 3 ) 
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Funcţiile de formă pentru rotiri, notate N; (^) în planul x-y, diferă 
doar prin semnul „minus" de funcţiile de formă pentru rotiri din planul x-z 
notate N i(p fe). 

Funcţii de formă pentru rotiri în planul x-z sunt: 



Deci 



N 3 <pft) = 
N 4v fe) = 



3$ 2 -4$ + l]=-ff 2 fe> 
i(-6^ 2 +6^) =-NM 
3^ 2 -2^]=-N 4 fe> 



(3.55) 



KJ=-LnJ 



(3.56) 
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3.1.1.2. Elementul Timoshenko 



Elementul cilindric de arbore Timoshenko are două noduri la 
capete, având fiecare câte patru grade de libertate: două translaţii şi două rotaţii. 
Teoria de bară Timoshenko se bazează pe următoarele ipoteze: 

a) Ipoteza secţiunilor plane (a lui Bernoulli): o secţiune plană şi normală la 
axa barei înainte de deformare rămâne plană şi după deformarea barei; 

b) Deplanarea secţiunii transversale se neglijează; 

c) Se consideră o tensiune de forfecare medie, independentă de y. 

Ipoteza secţiunilor plane introduce rigidităţi false relativ la 
forfecare. Astfel, deformaţiile de forfecare sunt subestimate şi de aceea relaţia 
pentru forţa tăietoare T nu este corectă. Atunci, soluţia este folosirea în calcule 
a unei „arii efective de forfecare" mai mică decât aria reală a secţiunii, A S <A. 
Din considerente energetice, se poate scrie: 




2 J L GA S 



(3.57) 



unde 



Jt dA 



(3.58) 



Pentru bara de secţiune inelară, având raza interioară r şi raza exterioară R 
factorul % are următoarea expresie [34]: 



% = 



(3.59) 



7 + 6v 20 + 12v 

+ 

6 + 6v 6 + 6v 



r/R 
1 + fr/R) 2 



In această expresie v reprezintă coeficientul de contracţie transversală 
(coeficientul lui Poisson). Pentru oţel el are valoarea v = 0.3, atunci % pentru 
bara de oţel, de secţiune inelară are următoarea expresie [34]: 



1 = 



1.13 + 3.03 



' r/R ^ 
l + (r/R) 2 



(3.60) 
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Ecuaţiile de mişcare în planul x-v 



Fie un element de bară de lungime dx, situat în planul x-y, izolat 
dintr-un arbore (Fig.3.6). Asupra acestui element acţionează, în plus faţă de 
elementul Bernoulli-Euler, cuplul de inerţie la rotaţie p, -\|> -dx . 



Mz +Mz' dx 




a. Ecuaţii de echilibru 



(2X(2) = O) M z - T y dx - M vdx Ş + £ydx - (V + M z 'dx ) = 0 



(lX=°) T y+ ^dx-[T y+T ;dx] = 0 



(3.61) 



unde 



A = pi 



.. d 2 w 
V = — y 
dt 2 



(3.62) 
(3.63) 



iar după simplificări rezultă: 



M z +T y -Av = 0 
-T y '+^v = 0 



(3.64) 
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b. Ecuaţiile deformatei elastice 



xy 



unde 7^ reprezintă lunecarea specifică în planul x-y, de forma: 



Yxy + 



du dv 

— + — 
dy 8x 



(3.65) 



(3.66) 



c. Relaţii cinematice 

Ţinând cont de figura 3.7 se poate scrie: 




Fig.3.7 



Jxy =V -\|/ 



(3.67) 



d. Ecuaţiile de mişcare 



Introducând (3.67) în (3.65) se obţine: 

T=GA s (v'- ¥ ) 



(3.68) 
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Derivând relaţia (3.68), o dată, în raport cu x , rezultă: 



L T y '=GA s (v"- ¥ ') 



(3.69) 



Relaţia (3.69) se introduce în (3.64), de unde se obţine: 

EI z \\f" + GA s (v' - \|/)+ |Iv|/ = 0 
-GA s (v"-\|/')+nv = 0 



(3.70) 



sau, altfel scris, se obţin ecuaţiile de mişcare: 



EI z y"-GA s (y-v')+Ay = 0 
GA s (\|/'-v")+|liv = 0. 



(3.71) 



într-adevăr, deplasările axiale 

u = -\|/y 



\ < 




Fig.3.8 

Atunci, lunecarea specifică 

du d\ _ d(-\|/y) d\ __ 



dy dx dy 



dx 



-\\r + v 



(3.72) 



(3.73) 
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Rezultă şi relaţiile: 



(3.74) 

dx dx 



M 



= -Ja x ydA (3.75) 



z 

A 



Neglijând termenul p,v|> în relaţia (3.71), aceasta devine: 

FJ z \|/ " - GA s (\|/ - v ') = 0 (3 .76) 



sau 



FT 

M>-v'= ' (3.77) 
(jA s 



adică 



FT 

v' = i|/- z (3.78) 
(jA s 

Aceasta introduce o relaţie statică (aproximativă) între v' şi \|/, (de fapt 
mărimi cinematice independente) care permite simplificarea calculelor 
ulterioare. Ea reprezintă principala aproximaţie introdusă în formularea 
elementelor finite de tip Timoshenko. 

Se poate introduce o coordonată adimensională ^ şi se poate 
exprima ecuaţia de mişcare în funcţie de aceasta: 

ţ> = ^ (3.79) 

Atunci 

d d d£ 1 d 

Deci 



(3.80) 
dx ă% dx L d£ 



v' = *=I.* (3.81) 
dx L d£ 
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dx L d^ 



(3.82) 



„ d\\r' d 

\|/ = — — = — 
dx dx 



'1 d V )_l l d^ = ^.d^ (3 83) 



Ld^J L L d^ 2 L d^ 2 



Introducând (3.81) şi (3.83) în (3.78) se obţine: 



1 dv 1 :i , 1 d 2 \|/ 



(3.84) 



Notând 



K=i2 n * 



GA S L 



(3.85) 



Ecuaţia (3.78) pentru elementul Timoshenko în coordonate adimensionale 
corespunzătoare planului x-y este: 



1 dv K d> 

= i|/ v 

Ld^ Y 12 d^ 2 



(3.86) 



Ecuaţiile de mişcare în planul x-z 

Fie un element de bară de lungime dx, situat în planul x-z, izolat 
dintr-un arbore (Fig.3.9). Asupra acestui element acţionează, în plus faţă de 
elementul Bernoulli-Euler, cuplul de inerţie la rotaţie ji-cp-dx . 

a. Ecuaţii de echilibru 

Pentru elementul de bară izolat se scriu sumă de momente 
încovoietoare faţă de axa y şi sumă de forţe tăietoare paralele cu axa z . 
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(X M y(2) = O) M y + T z dx - nwdx y + £q)dx - ^M v + M v 'dxj = 0 



(3.87) 



(£f z = O) T z + inwdx - (j z + T z 'dx) = 0 




Tz + Tz' dx 



Fig.3.9 



unde 



A = pi 

.. _ d 2 cp 



dt 



iar după simplificări rezultă: 



M y -T z -Acp = 0 
[T z '+^iw = 0 



b. Ecuaţiile deformatei elastice 



My + My' dx 



lT z =GA sYj 



unde reprezintă lunecarea specifică în planul x-z, de forma: 



(3.88) 



(3.89) 



(3.90) 



du dw 

OL dx 



(3.91) 
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c. Relaţii cinematice 

Ţinând cont de figura 3.10 se poate scrie: 



Yxz =<P + W 



(3.92) 




Fig.3.10 
d. Ecuaţiile de mişcare 



Introducând (3.92) în (3.90) se obţine: 



M y =EI y cp' 
[T z =GA s (cp + w') 



(3.93) 



Derivând relaţia (3.93) o dată în raport cu x, rezultă: 



M y =my 

T z '=GA s ((p' + v") 



(3.94) 



Relaţia (3.94) se introduce în (3.89), de unde se obţine: 
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EI y cp " - GA s (cp + w ' ) - £cp = 0 
GA s (q>' + w w )-|iw = 0. 



In ecuaţiile (3.87) - (3.95) s-au folosit următoarele relaţii: 



u = cp-z 



Yxz 



du d\\ d((pz) d\\ , 
— + — = v ^ 7 + — = (p + w 
dz dx dz dx 



du d((Dz) , 

8 X = = v ' = Z • (p 

dx dx 



M y = Jcy x z dA 

A 

M y =e y( p' 



(3.95) 



(3.96) 
(3.97) 

(3.98) 
(3.99) 




Fig.3.11 



Neglijând termenul jCtcp în relaţia (3.95) aceasta devine: 

EI cp" - GA s (cp + w') = 0 



(3.100) 
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sau 



adică 



EI 

(p + w' = — — (p" 
1 GA o 



w = -cp + 



EI 



GA 



Z y„ 



(3.101) 



(3.102) 



sau 



EI 



w 



'=(-<p)-^M-<p) 



GA, 



(3.103) 



Utilizând coordonata adimensională £ = — se poate exprima 
ecuaţia (3.103) în funcţie de aceasta. Atunci, din relaţia (3.80) 



w = 



dw 1 dw 



dx L d£ 



("<p) = 



_d(-<p)_l d(-q>) 



dx 



(3.104) 



(3.105) 



(-*) = 



" _d(-q>)' d ( 1 d(-(p)^_ 1 1 d 2 (-cp)_ 1 d 2 (-(p) 



dx dx 



L d£ J L L d^ 2 



d^ 2 



(3.106) 



Introducând (3.104) şi (3.106) în (3.103) se obţine: 



ldw = ( 0 ) 11 > 1 d 2 (-(p) 
Ld^ V V) GA S L 2 d^ 2 



(3.107) 



Notând 



EI 

K = 12 y 



GA S L 



(3.108) 



Ecuaţia (3.107) pentru elementul Timoshenko în coordonate adimensionale 
corespunzătoare planului x-z este: 

I d ^ = (_ (p )_K^P) (3 . 10 9) 
L d£ v Y/ 12 d£ 2 
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Funcţii de formă pentru elementul Timoshenko în planul x-y 



Pentru aproximarea deplasărilor şi rotirilor elementului de bară, se 
folosesc funcţii polinomiale de gradul trei de forma: 



v(0 = A3^+A 2 ^ 2 +A^ + A 
W0=B^ 3 +B 2 ^ 2 +B^ + B t 



(3.110) 



Se înlocuiesc relaţiile (3.110) în ecuaţia de mişcare a elementului 
din planul x-y (3.86). 



Dar 



Rezultă că: 



^ = 3A 3 4 2 +2A 2 4 + A 1 

^ = 363^+26^ + 1*! (3.111) 

^ = 66 3 ^ + 26 2 
d^ 2 3 



i(3A 3 ^ 2 +2A 2 ^ + A 1 )=B 3 ^ 3 +B 2 £ 2 +B^ + B 0 -^(6B 3 ^ + 2B 2 ) 

(3.112) 



Ordonând şi în membrul drept şi în membrul stâng după puterile lui ^ , rezultă: 



3A 



3£2 + —2 



2A, A^_ 3+B + 
L L 3 



KB 3 ^ 



B, 



KB 2 ^ 



(3.113) 



iar prin identificare se obţine: 
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B 3 =0 
3 A 



= B 



2A 2 KB 3 



L 

A* 

l L 



„ KB 2 



B 3 =0 

B 2 =-A 3 
2 L 3 

Bi=-A 2 
1 L 2 

■•=r 



Ai+yA 3 



astfel, pentru că B 3 = 0 , funcţia este de gradul doi: 



¥fe) = B 2^ 2 +B 1 4 + B 0 = 

= ^A 3 ^ 2 +^A 2 ^ + i| A x + Şa 3 



(3.114) 



(3.115) 



Deplasarea liniară v şi rotirea \\f se aproximează în funcţie de 
deplasările nodale, cu ajutorul funcţiilor de formă statice N, în felul următor: 



v(0=N 1 (4)v i +N 2 0;Vi +N 3 fe)v j +N 4 fe>Kj =Ln_ 
Yfe) = N t fejv, + N 2 feVi + N 3 fe)v j + N 4 fe>|/ j = Ln 



u. 



(3.116) 



unde am folosit notaţia (3.3) şi în plus 



LnJ=LNx 
[n]=[n 1 



N, N, N 



4. 



N, N, N 



4_ 



(3.117) 



Calculul funcţiilor de formă 



Funcţiile de formă au următoarele proprietăţi importante: 

1 . Suma funcţiilor de formă are valoarea egală cu 1 în orice nod; 

2. Corespunzător coordonatei n, funcţia de formă N n =l, iar celelalte 
funcţii de formă au valori nule. 

Pentru \ { = v(x = o)= v(^ = 0) -> N, (^) = 1;N, (£)= 1; 

Pentru Vl = i|/(x = 0) = = O) N 2 fe) = l;N 2 fe)= h 

Pentru Vj =v(x = L)=v(^ = l) ->N 3 (5) = 1;N 3 (£) = 1; (3.118) 

Pentru Vj = ¥ (x = L) = ¥ fe = l) -^N 4 (^)=l;N 4 fe) = l. 

Calculul funcţiilor de formă este prezentat pe larg în Anexa 2. 
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Funcţii de formă pentru elementul Timoshenko în planul x-z 

Ecuaţia de legătură între w şi (p pentru elementul Timoshenko în 
planul x-z este dată de relaţia (3.109): 

ldw, v K d 2 (-cp) 
Ld^ V } 12 d^ 2 

Pentru aproximarea deplasărilor şi rotirilor elementului de bară în 
acest plan, se folosesc tot funcţii polinomiale de gradul trei de forma: 



w(0=A34 3 +A 2 ^ 2 +A^ + A 0 
[-(pfe) = B3^ + B 2 ^ 2 + B^ + B 0 



(3.119) 



Se înlocuiesc relaţiile (3.119) în ecuaţia de mişcare a elementului 
din planul x-z . 



Dar 



^ = 3A 3 ^ 2 + 2A 2 ^ + A 1 
^P) = 3B 3 ^ + 2B 2 ^ + B 1 

^P) = 6B 3 ^ + 2B 2 
d£ 2 



(3.120) 



Rezultă că: 



i(3A 3 ^ 2 +2A 2 ^ + A 1 )=B 3 ^ 3 + B 2 £ 2 +B^ + B 0 -^(6B 3 ^ + 2B 2 ) 

Ordonând şi în membrul drept şi în membrul stâng după puterile lui ^ , rezultă: 



3A 



3 &2 



2A 2 . , A, 



KB 



¥+ H+ — = B 3 r +B 2 ^ Z + B x 



KB 
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iar prin identificare se obţin pentru constantele A-, Bj cu ( i = 1,2,3,4 ) aceleaşi 
valori ca şi în planul x-y: 



"3 

3 A 

L 

2A 

^-B ™ 2 B -Ua +î« ' 

lT -B 0 -— [ B «- L [ A i + 2 ^ y 



= 0 






B 3 


= 0 


>=B 2 






B 2 


= ^A 3 
L 3 




KB 3 
2 


=> - 


B x 




= B 0 - 


KB 2 
6 




B 0 





astfel încât, pentru că B 3 = 0 şi funcţia [ -cp(%) ] este de gradul doi: 



-<P(0 = B ^ 2 +B^ + B 0 = 



K * 1 

3 



De asemenea şi deplasarea w şi rotirea (- cp) se aproximează în 
funcţie de deplasările nodale, cu ajutorul funcţiilor de formă statice N: 



"feM^K -N^fo + N 3 ©Wj -M&i = LnJ k s } 

-cpfe) = N 1 fe)wi +N 3 fe)wj "N 4 fe>Pj =LnJK S ) 



unde am folosit notaţia (3.3) şi la fel 

LnJ=Ln x n 2 n 3 n 4 . 
[nJ=[n 1 n 2 n 3 n 4 _ 

Calculul funcţiilor de formă 

Folosind proprietăţile funcţiilor de formă, relaţiile (3.118), se 
calculează expresiile acestora şi în planul x-z. 

Calculul funcţiilor de formă este prezentat pe larg în Anexa 2. 
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Funcţii de formă pentru deplasări, valabile atât în planul x-y cât şi 
în planul x-z sunt: 



Ni ® = îiK[ 1 " 3 ^ 2+2 ^ +K(1 "^ 

N3ft)= î7K^ 2 " 2 ^ +K ^ 



(3.122) 



Funcţiile de formă pentru rotiri, notate Nj (%) în planul x-y, diferă 
doar prin semnul „minus" de funcţiile de formă pentru rotiri din planul x-z 
notate 5Lfe). 



Funcţii de formă pentru rotiri în planul x-z sunt: 





1 

1 + K 








1 

1 + K 


[ 3 ^-4^ + l + K(l-O]=-N 2 0a 


N 3(p ft) = - 


1 

1 + K 


"^(-6^+6^) 


= -N 3 fe> 




1 

1 + K 


fe 2 -2$ + K$]= 


= -N 4 fe> 



(3.123) 



Deci relaţia (3.56) este valabilă şi pentru funcţiile de formă corespunzătoare 
elementului Timoshenko: 



KJ=-LnJ 
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Este utilă introducerea unui al treilea set de funcţii de formă care se 
folosesc la calculul matricei de rigiditate şi se definesc astfel: 



Acestea sunt: 



H=H-r 



dN 



(3.123a) 



"N 3 fe)= i — 



N 2 fe)= N 4 fe) = 



K 

1 K 

2 1 + K 



(3.123b) 
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3.1.1.3. Matricea de masă, matricea giroscopică şi matricea de rigiditate 
pentru elementul de arbore cilindric 

In cele ce urmează, matricele de masă, giroscopică şi de rigiditate 
pentru un element de arbore se stabilesc pornind de la expresiile energiei 
cinetice şi potenţiale, folosind ecuaţiile Lagrange: 



dt 



dT av 



+ ^ = Qi, (i = l,2,-,n) (3.124) 



unde 

T - energia cinetică; 

V - energia potenţială; 

q; - coordonatele generalizate; 

q; - vitezele generalizate; 

Qi - forţele generalizate. 



Energia cinetică pentru un element de arbore se poate scrie: 
T s = im s (v 2 + w 2 )+i(j x co x 2 + J y Q) y 2 + J z co z 2 ) (3.125) 



unde 



J x = Jp şi J y = Jz = Jt sunt momentele de inerţie principale ale 
arborelui faţă de sistemul de axe mobil xyz, fixat pe arbore; 
m s - masa elementului de arbore; 

J P - momentul de inerţie masic polar al elementului de arbore; 

J T - momentul de inerţie masic transversal al elementului de arbore. 

co x ,0) y ,a) z - componentele vitezei unghiulare în sistemul de axe 

mobil xyz. Acestea au următoarea expresie: 



* > 


f 




> = < 











Q + (j>\|/ 
(pcosG + \j/sin9 
\j/cos0 - cpsinG 



m care 
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Q = viteza unghiulară cu care arborele se roteşte în jurul axei x; 
0 = Ot. 



Atunci, expresia energiei cinetice este: 



T s = im s (v 2 + w 2 )+ 1 J P S (O + cpy) 2 + 

+ ^ J x s [((pcos0 + \j/sin0) 2 + (\j/cos0 - cpsin©) 2 ] 



(3.126) 



iar după desfacerea parantezelor, devine: 

T s = im s (v 2 + w 2 )+ij P s Q 2 + J P s Qcpy + ij P s cpV 2 +^Jx S (cp 2 +¥ 2 ) 



v ' v v ' efect 

energia cinetica energie de giroscopic 
de translaţie rotatiect. 



=0 



energie cinetica 
de rotaţie 



(3.127) 



Pentru un element de arbore de lungime infinit mică dx, energia 



cinetică este: 



dT s = 



± P a(v 2 +w 2 )+ipI P (Q 2 +2Q(p ¥ )+ipl((p 2 +xj/ 2 ) 



dx (3.128) 



unde 



m s = \i dx = pA dx , 
J P =pl p dx, 
J T = pldx, 



|u = pA. 



T 1 I S 



59 



Dana-Codruţa Vişan Modelarea dinamică a sistemelor rotor-lagăre 



T 1 T S 

fX p =pl p =-J P 



Folosind aceste relaţii se poate scrie: 



dT s = 



M( V 2 + ^2) + A( ( p2 +vj/ 2) + ^_ Jp S Q2+ApQ(>x(/ 



dx (3.129) 



Integrând expresia energiei cinetice pentru elementul de arbore, se obţine: 
T s = ^ J (v 2 + w 2 )dx + & J («cp 2 + \j/ 2 )dx + -L J P S Q 2 J dx + ^ p qJ cpydx, 



Deci 



T s = ^J(v 2 + w 2 )dx + ^j(cp 2 +\j/ 2 )dx + ApQj(p\|/dx + - J P S Q 2 . 



(3.130) 



Energia de translaţie 



translaţie: 



Primul termen din relaţia (3.130) reprezintă energia cinetică de 



(3.131) 



In această relaţie, vitezele se pot scrie în funcţie de coordonatele nodale şi de 
funcţiile de formă: 



=LNJk s j 



w 



=LnJK s J 



v = LnJi/ }= v T 4 ' J LN f * = [Njfk.' }= = fa Jm 
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După înlocuire rezultă: 

t, s = \ k F • Jln j 1 ln> ■ K s K * k F • fJLN j 1 ln> • K s } 




(3.132) 

x 1 

Se trece la coordonata adimensională E = — . Atunci d£ = — dx, adică 

dx = Ld^ , astfel încât submatricea de masă de translaţie pentru elementul de 
arbore este: 

[m T s ]=^|_N_TI_Njdx, 

o 

sau 

[m T s ]= mL JLN(0J T LNfe)J d4 • (3.133) 

o 

Pentru că funcţiile de formă de translaţie au expresii identice atât în planul x-y 
cât şi în planul x-z, submatricea de masă de translaţie are aceeaşi formă în 
ambele plane. 

Energia de rotaţie 

Al doilea termen din relaţia (3.130) reprezintă energia cinetică de 

rotaţie: 




(3.134) 



în această relaţie, derivatele în raport cu timpul se pot scrie în funcţie de 
coordonatele nodale şi de funcţiile de formă: 
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v-lsJk 8 J 

v=LN| 1; s }=y T = K s RsJ r 



♦--L»V}-* T -^rLsJ 



După înlocuire rezultă: 

o o 

V / ✓ 

[ m R S ] [ m R S ] 

(3.135) 
x 1 

Se trece, de asemenea, la coordonata adimensională % = —, având d^ = — dx, 

L L 

adică dx = Ld^ , astfel încât submatricea de masă de rotaţie pentru elementul de 
arbore este: 

[m/^^jLNfLNjdx, 

o 

deci 

[m R s ]= AL/NOT Ln(^ . (3.136) 
o 

Submatricea de masă de rotaţie, pentru un element de arbore, are aceeaşi formă 
în ambele plane x-y şi y-z. 

Submatricea de masă pentru un element de arbore are forma 

generală: 

[m s ]=[m T s ]+[m R s ], (3.137) 

adică 

[m s ]= ^LjLN(0J r LN(0J^ + ALjLN(0f Ln(0J^ ( 3 - 138 > 
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Energia datorită efectului giroscopic (Coriolis) 

Al treilea termen din relaţia (3.130) reprezintă energia datorită 
efectului giroscopic: 

L 

T 3 S =^pQj(p\(/dx (3.139) 
o 

Introducând pe (p şi pe \|/ , rezultă: 

t/= -4.» W\ 

o 




Submatricea giroscopică are, deci, forma: 

o 

şi trecând la coordonata adimensională ^ , 

[g S ]= ApLjLNfâfLNfelte = 2 k1 ( 314 °) 
o 

Energia cinetică pentru un element de arbore cilindric se poate 

scrie: 

(3.141) 

înlocuind această expresie în ecuaţiile Lagrange, se obţine: 




(3.142) 
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sau 

[m S ]{u}+q[g S ]{ii}={f s } (3.143) 

unde matricea globală de masă J este simetrică, de ordin 8x8, iar matricea 
globală giroscopică G este antisimetrică, de ordin 8x8. 



Energia potenţială de deformaţie 

Pentru un element de arbore infinit mic, energia potenţială de 
deformaţie, în cazul general, are forma: 



dU = -<! 
2 



EI 



( d\|/Y ( dcp 



V dxy 



dx 



+ GA s (y xy 2 +y xz 2 ) 



dx 



(3.144) 



Integrând şi înlocuind deformaţiile de forfecare din relaţiile (3.67) şi (3.92), se 
obţine: 

U = — j(y' 2 +9' 2 )^ + ^j[(v'-V|/) 2 +(w' + 9) 2 K ( 3145 > 



a. Energia de încovoiere 



Primul termen din relaţia (3.145) reprezintă energia datorită 



încovoierii: 



U 1= ^J(v' 2 + cp' 2 )dx = 



(3.146) 



= ^' 2 dx + ^J(-cp') 2 dx 



EI 
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Se trece la coordonata adimensională £ folosită şi anterior şi se face notaţia: 



(■ ) V =|(- ) (3-147) 



dar 



Y dx Ldţ L L J " ' y ds Ldţ L L Jr ' ' 



¥ ' 2 =l |/'V= <p' 2 =q)'V = 



= LN>/j w = LnjV z s } 



dv = ldv = ll .Jfs) w -_*_I*_I| N vJLs) 
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Astfel încât, pentru energia de încovoiere, se ajunge la expresia: 



(3.148) 



Submatricea de rigiditate de încovoiere , utilizând coordonata adimensională \ 
este: 



dN 


i 


dN 









(3.149) 



b. Energia de forfecare 

Al doilea termen din relaţia (3.145) reprezintă energia datorită 



forfecării: 



GA C r / A2 , GA 1 



J(\|/-v') 2 dx + ^^J((p + w') 2 dx (3.150) 



u 2 = — 1V . , 

2 2 v, , 2 



unde 



¥ - v = InJ-ln'J^/ }= f [sJ4In v J]k s }= L4, s } 



<P + W = -CnJ- iN'Jk' }= # J" V J)K S }= -Ln|, z s }. 



Atunci 



U 2 a^f.GA s LiLNj I LNk-k S K'k S r-GA s LiLNj I LNk-k S } 



(3.151) 
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Submatricea de rigiditate de forfecare, utilizând coordonata adimensională \ 
este: 

1 . 

1 ^ II I /v I 

(3.152) 

o 



k s ]=GA s LjLNf[N.k 



generală: 



Energia totală de deformaţie pentru un element de arbore este: 
Submatricea de rigiditate pentru un element de arbore are forma 



adică 



:]♦[«< :]■ 



(3.154) 



dN 


T 


dN 









d^ + GAsLjLrvfLNji^ (3-155) 



înlocuind energia U în ecuaţiile Lagrange, se obţine: 



V] 



[0] 

FU 



u 



(3.156.a) 



sau 



[k s ]{»}=M 



(3.156.b) 



unde J este matricea globală de rigiditate, de ordin 8x8, pentru un element 

de arbore cilindric. 

Deci, combinând (3.143) şi (3.156.b) rezultă ecuaţiile de mişcare în 
funcţie de matricele globale: 

[m S ]{u}+q[g S ]{u}+[k S ] = {f s } (3.157) 
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Pentru stabilirea formelor explicite ale submatricelor de masă, 
giroscopică şi de rigiditate, este nevoie de calcularea integralelor următoare, 
pentru fiecare tip de element de arbore în parte: 



JLN(iOJ T LN(i0Jdi; 

o 



J 



dN 


T 


dN 









dl 



In Anexa 3 este prezentată pe larg rezolvarea acestor integrale 
pentru calculul submatricelor corespunzătoare elementului de arbore Bernoulli- 
Euler, iar în Anexa 4 este prezentat calculul submatricelor corespunzătoare 
elementului de arbore Timoshenko. 



Se obţin următoarele: 



Submatricea coerentă de masă pentru elementul de arbore 
Bernoulli-Euler este: 



156 22L 54 -13L 

4L 2 13L -3L 2 

156 -22L 

sim 4L 2 



(3.158) 
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Submatricea de rigiditate pentru elementul de arbore Bernoulli- 

Euler este: 



0 L 



dN 


T 


dN 









EI 

L 1 



12 6L 

i 4L 2 



-12 6L 
-6L 2L 2 
12 -6L 



sim 



4L 2 



(3.159) 



Submatricea coerentă de masă pentru elementul de arbore 
Timoshenko are următoarea expresie: 



k]= 



(1 + K) 2 



420 



156 22L 54 
4L 2 13L 
156 



-13L 
-3L 2 
-22L 
4L 2 



+ K 



120 



36 


-9L 


9L 


-2L 2 


84 


-11L 




2L 2 



+ K 



2 

120 



40 5L 


20 


-5L 






36 


3L 


-36 


3L 


; l 2 


5L 2 


-L 2 


+ 






4L 2 


-3L 


-L 2 




40 


-5L 

L 2 


30L 2 






36 


-3L 
4L 2 



+ 



+ K 



pL 

6Î7 



0 -3L 


0 


-3L 


; l 2 


3L 


-L 2 




0 


3L 






L 2 



K 



2 AL 
6L 2 



0 0 0 0 

1 2L 2 0 L 2 

0 o 

; 2l 2 



(3.160) 



69 



Dana-Codruţa Vişan Modelarea dinamică a sistemelor rotor-lagăre 



Sau altfel scrisă, 



m 



11 



sim 



, N 






111 12 




m s u 








m 22 


™23 


m s 24 




™ s 






m 33 


m s 34 












m s 44 



unde: 



= (l56 + 294K + 140K 2 )tt T + 36a R 
m^ 2 =(22 + 38.5K + 17.5K 2 )La T +(3-15K)La R 
m s 13 = (54 + 126K + 70K 2 )tx T - 36a R 
m^ 4 =-(l3 + 31.5K + 17.5K 2 )La T +(3-15K)La R 
m L = (4 + 7K + 3.5K 2 )L 2 a T + (4 + 5K + 10K 2 )L 2 a R 



___ S S 
m 23 = - 111 14 



m 24 = -(3 + 7K + 3.5K 2 )L 2 a T - (l + 5K - 5K 2 )l 2 cc r 



m 33 = m n 



___ S ^ S 

m 34 = ~ m i2 



s s 

111 44 = 111 22 



(3.161) 



a , = 



(j,L 
420(l + k) 2 



(3.162) 



a„ = 



JLlL 



R 30L 2 (l + K) 2 



(3.163) 
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Submatricea giroscopică pentru elementul de arbore Timoshenko: 



(1 + Kf 



30L 2 



36 3L 

i 4L 2 



-36 3L 
-3L -L 2 
36 -3L 
4L 2 



sau 



+ K 



6L 2 



0 - 



3L 0 
L 2 3L 
0 



-3L 
-L 2 
3L 

L 2 



+ K 



2 A L 



6L 2 



0 0 0 0 

1 2L 2 0 L 2 

o o 

; 2l 2 



(3.164) 



A 

§22 



A 

§23 
§33 



sim 



§14 
§24 
§34 
§44 



(3.165) 



în care, cu aceleaşi notaţii făcute anterior (3.162), (3.163) pentru a x şi a R , 
elementele submatricei giroscopice au următoarele valori: 



A 


= 72cx R 


§23 


— §12 


§12 


= 2(3-15K)La R 


§24 


= -2(l+5K 


A 


= "72a R 


§33 


= A 


§14 


= 2(3-15K)La R =gL 


§34 


~ §23 


§22 


= 2(4 + 5K + 10K 2 )L 2 a R 


§44 


= §22 
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Submatricea de rigiditate pentru elementul de arbore Timoshenko 



este: 



L J 1 + K 



EI 

V 



12 6L 

i 4L 2 



-12 6L 
-6L 2L 2 
12 -6L 
4L 2 



+ K 



0 0 0 

; l 2 o 
; o 



o 

-L 2 
0 
L 2 



(3.166) 
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3.1.1.4. Vectorul forţelor de dezechilibru masic pentru arborele cilindric 

Se consideră o distribuţie liniară a dezechilibrului masic în lungul 
elementului de arbore, între a CL şi a SL la capătul din stânga şi a CR şi a SR la 
capătul din dreapta, de forma: 



acfe) = a c L (l-0 +a cR^ 

as(0= a SL(l-0+ a SR^ 



(3.167) 



unde 



a CL =a L cosa L , 
a SL =a L sina L , 
a cR = a R cosa R , 
a sR =a R sina R , 



iar 



(3.168) 



a L , a R sunt razele de dezechilibru, a L , a R sunt fazele dezechilibrului. 



Sarcinile distribuite, corespunzătoare dezechilibrului masic sunt: 

; Py fet)=nn 2 a;fet) 

lp z fet)=^ 2 affet) 



(3.169) 



unde 



a* (ţ, t) = a* ©cos Qt - a* (^)sinQt 
a z 0 = a s fe) cos Qt + a c fe) sinQt 



(3.170) 



Lucrul mecanic virtual al acestor forţe este: 



5W = J(8v) T p y dx + J(8w) T p z dx = 

o o 

= fa F jLNjp y • L • d^ + £uŞ } T |LnJ t p z • L • d$ 

0 gj ' ,r TT" " 



(3.171) 
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Forţele nodale, echivalente cinematic sunt: 



{f J s }=HLn 2 )LNj r a^.dţ 

0 

{f z s }=nLQ 2 jLNj T ^ 



(3.172) 



Vectorul dezechilibrului masic are expresia: 



fel 



= |iLQ 2 



o 



cosQt + jiLQ 2 



-JLnJ t 4-«iş 

o 



sinQt . 



Exemplu de calcul 



J N, • d^ = J N, feXa CL (l " 0+ a CR ^ = 



=a CL jN 1 feXl-^ + a CR |N 1 (^ = 



(3.173) 



^— J(42 + 40K)a CL +(l8 + 20K)a CR ]. 



120(l + K) 
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Asemănător se calculează şi celelalte integrale corespunzătoare celorlalte funcţii 
de formă. Se ajunge la: 



* S 120(l + 



(1 + K) 



a 



CL 



a 



CR J 



a 



SL 



L a SR J 



cosQt + 



-M 



a 



SL 



L a SR j 



a 



CL 



a 



CR J 



sinQt 



(3.174) 



unde 



astfel încât 



W= 



42 + 40K 
(6 + 5K)L 
18 + 20K 
-(4 + 5K)l 



18 + 20K 
(4 + 5K)L 
42 + 40K 
-(6 + 5K)L 



(3.175) 







-te}' 




>cosQt + < 




.tel 




. tel' 



sinQt . 



(3.176) 
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Vectorul forţelor de dezechilibru masic pentru elementul de arbore 
Bernoulli-Euler 

Folosind în calcule funcţiile de formă pentru elementul Bernoulli- 
Euler, rezultă: 



42 


18 


6L 


4L 


18 


42 


-4L 


-6L 



(3.177) 



care conduce la: 







f 


Ui 














u 2 L 




-o 6 L 










u 3 




-o 7 






HLQ 2 


< 


u 4 L 


► cosQt + < 


-u 8 L 


► sinQt 


120 


















u 6 L 




u 2 L 










u 7 




u 3 








V 


u 8 L 




. u 4 L . 


) 



(3.178) 



unde 

u 1 =42a CL +18a, 
u 2 =6a CL +4a CR 
u 3 =18a CL +42a 
-u 4 =4a CL +6a ( 



w 5 =42a SL +18a SE 

«6= 6a SL+ 4a SR 

u 7 =18a SL +42a SE 
-u 8 = 4a sL +6a SR 
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Vectorul forţelor de dezechilibru masic pentru elementul de arbore 
Timoshenko 

După efectuarea calculelor se ajunge la următoarea formă: 



ls/ 120(1 + K) 







-w 5 


u 2 L 




-o 6 L 


» 3 




-u 7 




►cosQt + < 


-u 8 L 


u 6 L 




u 2 L 


o 7 






u 8 L 




u 4 L 



■sinQt 



(3.179) 



J 



unde 



u t = (42 + 40K)a CL + (l8 + 20K)a 1 
u 2 = (6 + 5K)a CL + (4 + 5K)a CR 
u 3 = (18 + 20K)a CL + (42 + 40K)a 
- u 4 = (4 + 5K)a CL + (6 + 5K)a CR 



o 5 =(42 + 401^ + (18 + 20K)a s 
u 6 = (6 + 5K)a SL + (4 + 5K)a SR 
u 7 =(l8 + 20K)a SL +(42 + 40K)a 
- u 8 = (4 + SK)^ + (6 + SKja^ 
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3.1.2. Elemente tronconice de bară cu două noduri 

In cele ce urmează se prezintă două tipuri de elemente finite conice, 
concepute pornind de la elementele finite cilindrice de tip Timoshenko, cu două 
noduri la capete, cu matrice de masă şi matrice giroscopică coerente. 

3.1.2.1. Elementul CON 12 

Greenhill, Bickford şi Nelson [70] au extins elementul finit 
cilindric de tip Timoshenko, la elementul finit conic şi au luat în calcul efectele 
inerţiei la rotaţie, la translaţie, efectul momentelor giroscopice, a deformaţiilor 
de forfecare şi de încovoiere, efectul încărcării axiale şi al amortizării interne. 

In descrierea elementului CON 12, inspirat din lucrarea [70], se 
neglijează efectul încărcării axiale şi al amortizării interne. 

Elementele conice înlocuiesc elementele cilindrice utilizate la 
modelarea în trepte a porţiunilor tronconice dintr-un arbore. Astfel, este posibilă 
compararea celor două tipuri de modelări. Rezultatele obţinute în urma utilizării 
elementului conic la modelarea rotorilor au fost mai bune în comparaţie cu cele 
obţinute în urma modelării în trepte a aceluiaşi rotor cu elemente cilindrice. 

Elementul CON 12 se bazează pe cel descris în [70]. El are 6 grade 
de libertate pe nod: deplasarea transversală, rotirea axei neutre şi unghiul de 
forfecare, în câte două plane ortogonale, rezultând matricele elementului de 
ordin 12x12. Acestor grade de libertate li se asociază următorul mod de variaţie: 
pentru deplasarea laterală o variaţie cubică şi pentru deformaţia de forfecare o 
variaţie liniară. 

In teoria de bară Timoshenko [140] se pune o condiţie suplimentară 
pentru neglijarea efectului inerţiei la rotaţie în ecuaţia diferenţială de echilibru a 
momentelor încovoietoare, relaţiile (3.76) şi (3.100). In absenţa termenului care 
ţine cont de inerţia la rotaţie, rotirea secţiunii transversale depinde de deplasarea 
laterală (altfel acestea sunt independente). Forţa tăietoare este constantă în 
lungul elementului. Ea este egală cu derivata întâi a momentului încovoietor , ca 
în cazul static. 

Cu funcţia de variaţie de gradul trei presupusă pentru deplasarea 
laterală, rotirea secţiunii transversale este o funcţie de gradul doi variind cu 
coordonata longitudinală şi momentul încovoietor variază liniar. 

Pentru ca elementul CON 12 să fie compatibil cu programele de 
calcul de dinamica rotorilor, se face o reducere Guyan [73] pentru condensarea 
matricelor elementului de la ordinul 12x12 la ordinul de mărime 8x8. 
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Pentru elementele conice, factorul de forfecare este calculat arbitrar 
folosind caracteristicile geometrice ale secţiunii transversale de la mijlocul 
elementului. Ca urmare, se observă că exactitatea analizei cu elemente finite 
conice depinde de alura elementului şi de alegerea unui diametru „efectiv" 
pentru calculul constantei de forfecare. O alternativă este utilizarea valorii medii 
pătratice a diametrelor de la capetele elementului, sau media geometrică a 
acestora [48]. Aceste variante par să se apropie mai mult de varianta ideală la 
care se doreşte utilizarea unui număr cât mai mic de elemente conice. 

Fie un element finit conic de lungime L, mărginit de nodurile i şi j, 
cu raza interioară r şi raza exterioară R, ca cel din figura 3.12. Acest element are 
12 grade de libertate, câte 6 pe fiecare nod: două translaţii, două rotaţii şi două 
deformaţii de forfecare. 



L 





X 








ce 


















® 

i— 




CU 


- q 




1 ce 



Fig.3.12 



Intr-o secţiune oarecare, raza interioară şi cea exterioară au 
următoarele forme: 

r = r i (l-^) + r j cî 

(3.180) 

R = R i (l-^) + R j ^. 

unde Ş, reprezintă coordonata adimensională de forma b= — . Sau, făcând 

L 

notaţiile: 

r R 

p = J si a = J . (3.181) 
r, R ; 
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(3.182) 



expresiile razelor devin următoarele: 

r = ii[l + (p-lfe] 

R = R i [l + (a-lX] 

Aria secţiunii transversale este: 

A = tt(r 2 - r 2 )= A, (l + ct^ + a 2 £ 2 ) (3.1 83) 



m care: 



A.^niR* -r*) (3.184) 



2 [ R-( CT -l)-r-(p-l) ] 

R= -r, 



^feţp. (3,86) 
R : -r : 



Asemănător se poate scrie şi expresia momentului de inerţie axial al 



secţiunii: 



i = ^^ = î^ = i,[ 1+5i ş +5 ^ +8 ^ +8 ^] 

(3.187) 

în care: 

t = lt (R, 4 -r, 4 ) (3l8g) 

s= 4 l R, 4 ( g -l)- r , 4 (p-l) J (31g9) 
R. -r : 
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S 2 = 



83 = 



84 = 



6\R*(a-if 


4 / 

~ r i (P 




Ri 4 


4 

-r. 
■ 




4|R. 4 (a-l) 


4 / 

~ r i (P 


"l)J 




4 

-r. 




k 4 (a-l) 4 


" r i 4 (p- 


Ol 



T» 4 4 

R : -r : 



(3.190) 



(3.191) 



(3.192) 



Vectorul deplasărilor pentru un astfel de element de arbore este: 

fc S } T ={ U l U 2 U 3 U 4 U 5 U 6 U 7 U 8 U 9 U 10 U ll U 12 f = 

= {v, 



w i «Pi ^ Vj Wj 9j Vj y xyi Yxzi y xyj y xzj } T 

(3.193) 



unde: 



cp = -w -y xz 
\|/ = v' + y xy . 

iar y xy si y xz reprezintă deformaţiile de forfecare, ca în figura 3.13: 



(3.194) 



y(v) 



\ 



\ 




Fig.3.13 
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Funcţiile de formă pentru deplasări sunt: 



[N] = 



unde 



"Nj o o 
0 N, -N 2 



N 2 N 3 



0 0 N, -N. 



0 -N, 



0 



0 



0 N 3 -N 4 0 



0 -N, 



0 -N, 



(3.195) 



^(0 = 1-3^+2^ 
N 2 (0 = L(^-2^ 2 + ^ 3 ) 
N 3 (0=3^ 2 -2^ 

n 4 ©=l(-^ + ^ 3 ) 



(3.196) 



sunt, de fapt, funcţiile de formă pentru deplasări din rel.(3.54) deduse pentru 
elementul cilindric Bernoulli-Euler. 



Funcţiile de formă pentru rotiri: 



[n]= 



unde 



0 N t N 2 



0 



0 -N 1 N 3 



0 



0 -N 4 0 -N< 



-N t 0 0 N 2 Nj 



0 



0 N, N, 



0 N„ 



0 



(3.197) 



N 2 =l-4^ + 3£ 2 
N 3 =-2^ + 3^ 2 
[N 4 =3^-3^ 2 



(3.198) 



sunt, de fapt, funcţiile de formă pentru rotiri din rel.(3.55) deduse pentru 
elementul cilindric Bernoulli-Euler. 



Funcţiile de formă pentru deformaţiile de forfecare: 



0 0000000 Xi o 
00000000 0 Xi 



%2 

o 1 



o 



(3.199) 
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unde 



(3.200) 



Dacă se foloseşte un sistem de referinţă mobil, care se roteşte cu 
viteza unghiulară de precesie co, atunci, folosind transformarea de coordonate 



ecuaţia de mişcare pentru un element de arbore conic se scrie: 



(K ] + [M R ]){p s } + |co([m t ]+ (1 - X\G^ }- 



{ [K B ]-« 2 ([M T ] + (l-2^Mj }{p s }={q c S } 



(3.201) 



(3.202) 



în care toate matricele sunt de ordinul 12x12, simetrice, în afară de matricea [G] 
care este antisimetrică, 



CO 



(3.203) 



iar [R] este o matrice bloc diagonală, în care fiecare element de dimensiuni 2x2 
este o matrice de rotaţie: 

cos cot - sin cot 
sin cot cos cot 



[m t ] este o matrice antisimetrică a cărei expresie se dă mai jos. 
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Matricele elementului 



Matricele elementare din această ecuaţie de mişcare au următoarele 



forme: 



Matricea de rigiditate la încovoiere: 



[K B ] = 



EL 



5040L 3 



(1,1) 




















0 


(2,2) 


















0 


(3,2) 


(3,3) 












sim 




(4,1) 


0 


0 


(4,4) 














(5,1) 


0 


0 


(5,4) 


(5,5) 












0 


(6,2) 


(6,3) 


0 


0 


(6,6) 










0 


(7,2) 


(7,3) 


0 


0 


(7,6) 


(7,7) 








(8,1) 


0 


0 


(8,4) 


(8,5) 


0 


0 


(8,8) 






(9,1) 


0 


0 


(9,4) 


(9,5) 


0 


0 


(9,8) 


(9,9) 




0 


(10,2) 


(10,3) 


0 


0 


(10,6) 


(10,7) 


0 


0 


(10,10) 


(11,1) 


0 


0 


(11,4) 


(11,5) 


0 


0 


(11,8) 


(11,9) 


0 


0 


(12,2) 


(12,3) 


0 


0 


(12,6) 


(12,7) 


0 


0 


(12,10) 



(11,11) 
o 



(12,12) 



(3.204) 

în care 

(l,l) = 60480 + 302408! + 241928 2 + 211688 3 + 190088 4 
(4,1) = L(30240 + 100808! + 70568 2 + 60488 3 + 54728 4 ) 
(5,1) = (6,2) = -(2,2) = -(5,5) = -(6,6) = -(1,1) 
(8,1) = L(30240 + 201608! + 1717368 2 + 151208 3 + 135368 4 ) 

(9.1) = -L(30240 + 151208! + 120968 2 + 10584S 3 + 95048 4 ) 
(11,1) = (10,2) = (12,2) = -(9,5) = -(11,5) = -(10,6) = -(l2,6) = (9,l) 

(3.2) =(5,4) =-(6,3) =-(4,1) 
(7,2)= (8,5)= -(7,6)= -(8,1) 
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(3,3) = (4,4) = L 2 (20160 + 50408! + 26888 2 + 20168 3 + 17288 4 ) 

(7,3) = (8,4) = L 2 (10080 + 50408! + 43688 2 + 40328 3 + 37448 4 ) 

(10,3) = L 2 (15120 + 50408! + 35288 2 + 30248 3 + 27368 4 ) 
(12,3) = -(9,4) = -(11,4) = (10,3) 

(7,7) = (8,8)= 1/(20160 + 151208! +127688 2 +110888 3 +97928 4 ) 

(10,7) = V (15120 + 100808! + 85688 2 + 75608 3 + 67688 4 ) 

(12,7) = -(9,8)= -(11,8)= (10,7) 

(9,9) = L 2 (15120 + 75608! + 60488 2 + 52928 3 + 47258 4 ) 

(11,9) = (10,10) = (12,10) = (ll,ll) = (12,12) = (9,9) 



[K s ] = 



Matricea de rigiditate la forfecare: 
0 

0 0 
0 0 0 
0 0 0 0 
0 0 0 0 0 
KLGA . 0 0 0 0 0 0 
120 0 0 0 0 0 0 0 

00000000 



00000000 
00000000 



(9,9) 

0 (10,10) 
00000000 (11,9) 0 

00000000 0 (12,10) 



sim 



(11,11) 

0 (12,12) 



(3.205) 
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unde 

(9,9) = (10,10) = (40 + 10 ttl + 4a 2 ) 
(11,9) = (12,10) = (20 + lOcti + 6<x 2 ) 
(ll,ll) = (12,12) = (40 + 30 ai + 24a 2 ) 

• Matricea de masă de translaţie este: 



c TJ 362880 



(1,1) 




















0 


(2,2) 


















0 


(3,2) 


(3,3) 














sim 


(4,1) 


0 


0 


(4,4) 














(5,1) 


0 


0 


(5,4) 


(5,5) 












0 


(6,2) 


(6,3) 


0 


0 


(6,6) 










0 


(7,2) 


(7,3) 


0 


0 


(7,6) 


(7,7) 








(8,1) 


0 


0 


(8,4) 


(8,5) 


0 


0 


(8,8) 






(9,1) 


0 


0 


(9,4) 


(9,5) 


0 


0 


(9,8) 


(9,9) 




0 


(10,2) 


(10,3) 


0 


0 


(10,6) 


(10,7) 


0 


0 


(10,10) 


(11,1) 


0 


0 


(11,4) 


(11,5) 


0 


0 


(11,8) 


(11,9) 


0 


0 


(12,2) 


(12,3) 


0 


0 


(12,6) 


(12,7) 


0 


0 


(12,10) 



(11,11) 

0 (12,12) 



(3.206) 



unde 

(l,l) = (2,2) = 134784 + 31104a! + 10944a 2 

(4,1) = -(9,1) = -(3,2) = -(10,2) = L(19008 + 6048 ai + 2448a 2 ) 

(5,1) = (6,2) = 46656 + 23328a! + 13248a 2 

(ll,l) = (7,2) = (12,2) = -(8,1) = L(ll232 + 5184 ai + 2736a 2 ) 

(3,3) = L 2 (3456 + 1296a! + 576a 2 ) 
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(10,3) = (4,4) = -(9,4) = (9,9) = (l0,10) = (3,3) 

(6,3) = -(5,4) = (9,5) = (10,6) = -L(ll232 + 6048a! + 3600a 2 ) 

(7,3) = -L 2 (2592 + 1296 ttl + 720a 2 ) 

((12,3) = (8,4) = -(11,4) = (10,7) = -(9,8) = (11,9) = (12,10) = (7,3)) 

(5,5) = (6,6) = 134784 + 103680a! + 83520a 2 

(8,5) = -(11,5) = -(7,6) = -(12,6) = -L(l9008 + 12960a! + 9360a 2 ) 

(7,7) = L 2 (3456 + 2160a! + 1440a 2 ) 

(12,7) = (8,8) = -(11,8) = (11,11) = (12,12) = (7,7) 



• Matricea de masă de rotaţie este de forma: 



[m r ] = 



Ml, 



362880L 



(u) 














0 


(2,2) 












0 


(3,2) 


(3,3) 








sim 


(44) 


0 


0 


(4,4) 








(54) 


0 


0 


(5,4) (5,5) 








0 


(6,2) 


(6,3) 


0 0 


(6,6) 






0 


(7,2) 


(7,3) 


0 0 


(7,6) (7,7) 






(8,1) 


0 


0 


(8,4) (8,5) 


0 0 


(8,8) 




(9,1) 


0 


0 


(9,4) (9,5) 


0 0 


(9,8) (9,9) 




0 


(10,2) 


(10,3) 


0 0 


(10,6) (10,7) 


0 0 


(10,10) 


(IU) 


0 


0 


(11,4) (11,5) 


0 0 


(11,8) (11,9) 


0 


0 


(12,2) 


(12,3) 


0 0 


(12,6) (12,7) 


0 0 


(12,10) 



(11,11) 

o (12,12) 



(3.207) 
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unde 

(l,l) = 435456 + 2177288! +1244168 2 + 777608 3 + 518408 4 
(4,1) = L(36288 + 362888! + 259208 2 + 181448 3 + 129608 4 ) 
(5,1) = (6,2) = -(2,2) = -(5,5) = -(6,6) = -(1,1) 
(8,1) = L(36288 - 103688 2 - 129608 3 - 129608 4 ) 

(9.1) = -L(217728 + 1088648! + 622088 2 + 388808 3 + 259208 4 ) 
(11,1) = (10,2) = (12,1) = -(9,5) = -(11,5) = -(10,6) = -(12,6) = (9,1) 

(3.2) =-(6,3) = (5,4) = -(4,l) 

(7.2) =(8,5) =-(7,6) =-(8,1) 

(3.3) = (4,4) = L 2 (48684 + 120968! + 69128 2 + 47528 3 + 3456S 4 ) 
(7,3) = (8,4) = -L 2 (12096 + 60488j + 51848 2 + 47528 3 + 43208 4 ) 
(10,3) = V (18144 + 181448J + 129608 2 + 90728 3 + 64808 4 ) 
(12,3) = (9,4) = -(11,4)= (10,3) 

(7,7) = (8,8) = L 2 (48384 + 362888! + 311048 2 + 280808 3 + 259208 4 ) 
(10,7) = V (18144 - 51848 2 - 64808 3 - 64808 4 ) 
(12,7) = -(9,8)= -(11,8)= (10,7) 

(9,9) = L 2 (108864 + 544328! + 311048 2 + 194408 3 + 129608 4 ) 
(11,9) = (10,10) = (12,10) = (ll,ll) = (12,12) = (9,9) 
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• Matricea giroscopică este: 
o 

(2,1) 0 

(3,l) 0 0 anti sim 

0 (4,2) (4,3) 0 

0 (5,2) (5,3) 0 0 

(6,l) 0 0 (6,4) (6,5) 0 

(7,1) 0 0 (7,4) (7,5) o o 

0 (8,2) (8,3) 0 0 (8,6) (8,7) 0 
0 (9,2) (9,3) 0 0 (9,6) (9,7) 0 0 

(10,1) o o (10,4) (10,5) o o (10,8) (10,9) o 

0 (ll,2) (ll,3) 0 0 (ll,6) )(ll,7)| 0 0 (ll,10) 0 

(12,1) o o (12,4) (12,5) o o (12,8) (12,9) o (12,11) o 

(3.208) 
unde 

(2,l)= 879012 + 4354568! +2488328 2 +1555208 3 +1036808 4 
(3,1) = -(72576 + 725768! + 518408 2 + 362888 3 + 259208 4 ) 

(6,1) = -(5,2) = -(6,5) = -(2,1) 

(7.1) = L(- 72576 + 207368 2 + 259208 3 + 259208 4 ) 

(10,1) = -L(435456 + 2177288! + 1244168 2 + 77760S 3 + 518408 4 ) 
(12,1) = -(9,2) = -(11,2) = -(10,5) = -(12,5) = (9,6) = (11,6) = (10,1) 

(4.2) =(5,3) =(6,4) = (3,1) 

(8.2) = -(7,5) = -(8,6) = (7,1) 

(4.3) = L 2 (96768 + 241928! + 13824S 2 + 95048 3 + 69128 4 ) 

(8,3) = -(7,4) = -L 2 (24192 + 120968! + 103688 2 + 95048 3 + 86408 4 ) 
(9,3) = -L 2 (36288 + 362888! + 259208 2 + 181448 3 + 12960S 4 ) 



[G] = 



ui, 



362880L 
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(11,3) =(10,4) =(12,4) = (9,3) 

(8,7) = 1^(96768 + 725768! +622088 2 + 561608 3 + 518408 4 ) 
(9,7) = L 2 (- 36288 + 103688 2 + 129608 3 + 129608 4 ) 
(11,7) =(10,8) =(12,8) = (9,7) 

(10,9) = L 2 (217728 + 1088648! + 622088 2 + 388808 3 + 259208 4 ) 
(12,9)= -(11,10)= (12,11)= (10,9) 



• Matricea circulatorie a elementului: 



[K c ] = 



EL 



5040L 3 



0 

(24) 


0 










(34) 


0 


0 






anti sim 


0 


(4,2) 


(4,3) 


0 






0 


(5,2) 


(5,3) 


0 0 






(6,1) 


0 


0 


(6,4) (6,5) 


0 




(7,1) 


0 


0 


(7,4) (7,5) 


0 0 




0 


(8,2) 


(8,3) 


0 0 


(8,6) (8,7) 


0 


0 


(9,2) 


(9,3) 


0 0 


(9,6) (9,7) 


0 0 


(10,1) 


0 


0 


(10,4) (10,5) 


0 0 


(10,8) (10,9) 0 


0 


(H,2) 


(11,3) 


0 0 


(11,6) (11,7) 


0 0 (11,10) 


.(12,1) 


0 


0 


(12,4) (12,5) 


0 0 


(12,8) (12,9) 0 



(12,11) o 



(3.209) 
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unde 

(2,1) = -(60480 + 302408! + 241928 2 + 211688 3 + 190088 4 ) 
(3,1) = L(30240 + 100808J + 70568 2 + 60488 3 + 54728 4 ) 
(6,1) = -(5,2) = -(6,5) = -(2,1) 

(7.1) = 30240 + 201608! + 171368 2 + 151208 3 + 135368 4 
(10,1) = L(30240 + 151208! + 120968 2 + 105848 3 + 95048 4 ) 
(12,1) = -(9,2) = -(11,2) = -(10,5) = -(12,5) = (9,6) = (11,6) = (10,1) 

(4.2) = (5,3) =(6,4) =(3,1) 

(8.2) = -(7,5) = -(8,6) = (7,l) 

(4.3) = -L 2 (20160 + 50408! + 26888 2 + 20168 3 + 17288 4 ) 

(8,3) = -(7,4) = -L 2 (10080 + 50408j + 43688 2 + 40328 3 + 37448 4 ) 
(9,3) = V (15120 + 50408! + 35288 2 + 30248 3 + 27368 4 ) 
(11,3) =(10,4) =(12,4) = (9,3) 

(8,7) = -L 2 (20160 + 151208i + 127688 2 + 110888 3 + 97928 4 ) 
(9,7) = L 2 (15120 + 100808! + 85688 2 + 75608 3 + 67688 4 ) 
(11,7) =(10,8) =(12,8) = (9,7) 

(10,9) = -L 2 (15120 + 7560&! + 60488 2 + 52928 3 + 47528 4 ) 
(12,9)= -(11,10)= (12,11)= (10,9) 
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• Matricea de masă [M T 1 are forma: 



anti sim 



(10,9) 

o (11,10) 
(12,9) o (12,11) o 

(3.210) 
unde 

(2,1) = (6,5) = 269568 + 62208a! + 21888a 2 

(3,1) = (10,1) = (4,2) = -(9,2) = -L(38016 + 12096a! + 4896a 2 ) 

(6,1) = -(5,2) = 93312 + 46656a x + 26496a 2 

(7,1) = (12,1) = (8,2) = -(11,2) = L(22464 + 10368 ai + 5472a 2 ) 

(4,3) = -(9,3) = -(10,4) = (10,9) = L 2 (6912 + 2592 ai + 1152a 2 ) 

(5,3) = (6,4) = -(10,5) = (9,6) = L(22464 + 12096a! + 7200a 2 ) 

(8.3) = (12,9) = -V (5184 + 2592 ttl + 1440a 2 ) 
(11,3) = -(8,3) 

(7.4) = (12,4) = (9,7) = (10,8) = (11,10) = (11,3) 

(7.5) = (12,5) = (8,6) = -(11,6) = L(38016 + 25920a! + 18720a 2 ) 

(8,7) = -(ll,7) = -(12,8) = (l2,ll) = L 2 (6912 + 4320a! + 2880a 2 ) 



[m t ]= 



362880 



A 
U 








n 






A 

u 


A 

u 


n 
U 


(d 1\ 


(d *\ 


Q 


[5,2) 


(5,3; 


(6,1) 


0 


0 


(74) 


0 


0 


0 


(8,2) 


(8,3) 


0 


(9,2) 


(9,3) 


(10,1) 


0 


0 


0 


(H,2) 


(11,3) 


(12,1) 


0 


0 



o 
o 

(6,4) 
(7,4) 

0 

0 



o 

(6,5) 

(7,5) 
0 

o 



(10,4) (10,5) 

o o 
(12,4) (12,5) 



o 
o 

(8,6) 

(9,6) 
0 



0 

(8,7) 
(9,7) 
0 



(11,6) (11,7) 
0 o 



o 

(10,8) 
o 

(12,8) 
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• Vectorul forţelor de dezechilibru masic: 

{Qc S }=^^{Q1 Q2 Q3 Q12} T (3.211) 

unde 

Ql = e i (1764 + 420a! + 156a 2 )+ e i (756 + 336a! + 156a 2 ) 

Q2 = £, (1764 + 420a! + 156a 2 )+ ^ (756 + 336aj + 156a 2 ) 

Q3 = -Li;, (252 + 84a! + 36a 2 )- (l68 + 84a! + 48a 2 ) 

Q4 = Le, (252 + 84ai + 36a 2 )+ Lej (168 + 84a! + 48a 2 ) 

Q5 = Ej (756 + 420 ttl + 264a 2 )+ e i (1764 + 1344 ai + 1080a 2 ) 

Q6 = Ci(756 + 420a! +264a 2 )+Cj(l764 + 1344a! +1080a 2 ) 

Q7 = L£. (168 + 84aj + 48a 2 )+ (252 + 168a! + 12 <*<*2 ) 

Q8 = -LEj (168 + 84ai + 48a 2 )- Lej (252 + 168a! + 12 °a 2 ) 

Q9=-Q4 
Q10=Q3 
Qll=-Q8 
Q12=Q7 

în care 

e şi £ reprezintă coordonatele centrului de masă excentric în planele 
perpendiculare pe nodurile i şi j ale elementului. 
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3.1.2.2. Elementul CON 8 

In lucrarea [58], Genta şi Gugliotta prezintă un element de grindă 
conic axial simetric, având secţiunea transversală circulară sau inelară, bazat pe 
teoria de bară Timoshenko. Acest element are două grade de libertate complexe 
pe flecare nod şi autorii îl consideră mai eficient şi preferabil de utilizat. Cu 
toate că notaţia complexă folosită aici devine foarte greoaie, dacă întregul sistem 
nu este axial simetric, metoda descrisă poate fi uşor extinsă la rotori axiali 
simetrici pe suporţi axial nesimetrici, dublând numărul ecuaţiilor şi rearanjând 
matricele. 

La elementul conic CON 8, bazat pe [58] deformaţiile de 
încovoiere şi de forfecare sunt cuplate şi astfel sunt incluse în analiza arborelui 
şi gradele de libertate suplimentare. Aceasta este baza în descrierea elementului 
conic CON 8. El are doar patru grade de libertate pe nod, deplasarea laterală şi 
rotirea secţiunii transversale în două plane perpendiculare şi astfel, acest element 
se poate implementa direct într-un program de calcul cu elemente finite 
cilindrice, ceea ce constituie un real avantaj. De fapt elementul CON 8 poate 
înlocui elementul cilindric, acesta din urmă devenind un caz particular. 

Geometria elementului rezultă din expresii polinomiale. Acestea 
sunt utilizate pentru determinarea matricelor elementului şi includ efectele 
inerţiei la rotaţie şi la translaţie, a cuplului giroscopic, a deformaţiilor de 
forfecare şi de încovoiere. Matricea de masă şi matricea giroscopică sunt matrice 
coerente de ordin 8x8. Matricea de rigiditate pentru elementul CON 8 este 
aceeaşi cu matricea obţinută după condensarea elementului CON 12. Utilizând o 
aranjare convenabilă a gradelor de libertate în vectorul coordonatelor nodale, se 
obţin submatrici simple de ordin 4x4 [140]. 

Funcţiile de formă, bazate pe funcţiile de formă pentru elementul de 
grindă Timoshenko, conţin constanta de forfecare, care este proporţională cu 
unghiul de forfecare, acesta fiind presupus a fi constant în lungul elementului. 
Constanta de forfecare K este funcţie de modulul de elasticitate longitudinal, de 
modulul de elasticitate transversal, de caracteristicile geometrice ale elementului 
şi de factorul de forfecare %. Factorul de forfecare este funcţie de coeficientul lui 
Poisson u şi de raza interioară r şi de raza exterioară R a arborelui. 

Elementul conic CON 8 este prezentat în Figura 3.14. Acesta este 
un element de secţiune circulară sau inelară având grosimea liniar variabilă, 
lungimea L şi nodurile de la capete notate cu i şi j. 

Cu notaţia arătată în figură, expresiile pentru aria A şi momentele 
de inerţie transversal I t şi polar I p ale secţiunii transversale sunt: 

A(^) = A i (l + a^ + p i ^ 2 ), (3.212) 
I t (0 = I ti(l + ^ + P2^ 2 +Y 2 4 3 +« 4 ), (3.213) 
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(3.214) 




Fig.3.14 



unde coordonata adimensională £ este de forma: L = — . 

* L 



(3.215) 



Coeficienţii a 1 ,p i ,a 2 ,P 2 ,y 2 ,8 2 depind de geometria elementului şi au 
următoarele expresii: 



a, = 



27c(Ri AR -r, Ar) 



(3.216) 



_ rc(AR 2 -Ar 2 ) 
Pl " A : 



(3.217) 



a, = 



rc^AR-r^Ar) 



(3.218) 
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P 2 = 



y 2 = 



8,= 



37r(Rj 2 AR 2 -r ; 2 Ar 2 ) 



2I ti 



^(RjAR'-rjAr 3 ) 
I ti 

7i(AR 4 -Ar 4 ) 
4L. 



(3.219) 



(3.220) 



(3.221) 



unde 



AR = Rj -R;, 



(3.222) 



Ar = r j — T; , 

în care R şi r reprezintă raza exterioară, respectiv raza interioară a secţiunii 
transversale. 

Matricele elementului 

Deplasările elementului pot fi exprimate, în notaţii complexe astfel: 



{**}= 



f v + iwl 


"N" 




_N_ 



{ti 



(3.223) 



Vectorul {q s }reprezintă coordonatele generalizate complexe şi este de forma: 



{5 S }= 



v j +iw j 
IVj - i<Pj J 



(3.224) 



Funcţiile de formă |_Nj şi |_nJ sunt aceleaşi cu cele folosite pentru elementul de 
grindă cilindric tip Timoshenko: 

LnJ=Ln, n 2 n 3 n 4 J, 
LnJ=Ln 1 n 2 n 3 n 4 J. 
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Expresiile acestora sunt date în relaţiile (3.122) şi (3.123). Ca şi în aceste relaţii, 
valoarea constantei K este 



K = 



12 EI 



şi poate fi calculată utilizând caracteristicile geometrice A, I şi % ale secţiunii 
transversale de la mijlocul elementului. Coeficientul de forfecare % este dat de 
relaţiile(3.59),(3.60) [48]. 

Submatricea de rigiditate [k ], submatricea de masă [m ] şi 

c 

submatricea giroscopică [g ] se definesc astfel: 



[ kS ]=!K4 LNjr 4 LsJ l d5+GL î xA 



H--LnJ 



-iTr 



LnJ--LnJ 

J d£ L J 



(3.225) 



k]=[m/] + [m/]=p L f^^^^ 



(3.226) 



[g S ] = pL) I pLNfLN^ = 2[m R s ], 



(3.227) 



unde p este densitatea materialului. 

In această situaţie submatricea giroscopică este simetrică, ceea ce 
constituie un avantaj. 

Vectorul forţelor pentru distribuţia dezechilibrului static e(%) şi 

pentru distribuţia dezechilibrului dinamic 0(^) este 



f L L 

{ f s }= (o 2 1 pL J A|_N J t ed£ - { I [_N J 0d£ 



(3.228) 



unde £i reprezintă viteza unghiulară a arborelui. 
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Introducând funcţiile de formă în expresiile de mai sus şi rezolvând 
integralele, se găsesc următoarele expresii pentru submatricele elementului 
conic CON 8: 



EI 



105L 3 (l + K) 2 



G/AjK/ 



12L(l + K) 2 



unde 

k, = 1260 + 630a 2 + 5O40 2 + 441y 2 + 3968 2 , 



X 


Lk 2 


"k, 


Lk 3 




L 2 k 4 


-Lk 2 


Lk 5 








-Lk 2 




sim 




L 2 k 6 


Lk 8 


-k 7 


Lk 8 " 




L 2 k 9 


-Lk 8 


L 2 k 9 








Lk 8 




sim 




L 2 k 9 





(3.229) 



k 2 =630 + 210a 2 +147p 2 +126y 2 +1148 2 -k(l05a 2 +1050, + 94 5y 2 +848 2 ) 
k 3 = 630 + 420a 2 + 357p 2 + 315y 2 + 2828 2 + k(l05a 2 + 105p 2 + 94 • 5y 2 + 848 2 ) 

k 4 = 420 + 210k + 105k 2 + a 2 (l05 + 52 • 5k 2 )+ P 2 (56 - 35k + 35k 2 )+ 
+ y 2 (42 - 42k + 26 • 25k 2 )+ 8 2 (36 - 42k + 21k 2 ), 

k 5 = 210 - 210k - 105k 2 - a 2 (l05 - 105k - 52 • 5k 2 )+ P 2 (91 - 70k - 35k 2 )+ 
+ y 2 (84 - 52 • 5k - 26 • 25k 2 )+ 8 2 (78 - 42k - 21k 2 ), 

k 6 = 420 + 210k + 105k 2 + a 2 (315 + 210k + 52 • 5k 2 )+ p 2 (266 + 175k + 35k 2 )+ 
+ y 2 (231 + 147k + 26 • 25k 2 )+ 8 2 (204 + 126k + 21k 2 ), 

k 7 =(12 + 6^ +4pJ, 

k 8 = 6 + 3^ +2p 15 

k 9 =3 + 1-5^ +p t . 
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[m T s ]= 



PA;L 



1260(l + K) 2 



m, Lm 2 m 3 -Lm 4 



L 2 m 



'3 



sim 



s —6 -Lm 7 
m 8 - Lm 9 
L 2 m 10 



(3.230) 



unde 



m 



m 



m 



m 



m 



m 



m 



m 



m 



1 = (468 + 882k + 420k 2 )+ ttl (l08 + 210k + 105k 2 )+ (38 + 78k + 42k 2 ), 
= (66 + 115- 5k + 52-5k 2 )+ot 1 (21 + 40-5k + 21k 2 )+p i (8-5 + 18k + 10-5k 2 ) 
= (l62 + 378k + 210k 2 )+ ttl (81 + 189k + 105k 2 )+ (46 + lllk + 63k 2 ), 
= (39 + 94. 5k + 52-5k 2 )+a 1 (l8 + 40-5k + 21k 2 )+p i (9-5 + 21k + 10-5k 2 ) 
= (l2 + 21k + 105k 2 )+a 1 (4.5 + 9k + 5.25k 2 )+p i (2 + 45k + 3k 2 ), 
= (39 + 94-5k + 52-5k 2 )+a 1 (21 + 54k + 31-5k 2 )+p i (l2-5 + 34 5k + 21k 2 ) 
7 =(9 + 21k + 10-5k 2 )+a 1 (4-5 + 10-5k + 5-25k 2 )+p i (2-5 + 6k + 3k 2 ), 
= (486 + 882k + 420k 2 )+ a, (360 + 672k + 315k 2 )+ ^ (290 + 540k + 252k 2 ) 
= (66 + 115- 5k + 52 5k 2 )+a 1 (45 + 75k + 31-5k 2 )+p i (32-5 + 52-5k + 21k 2 ) 
= (l2 + 21k + 10.5k 2 )+a 1 (7.5 + 12k + 5.25k 2 )+p 1 (5 + 7.5k + 3k 2 ). 



[ m R S ] = 



Pli 



210L(l + K) 2 



Lm 12 


-m n 


Lm 13 


L 2 m 14 


-Lm 12 


-L 2 m 15 






-Lm 13 


sim 




L 2 m 16 



(3.231) 



unde 
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m n = 252 + 126a 2 + 72p 2 + 45y 2 + 308 2 , 

m 12 = 21 - 105k + a 2 (21 - 42k)+ p 2 (15 - 21k) + y 2 (10 • 5 - 12k) + 8 2 (7 • 5 - 7 • 5k) 
m 13 = 21 - 105k - 63a 2 k - p 2 (6 + 42k)- y 2 (7 • 5 + 30k)- 8 2 (7 • 5 + 22 • 5k), 

m 14 =28 + 35k + 70k 2 +a 2 (7-7k + 17-5k 2 )+P 2 (4-7k + 7k 2 )+ 
+ y 2 (2-75-5k + 3-5k 2 )+8 2 (2-3 5k + 2k 2 ), 

m 15 =7 + 35k-35k 2 +a 2 (3 5 + 17 5k-17-5k 2 )+p 2 (3 + 10-5k-10-5k 2 )+ 
+ y 2 (2-75 + 7k-7k 2 )+8 2 (2-5 + 5k-5k 2 ), 

m 16 = 28 + 35k + 70k 2 + a 2 (21 + 42k + 52 • 5k 2 )+ P 2 (l8 + 42k + 42k 2 )+ 
+ y 2 (l6 • 25 + 40k + 35k 2 )+ 8 2 (l5 + 37 • 5k + 30k 2 ). 



(3.232) 



Presupunând o distribuţie liniară a dezechilibrului, vectorul forţelor 
nodale devine: 



{F}= 



120(l + K) 



A,L[a} 



£ y +i8 r,jj 



(3.233) 



unde elementele matricelor [a] şi [b] sunt: 



a u =42 + 40K + 10a 1 (l + K)+p i 



26 



+ 4K , 



a 12 = 18 + 20K + a, (8 + 10K)+ ^ \ y + 6K 



a 21 =L 



6 + 5K + 2a 1 (l + K)+p 1 



+ K 
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a 22 - L 



4 + 5K + a 1 (2 + 3K)+p i 



8 



+ 2K 



a 31 =18 + 20K + 10a 1 (l + K)+p 1 



'44 ^ 
— + 6K 

V7 j 



a 32 =42 + 40K + a 1 (32 + 30K)+p i ^î^ ) + 24K 



a 41 =-L 



4 + 5K + 2a 1 (l + K)+p i 



- + K 
V7 J 



a 42 =-L 



6 + 5K + a 1 (4 + 3K)+p i 



— + 2K 



V 



b n =" 



A 60 + 24a 2 +12p 2 + Ş Y2+ Ş5 2 > 



V 



J 



b 12 =- 



' 120 90 ^ 

60 + 36a 2 + 24p 2 + — y 2 + — 8 2 



V 



V 



b 21 =L 



'-10 „ ( . 8 



10 + 40K + a 2 (-2 + 10K)+p 2 (-2 + 4K)+y 2 — - + 2K +8 



-1 + -K 

v 7 , 



b 22 =L 



/ ^2 A ( 25 20 A" 

-10 + 20K + a 2 (-8 + 10K)+P 2 (-6 + 6K)+y 2 + 4K + 8 2 + — K 

^7 ) \ 7 7 J_ 



b 31 = -b n 



^32 _ k 12 



b 41 =L 



-10 + 20K + a 2 (-2 + 10K)+6p 2 K + y. 



+ 4K +5 ^ + f K 
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b 42 =L 



10 + 40K + a 2 (12 + 30K)+ 0 2 (l2 + 24K)+ y . 



80 



+ 20K 

7 j 




Elementul conic de tip CON 8 poate fi folosit şi pentru rotorii axial 
nesimetrici, sau când prezenţa lagărelor introduce elemente nesimetrice în 
matricele de rigiditate şi de amortizare. Dacă se foloseşte o notaţie 
convenţională bazată pe coordonate reale şi gradele de libertate sunt descrise în 
următoarea ordine 

{u s }={ Vi \\f i Vj Wi -Cpi Wj -cpjf (3 234) 

atunci ecuaţia de mişcare se poate scrie astfel: 



K] [0] 

L M [m,U 



+ C0 



[0] [G]\ 
L"[G] [0]_ 



[c,] [cj\ 
[cj [c,] 



[k,] [k„J1 

K] [Kj 



+ (0 



[0] [c,Jl 
,-[c t ] [0]. 



Im^je 1 



(3.235) 



Când elementul are o simetrie axială, matricele cu indicele y şi z 
sunt egale şi coincid cu cele obţinute aici, iar matricele cu indicele yz dispar. 
Matricele asamblate nu au aceste proprietăţi. Matricea globală giroscopică este 
acum antisimetrică. 

Vectorul deplasărilor |u s } este uşor diferit de cel utilizat în mod 
normal: (-cp) a fost, de fapt, ales ca rotire în jurul axei y în loc de (cp). 

Concluzie:Utilizarea elementului finit conic de tip CON 8 în calculele de 
dinamica rotorilor, construit cu două grade de libertate complexe pe fiecare nod, 
este mai eficientă decât utilizarea altor variante de calcul, datorită aplicabilităţii 
mai generale. 

Elementele de arbore tronconic descrise în paragrafele 3.1.2.1. 
(elementul CON 12) şi 3.1.2.2. (elementul CON 8), deoarece ţin cont de inerţia 
la rotaţie şi de deformaţiile de forfecare prin intermediul constantei 



K = 12 



EI 



GA S L 



se numesc, în general, elemente tronconice de bară 



Timoshenko. Dacă în expresiile matricelor acestor elemente constanta K se 
egalează cu zero, adică dacă nu se ţine cont de forfecarea elementului, ci numai 
de inerţia la rotaţie, atunci avem de a face cu elemente tronconice de bară de tip 
Rayleigh. 

Dacă în expresiile matricelor elementelor descrise nu se ţine cont 
nici de inerţia la rotaţie, nici de forfecarea elementului, ci numai de încovoierea 
acestuia, atunci avem de a face cu elemente tronconice de bară de tip Bernoulli- 
Euler. 
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3.2. Discuri rigide 



în studiul mişcării laterale a rotorului, deplasarea oricărui punct, 
incluzând şi centrul geometric al discului, se defineşte prin două translaţii (v, w) 
şi două rotaţii (cp, \|/). 

în cele ce urmează rotorul se consideră axial simetric având axele 
ox, oy şi oz axe principale de inerţie, figura 3.15. 





Fig.3.15 

• Matricea de inerţie şi matricea giroscopică 



Neglijând orice dezechilibru, energia cinetică a discului axial 



simetric este: 



T d =im d (v 2 + w 2 )+i(j x o) x 2 + J y (D y 2 + J z to z 2 ) 



(3.236) 



unde 

Jx=J P şi J y =Jz=JT sunt momentele de inerţie masice principale ale discului faţă 
de sistemul de axe xyz fixat pe disc; 
J p =moment de inerţie masic polar al discului; 
J T =moment de inerţie masic transversal al discului; 
m d =masa discului. 
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Componentele vitezei unghiulare în sistemul de axe mobil xyz sunt: 







Q + cp\|/ 


«co y 


► = < 


(pcosG + \j/sinG ► 






vj/cosG - cpsinG 



(3.237) 



m care 



Q=viteza unghiulară cu care discul se roteşte în jurul axei x. 

înlocuind aceste componente în expresia energiei cinetice a 
discului, se obţine: 



T d =im d (v 2 + w 2 ) + ij p (Q + (p M /) 2 + 



+ — J T [(cpcosG + \j/sinG) 2 + (vj/cosG - cpsinG) 2 ] 



(3.238) 



T d = im d (v 2 + w 2 )+ 1 J p Q 2 + J p Qcj>y + \ J p cj> V 2 + |j T (q> 2 + ij/ 2 ) 



translaţia 



efect 

energie giroscopic 
rotaţie. 



rotaţie 



(3.239) 



Aplicând ecuaţiile Lagrange, se obţine: 



d_ 
dt 



(3.240) 



=0 



dt 



5T l 



K dy ) d\|/ 



= J T y-QJ p cp = M z , 



(3.241) 
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d_ 
dt 



dw 

=o 



= m H w = T, 



(3.242) 



d 5T d _ dT d 
dt <p)J d(-(p) 



= J T (-(p)-QJ P \i/ = -M 



(3.243) 



în formă matriceală 



m d 0 
0 J T 



0 
0 



0 
0 



0 0 m d 0 
0 0 0 J T 



v 

¥ 

w 



+ Q 



0 


0 


0 


0" 


V 






0 


0 


0 


Jp 


¥ 




M, 


0 


0 


0 


0 


w 


> = < 


> 


0 


-Jp 


0 


0 


-9. 




-M y 



(3.244) 



Introducând vectorii deplasare de ordin 2x1 din planele X-Y, respectiv X-Z 

{<}-[ 



v 

LVJ 



(3.245) 



K1=| 



w 



şi vectorii corespunzători sarcinilor ce acţionează asupra discului în cele două 
plane 



M 



(3.246) 



M 



y J 



105 



Dana-Codruţa Vişan Modelarea dinamică a sistemelor rotor-lagăre 



ecuaţia matriceală devine: 




(3.247) 



unde membrul drept al forţelor care acţionează asupra discului include 
dezechilibru masic, dezechilibru datorat lipsei de aliniere, forţe de legătură şi 
alte efecte. 

Submatricea 



kl= 



m 
0 



0 

J f 



(3.248) 



reprezintă submatricea de inerţie pentru un disc axial simetric, iar submatricea 



0 0 
0 J P 



(3.249) 



reprezintă submatricea giroscopică pentru un disc axial simetric. 



Vectorii forţelor de dezechilibru 



Pentru analiza răspunsului unui rotor, este important să se 
determine vectorii dezechilibrului masic şi cei ai dezechilibrului datorat montării 
sub un unghi a discului. 



• Vectorul dezechilibrului masic 

In figura 3.16 centrul de masă G al discului este diferit de centrul 
geometric C al său, aşa încât CG se numeşte excentricitate şi se notează cu a. Se 
consideră sistemul de axe fix OXYZ, cu axa OX în lungul liniei lagărelor şi 
sistemul de axe mobil, legat de disc, Cxyz. In acest sistem de axe mobil, la 
timpul t=0, se notează a c =a cos a, a s =a sin a, unde a este unghiul făcut de axa 
Cy cu excentricitatea CG. 
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Forţa centrifugă de dezechilibru care acţionează în centrul de masă 
G, în lungul lui CG este m d aQ . Componentele acesteia în sistemul de axe 

1 1 

rotitor Cxyz sunt: m d a c Q şi m d a s Q , iar componentele în sistemul de axe fix 
sunt: 

f y = m d a c Q 2 cosOt - m d a s Q 2 sinQt , 

(3.254) 

f z = m d a s Q 2 cosQt + m d a c Q 2 sinQt . 



Cu ajutorul acestor componente se formează vectorul forţelor de 

dezechilibru. 




Fig.3.16 

Vectorul forţelor de dezechilibru (neglijând înclinarea discului) 

este: 



{,-}= 



a c Q 2 




-a s Q 2 


0 




0 


►cosQt + m d < 


a s Q 2 




a c Q 2 > 


0 




0 



sinDt. 



(3.251) 
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• Vectorul datorat înclinării permanente a discului 

Un disc montat oblic pe arbore introduce momente giroscopice 
active care fac ca forţele din disc să-şi schimbe intensitatea. 

Se consideră un disc rigid care nu este perpendicular pe arbore, ca 
în figura 3.17. Mărimea vectorială x reprezintă abaterea unghiulară statică dintre 
axa principală de inerţie a discului şi axa arborelui, în poziţia iniţială, adică la 



componentelor unui vector rotitor, x, normal pe axa de rotaţie şi pe linia 
unghiului maxim de înclinare. Vectorul x face un unghi p cu axa rotitoare z, 
astfel încât vectorul discului înclinat are componentele 



t=0. 



Discul montat oblic pe arbore poate fi definit cu ajutorul 



x c = x cos (3 



(3.252) 



x s = x sinp 



(3.253) 



în sistemul de axe mobil xyz. 



4 



y 




Fig.3.17 



Componentele deplasărilor unghiulare devin: 




(3.254) 



unde (p şi \|/ sunt rotiri elastice. 
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Vitezele şi acceleraţiile unghiulare corespunzătoare sunt: 
d y = (j> - Qxcos(Qt + p) 

d z = \j/ - Qxsin(Qt + p) 



(3.255) 



d y = (p + Q 2 xsin(Qt + p) 
d z = \|> - Q 2 xcos(Qt + p) 



(3.256) 



Utilizând ecuaţia a doua şi a patra a lui Lagrange şi înlocuind în 
acestea pe 9 cu a y şi pe \\f cu a z , se obţine: 

f J T d z +OJ P (-d y )=M z 

' J T (-d y )-QJ P d z =-M y 



(3.257) 



Utilizând expresiile vitezelor şi acceleraţiilor unghiulare stabilite 
anterior, se obţin componentele cuplului giroscopic: 

M z = J T \|> - J p Q(p + (j p - J T )Q 2 TCOs(Qt + p) 
-M y = J T (-(p)-J p Q\j/ + (jp-J T )Q 2 Tsiii(Qt + p) 



(3.258) 



Vectorul momentelor giroscopice active care acţionează asupra 



discului este: 
f-M. 



M 



y J 



= (J P -J T )Q : 



( 



TCOSP 

xsiiip 



cosQt + 



-xsiiip 

TCOSP 



sinQt 



(3.259) 



Vectorul forţelor şi momentelor de dezechilibru este: 
{f d }=Q 





m d a £ 




- m d a s 


Q 2 < 


-(j P - JtK 


>cosQt + Q 2 < 


(Jp-JtK 




m d a s 








— (jp — Jt^s. 




— (jp — Jt)^c. 



sinQt. (3.260) 



Diferenţa dintre momentele de inerţie ale secţiunii J P -J T poate fi 
pozitivă (disc subţire), sau poate fi negativă (disc gros). Când Jp-Ji>0, 
momentul aplicat are efect de redresare, el acţionează asupra discului în direcţie 
opusă înclinării iniţiale. Când Jp-Jx<0, momentul aplicat creşte dezechilibrul 
unghiular acţionând asupra discului gros în continuare în direcţia înclinării 
iniţiale. 
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3.3. Lagăre 

Maşinile cu rotor sunt întotdeauna rezemate pe lagăre. Acestea fac 
parte din sistemul rotor-lagăre şi se modelează în funcţie de tipul reazemului. Se 
prezintă în cele ce urmează modelarea lagărelor de rostogolire rezemate elastic 
şi modelarea lagărelor cu alunecare. 



3.3.1. Lagăre de rostogolire rezemate elastic 

Lagărele radiale sunt modelate de obicei luând în considerare 
rigiditatea la translaţie şi coeficienţii de amortizare. In cele ce urmează nu sunt 
luate în considerare efectele inerţiei. 

Matricea de rigiditate pentru un lagăr radial 

Relaţia între forţe şi deplasări pentru un lagăr radial, în sistem de 
axe yz, se poate scrie: 



k yy k yz 

^zy ^zz 



V 
W 



v 



w 



(3.261) 



unde elementele matricei [k ] se numesc coeficienţi de influenţă, astfel încât ky 
reprezintă forţa elastică pe direcţia "i" datorată deplasării unitare de pe direcţia 

"j" 

Se consideră sistemul de axe y'z', ca în figura 3.18, rotit cu un 
unghi a faţă de sistemul de axe yz. 

Vectorul deplasărilor în noul sistem de axe se poate scrie: 



f ^^ 
V 




cosa sin a 




* > 

V 








-sin a cosa 


< 






' w' 

V J 






w 





-N 



v 
w 



(3.262) 



Vectorul forţelor în noul sistem de axe este: 



cosa sin a 
-sina cosa 



-N 



(3.263) 
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Fig.3.18 

Renunţând la indicele "b" de la notaţia [k ], noua relaţie între forţe 
şi deplasări devine: 



-N 



= [tM 



V 

w 



-frHrr 



W 



(3.264) 



1 X 

Pentru că [TŢ =[T] , matricea de rigiditate transformată este: 



M-frMrr- 



c s 


k yy 


k yz 


C 


-s 


-s c 






s 


c 



(3.265) 



unde c=cos a şi s=sin a. 

Matricea de rigiditate originală, nesimetrică, se poate descompune 
într-o sumă de două matrici, una simetrică [k s ] şi cealaltă antisimetrică [k a ]: 



k yy k yz 



k s k 

s zz 



+ 



0 k. 



" k a 0 



(3.266) 



[K] 



unde 



k a = 



2 

k yz -k zy 



(3.267) 
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Matricea de rigiditate transformată devine: 



[k']=[T]([k,] + [kjTr =["/] + ["•'} (3-268) 



unde 



k yy c 2 +k zz s 2 +2k s cs (k^-k^c + k^-s 2 )" 
(k^-k^c+k^-s 2 ) k^s'+k.c 2 -2k.es 



(3.269) 



şi [k a ']=[k a ]. 



Aceste rezultate implică următoarele: 



(3.270) 



Componenta antisimetrică [k a ] este independentă de matricea de rotaţie 
[T], deci şi de unghiul a. 

2k 

Pentru tg2a* =- * — , componenta simetrică [k s ] devine diagonală 



k yy k zz 



şi unghiurile a* si a* + 90° definesc direcţiile principale de 
elasticitate. 



Matricea simetrică [k s ] se poate diagonaliza prin transformarea de 
coordonate [T], în timp ce matricea antisimetrică [k a ] rămâne neschimbată. 
Aceste observaţii implică faptul că matricea de rigiditate nesimetrică se poate 
transforma într-o matrice în care elementele de pe diagonala secundară sunt 
antisimetrice. 



In cazul în care elementele diagonale sunt identice, apar două 
cazuri particulare interesante: 



Când elementele de pe diagonala secundară sunt antisimetrice, matricea 
transformată este aceeaşi cu cea originală: 



c 

-s 



-k 



c 

-s 



-k 



prin urmare lagărul este izotrop . 
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Când elementele de pe diagonala secundară sunt identice 



[ k -]= 




\a l +k 2 sin2a 



k 2 cos 2a 



k 2 cos 2a 
-k 2 sin2a 



şi pentru a* = 45° 



[k'] = 



k 1+ k ; 
0 



0 

k,-k : 



astfel, lagărul este ortotrop , având axele principale de elasticitate la 45° şi 
la 135°. 



Concluzii 



1 . Dacă | k yz | ^ | k zy | - lagăr anizotrop; 

2. Dacă k yz =-k zy şi k yy =k zz - lagăr izotrop; 

3. Dacă k yz =k zy şi k yy ^k zz - lagăr ortotrop, având direcţiile principale 
definite de unghiurile a* si a* + 90° ; 

4. Dacă k yz =k zy şi k yy =k zz - lagăr ortotrop, având direcţiile principale la 
45° şi la 135°; 

5. Dacă k yz =k zy =0 şi k yy ^k zz - lagăr ortotrop având direcţiile principale pe 
orizontală şi verticală; 

6. Dacă k yz =k zy =0 şi k yy =k zz - lagăr izotrop; 

7. Cuplajul antisimetric în componenta [k a ], adică k yz ^k zy , este cauza 
principală a instabilităţii într-un sistem anizotrop. 



3.3.2. Lagăre cu alunecare 

Proiectarea, construcţia şi exploatarea lagărelor cu alunecare 
necesită mai întâi o cunoaştere a ecuaţiilor, relaţiilor de calcul, proprietăţilor şi 
particularităţilor calitative ale funcţionării acestora. Marea majoritate a 
proprietăţilor, particularităţilor şi caracteristicilor de funcţionare ale lagărelor cu 
alunecare se pot deduce plecând de la studiul mişcării din pelicula de lubrifiant 
fluid, care desparte cele două suprafeţe în contact. 
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Coeficienţi dinamici ai lagărelor cu alunecare 

Se consideră că fusul se mişcă în cuzinet de-a lungul unei orbite, în 
jurul poziţiei de echilibru, figura 3.19. 

Zo 




orbita centrului 
fusului 



poziţia de 
echilibru 
staţi c 



turaţia 
-H_ fusului 



Fig.3.19 



Rezultanta forţei de reacţiune din filmul de lubrifiant are 
amplitudine şi direcţie variabilă (nu mai este verticală). Componentele sale F y şi 
F z sunt funcţii neliniare de deplasările y şi z ale centrului fusului şi de 
componentele vitezei acestuia y si z : 



F y =F y (y,z,y,z), 
F z =F z (y,z,y,z). 



(3.271) 



Pentru o amplitudine mică a mişcării, reacţiunea cuzinetului poate 
fi exprimată prin termenii de ordinul întâi dintr-o dezvoltare în serie Taylor în 
funcţie de componentele poziţiei de echilibru static: 



F y= F y„ +k yyy +k y Z Z + C yyy + V> 



(3.272) 



F z=Fz 0 +k zy y + k zz z + c zv y + c zz z. 



ZZ zv • 



114 



Capitolul 3. Modele cu elemente finite ale sistemelor rotor-lagăre 



Componentele forţei de reacţiune statică sunt: 



F v = W 



F z =0. 



(3.273) 



Cei opt coeficienţi ai componentelor forţelor liniarizate se 
calculează ca gradienţi corespunzători poziţiei de echilibru static 
( y = z = y = z): 



k. = 



yy 



k = 



zy 



C yy 



C zy = 



fdF^ 



V W J J 

z 



k = 



yz 



J 



k = 



v az 



C yz = 



C„ = 



az 



(3.274) 



Aceştia rezultă din soluţia ecuaţiei lubrifiantului. Liniarizarea 
forţelor de reacţiune din cuzinet are avantajul decuplării ecuaţiei mişcării 
rotorului de cea a lagărelor. Pe de altă parte, ecuaţiile rotorului trebuie să se 
integreze simultan cu cea a lubrifiantului. 

In formă matriceală, expresiile forţelor F y şi F z se pot scrie: 







fwl 






l*.J 


► = < 


0 





(3.275) 



Creşterile forţelor din filmul de lubrifiant, datorate mişcărilor mici 
în jurul poziţiei de echilibru static, se exprimă în funcţie de coeficienţii de 
influenţă ky şi în funcţie de coeficienţii de amortizare Cy: 
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AF 

y 

AF 



k yy k yz 

k zy k zz 



y+ 



c c 

yy yz 

C C 

zy zz 



>=[k b Fi+[c b } 



.(3.276) 



Cei opt coeficienţi liniarizaţi, de influenţă şi de amortizare, depind 
de condiţiile de funcţionare ale fusului, prin urmare depind de viteza de rotaţie. 



Coeficienţii de influenţă adimensionali se definesc după cum 



urmează: 



K yy K yz 

K zy K zz 



(3.277) 



unde C reprezintă jocul din lagăr, adică diferenţa între raza cuzinetului şi cea a 
fusului, C=Rc-Rf« 

Coeficienţii de amortizare adimensionali se definesc astfel: 



C C 

yy yz 

c c 

^ zy ^zz 



w L ' 



(3.278) 



unde Q reprezintă viteza unghiulară a fusului, exprimată în rad/sec. 



Pentru un set de parametri geometrici daţi şi vâscozitatea 
lubrifiantului cunoscută, aceşti opt coeficienţi sunt funcţie de numărul 
Sommerfeld, S, sau funcţie de excentricitatea fusului în cuzinet. Valorile pentru 
aceşti coeficienţi sunt daţi în cartea lui Someya [158]. 

Pelicula de lubrifiant fluid din lagăr, se modelează prin intermediul 
variaţiilor forţelor T y şi T z , obţinute în urma unei dezvoltări în serie Taylor în 
raport cu relaţiile pentru forţă, deplasare şi viteză, corespunzătoare configuraţiei 
fusului în lagăr la echilibru: 



yy 

b 



K yz 

b 



kzy 



b 

yy 

b 

zy 



(3.279) 



FT 



In formă expandată: 
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0 
u 




o 

U 




V 




C yy 


u 














0 


0 


0 


0 








0 


0 










< 










0 


k b 


0 




w 




c b 


0 




zz 










0 


0 


0 


0 




"9. 




0 


0 



yz 



0 

o 
o 
o 



c . * 

v 

V 




T 

A y 


< 


► = < 


> 


w 






1-9. 




-M,_ 



. (3.280) 



In formă partiţionată: 



[ K yy] [ k îzJ 


te} 




"fel 


te} 






'tel 


► + 


„[<] kl 


'tef 





(3.281) 



unde: 



0 



*5 
0 0 



(3.282) 



0 

0 0 



ij=y,z. 



(3.283) 



In aceste relaţii s-au considerat numai coeficienţii de translaţie. 
Efectele inerţiei lagăr- fundaţie se neglijează. 
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4. Răspunsul dinamic al rotorilor modelaţi cu elemente 
de bară Timoshenko cilindrice şi tronconice 

4.1. Introducerea elementului tronconic 
în programul DRACULA 

In scopul obţinerii răspunsului dinamic al rotorilor puşi în discuţie, 
s-au creat în MATLAB următoarele subrutine, utilizate în programul 
DRACULA, în care s-au introdus elementele tronconice de tip CON 12 şi CON 
8: 

Subrutină pentru crearea matricei de masă a elementului conic CON 12 



Subrutină pentru crearea matricei de rigiditate a elementului conic CON 12 



Subrutină pentru crearea matricei giroscopice a elementului conic CON 12 

Subrutina pentru crearea matricei de masă a elementului conic CON 8 
Subrutina pentru crearea matricei de rigiditate a elementului conic CON 8 

Subrutina pentru crearea matricei giroscopice a elementului conic CON 8 

Cu ajutorul acestui program s-au realizat câteva simulări numerice. 
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4.2. Simulări numerice 

Arborele unui rotor se modelează ca grindă deformabilă, iar 
celelalte componente ataşate acestuia se modelează ca fiind corpuri rigide. 

Rotorii din exemplele prezentate ulterior includ porţiuni cilindrice 
şi conice tubulare. După o evaluare critică a mai multor tipuri de elemente finite 
de arbore, pentru analiza în Dinamica rotorilor s-au implementat în programul 
DRACUL A, în MATLAB, cele două tipuri de elemente conice prezentate pe 
larg în paragrafele precedente, elementul conic CON 12 şi elementul conic CON 
8. 

Simularea numerică realizată cu programul DRACUL A a arătat 
diferenţe mari între valorile frecvenţelor proprii obţinute pentru modelele bazate 
pe elemente tip Timoshenko şi pentru modelele bazate pe elemente Bernoulli- 
Euler, ceea ce îndreptăţeşte pe deplin modelarea care ţine cont de inerţia la 
rotaţie şi de deformaţiile de forfecare, adică modelarea cu elemente tip 
Timoshenko. 

Exemplele prezentate în cele ce urmează au fost alese pentru a 
compara performanţele modelării unui rotor tronconic prin elemente finite 
conice cu rezultatele modelării convenţionale prin elemente finite cilindrice în 
trepte ale aceluiaşi rotor. 

4.2.1. Exemplul 1 

Acest rotor, figura 4.1, este un trunchi de con tubular, cu unghiul de 
înclinare de 15 grade, având la capete porţiuni cilindrice tubulare sprijinite pe 
reazeme rigide. 




Fig.4.1 
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In figura 4.2 se prezintă diferite moduri de discretizare a secţiunii 
conice a acestui rotor, utilizând 1, 2, 4, 8, 16 şi 32 de elemente cilindrice sau 
conice. 




0.5 ■ 
0 4 - 
0 3 - 
0.2 ■ 
0.1 ■ 
o - 
-0 1 - 
-0.2 - 
-0 3 - 
-0.4 ■ 
-0 5 - 



D 5 ■ 

LI 4 - 

[I 3 - 

0 2 - 

0 1 - 

0 - 

-0 1 - 

-0 2 - 

-0 3 - 

-0.4 - 

-0.5 ■ 





Fig.4.2 
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Fig.4.2 (continuare) 

Rezultatele obţinute în urma acestor diferite discretizări sunt 
prezentate în tabelele care urmează. 
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In figura 4.3 este prezentată o modalitate generală de discretizare 
utilizând 4 elemente finite cilindrice în trepte. 




Fig.4.3 

Examinând Tabelul 4.1 se pot compara rezultatele obţinute la 
turaţia de 1200 rot/min modelând porţiunea tronconică prin 1, 4 şi 32 de 
elemente de arbore Timoshenko cilindrice, de tip CON 12 şi de tip CON 8. S-au 
comparat primele 8 moduri de precesie. Elementele conice produc frecvenţe 
proprii convergente, monoton descrescătoare odată cu mărirea numărului de 
elemente, pe când elementele cilindrice produc o convergenţă nemonotonă. 



Tabel 4.1. Frecvenţe proprii obţinute utilizând elemente de arbore Timoshenko 





Element cilindric 
Timoshenko 


Element conic CON 12 
de tip Timoshenko 


Element conic CON 8 
de tip Timoshenko 




Frecvenţa proprie [Hz] la 1200 rot/min 


Mod 
Nr. 


4elem. 


32elem. 


lelem. 


4elem. 


32elem. 


lelem. 


4elem. 


32elem. 


1B 


12.400 


12.324 


12.814 


12.350 


12.323 


15.752 


12.581 


12.326 


1F 


17.323 


17.238 


19.325 


17.333 


17.239 


20.297 


17.457 


17.241 


lBbis 


29.891 


29.609 


34.295 


29.916 


29.609 


32.721 


29.827 


29.606 


lFbis 


65.802 


65.463 


67.844 


65.718 


65.461 


69.441 


65.771 


65.460 


2B 


42.677 


41.624 


48.742 


42.742 


41.625 


48.862 


42.875 


41.628 


2F 


44.517 


43.280 


49.832 


44.584 


43.281 


50.293 


44.657 


43.282 


3B 


72.251 


68.274 


82.136 


72.409 


68.276 


77.273 


72.434 


68.277 


3F 


73.731 


69.524 


112.33 


73.937 


69.526 


106.07 


73.928 


69.527 
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în Tabelul 4.2 se prezintă o comparaţie asemănătoare utilizând 
elemente de arbore Bernoulli-Euler, adică elemente pentru care se neglijează 
efectul inerţiei la rotaţie, deformaţia de forfecare, dar se consideră efectul 
giroscopic. In acest caz, convergenţa nemonotonă a frecvenţelor proprii, în cazul 
modelării cu elemente cilindrice, este mult mai pronunţată. 



Tabel 4.2. Frecvenţe proprii obţinute utilizând elemente de arbore Bernoulli-Euler 





Element cilindric 


Element conic CON 12 


Element conic CON 8 




Bernoulli-Euler 


de tip Bernoulli-Euler 


de tip Bernoulli-Euler 




Frecvenţa proprie [Hz] la 1200 rot/min 


Mod 
Nr. 


4elem. 


32elem. 


lelem. 


4elem. 


32elem. 


lelem. 


4elem. 


32elem. 


1B 


17.824 


17.857 


17.888 


17.858 


17.857 


17.888 


17.858 


17.857 


1F 


29.241 


29.300 


29.444 


29.303 


29.302 


29.445 


29.304 


29.303 


2B 


81.426 


81.208 


93.932 


81.248 


81.208 


93.931 


81.247 


81.207 


2F 


119.11 


118.75 


137.75 


118.80 


118.75 


137.75 


118.80 


118.75 


3B 


183.54 


183.80 


247.04 


184.35 


183.81 


247.04 


184.35 


183.81 


3F 


265.58 


265.37 


323.99 


266.05 


265.37 


323.99 


266.05 


265.37 


4B 


323.51 


328.64 


617.54 


332.25 


328.71 


617.51 


332.25 


328.71 


4F 


473.88 


433.33 


914.26 


478.53 


473.36 


914.30 


478.53 


473.36 



în figurile 4.4.a şi b, se prezintă forma câtorva moduri de precesie a 



rotorului discretizat în patru elemente tronconice, pe porţiunea de arbore conic, 
respectiv în 32 de elemente tronconice. 

în figura 4.5 se prezintă orbitele eliptice ale nodului 33, ales 
arbitrar, de pe arbore, pentru cele 8 moduri de precesie analizate. 

Turaţiile critice amortizate sunt date de abscisele punctelor de 
intersecţie dintre curbele frecvenţelor proprii şi linia de sincronism (desenată 
punctat). 

Intersectând liniile din diagramă cu linia verticală de la 1200 
rot/min rezultă frecvenţele proprii din coloana 5, Tabelul 4.1. 

Diagrama Campbell din figura 4.6 prezintă dependenţa de turaţie a 
primelor 14 frecvenţe proprii amortizate, calculate folosind elemente conice 
Timoshenko de tip CON 12 pe porţiunea tronconică. Curbele corespunzătoare 
modurilor „spectrului de ordinul întâi" apar în perechi relativ apropiate, efectul 
giroscopic fiind redus, pentru că rotorul este un arbore fără discuri sau alte 
componente care să introducă acceleraţie Coriolis mare. 
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Mode # 1 B 
17.86 Hz 




Spin speed 1200 rpm 



Mode # 1 F 
29.3 Hz 




Spin speed 1200 rpm 



Mode # 2 B 
81.25 Hz 




Mode # 2 F 
118.8 Hz 




Spin speed 1200 rpm 



Spin speed 1200 rpm 



Mode #3 B 
184.35 Hz 




Mode # 3 F 
266.05 Hz 




Spin speed 1200 rpm 



Spin speed 1200 rpm 



Mode # 4 B 
332.25 Hz 




Mode # 4 F 
478.52 Hz 




Spin speed 1200 rpm 



Spin speed 1200 rpm 



Fig.4.4.a 
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Mode # 1 B 
17.86 Hz 




Spin speed 1200 rpm 



Mode # 1 F 
29.3 Hz 




Spin speed 1200 rpm 



Mode # 2 B 
81.21 Hz 




Mode # 2 F 
118.75 Hz 




Spin speed 1200 rpm 



Spin speed 1200 rpm 



Mode #3 B 
183.81 Hz 




Mode # 3 F 
265.37 Hz 




Spin speed 1200 rpm 



Spin speed 1200 rpm 



Mode # 4 B 
328.71 Hz 




Mode#4F 
473.35 Hz 




Spin speed 1200 rpm 



Spin speed 1200 rpm 



Fig.4.4.b 
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Campbell Diagram 




1 1 ! 1 1 1 

0 1000 2000 3000 4000 5000 6000 

Spin speed, rpm 



Fig.4.6 



Campbell Diagram 




0 1000 2000 3000 4000 5000 6000 



Spin speed, rpm 
Fig.4.7 
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Modurile corespunzătoare „spectrului de frecvenţe de ordinul doi" pot fi 
recunoscute uşor, ele apar în perechi divergente în diagrama Campbell datorată 
efectului giroscopic mare. Aceste moduri apar atunci când rotirea secţiunii 
transversale datorită forfecării este de sens contrar rotirii produse de încovoiere 
(Fig.4.7.1,a). Când rotirile au acelaşi sens ia naştere un mod de precesie obişnuit 
(Fig.4.7.2,b). 

Deformaţia 
totală 

Deformaţia datorită 
încovoierii 




Deformaţia datorită 
forfecării 



Fig.4.7.1,a 




Deformaţia 
totală 



Deformaţia datorită 
încovoierii 

Deformaţia datorită 
forfecării 



Fig.4.7.1,b 



Diagrama Campbell din figura 4.7 prezintă dependenţa de turaţie a 
primelor 12 moduri de precesie, calculate utilizând elemente conice Bernoulli- 
Euler pentru porţiunea tronconică, adică elemente pentru care se neglijează 
efectul inerţiei la rotaţie şi efectul deformaţiei de forfecare. Comparând această 
diagramă cu cea din figura 4.6 rezultă diferenţe relativ mari pentru valorile 
pulsaţiilor proprii. 
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Mode#1 B 
12.324 Hz 




Spin speed 1200 rpm 



Mode#1 F 
17.2433 Hz 




Spin speed 1200 rpm 



Mode#1 B bis 
29.6241 Hz 




Spin speed 1200 rpm 




Spin speed 1200 rpm 



Mode#1 Fbis 
65.4742 Hz 




Spin speed 1200 rpm 




Spin speed 1200 rpm 



Mode #2 B bis 
74.1135 Hz 




Spin speed 1200 rpm 



Mode #3 B 
68.489 Hz 




Spin speed 1200 rpm 



Mode #2 Fbis 
113.1461 Hz 




Spin speed 1200 rpm 




Spin speed 1200 rpm 



Mode #3 B bis 
111.3211 Hz 




Spin speed 1200 rpm 




Spin speed 1200 rpm 



Fig.4.7.2 



Mode #3 Fbis 
149.6711 Hz 




Spin speed 1200 rpm 




Spin speed 1200 rpm 
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In figura 4.7.2 s-au ilustrat formele modurilor proprii de precesie 
ale rotorului discretizat cu elemente de arbore tip Timoshenko, care permit 
vizualizarea spectrului de frecvenţe de ordinul doi corespunzătoare diagramei 
Campbell din figura 4.6. Aceste moduri au fost notate cu „bis". 

In acest prim exemplu a fost ales pentru modelare un caz real de 
rotor conic, având un unghi relativ mic de înclinare. 

Modelarea cu elemente finite conice duce la obţinerea de valori mai 
mari ale pulsaţiilor proprii, în fiecare caz studiat, acestea convergând monoton 
de la valori superioare către valorile „exacte" (considerate a fi foarte apropiate 
de cele obţinute după discretizarea rotorului în 32 de elemente). 

Rezultatele obţinute în urma discretizării rotorului în elemente 
cilindrice în trepte sunt convergente, dar nu numai de la valori superioare, ci 
chiar şi de la valori inferioare către valori apropiate de cele „exacte". 
Convergenţa nemonotonă a acestui tip de modelare este mai evidentă pentru 
elementele Bernoulli-Euler. 

Elementul conic CON 12, cu reducerea Guyan a gradelor de 
libertate pentru deformaţiile de forfecare, conduce la o precizie mai mare a 
valorilor pulsaţiilor proprii ale primelor două moduri de precesie decât 
elementul conic CON 8. Elementul CON 8 conduce la o exactitate mai mare în 
calcul pentru modurile superioare de precesie. 

Modelarea porţiuniii tronconice printr-un singur element conic 
produce erori mari, astfel încât nu este recomandată. 

Dacă inerţia la rotaţie şi forfecarea transversală se neglijează, atunci 
elementele CON 12 şi CON 8 conduc, în esenţă, la aceleaşi rezultate. 

Utilizarea elementelor finite conice pentru porţiunile de arbore 
tronconic conduce la creşterea preciziei calculelor şi la simplificarea procedurii 
de modelare. 
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4.2.2. Exemplul 2 

Rotorul din figura 4.8 este un rotor real produs de firma „Brown 
Boveri", existent în Compartimentul de prezentare al echipamentelor mecanice 
pentru turbinele de abur, având numărul de înregistrare CH-T 1 10 273 D, [214]. 




Fig.4.8 



Rotorul este alcătuit din câteva porţiuni de arbore cilindric în trepte 
şi o porţiune de arbore tronconică. Porţiunile cilindrice au fost modelate 
utilizând elemente de arbore cilindric, iar porţiunea tronconică a fost modelată 
cu elemente conice, discretizând porţiunea tronconică în 4 elemente conice 
(Fig. 4.9), de tip CON 12, apoi de tip CON8 bazate pe elemente Timoshenko. 
Calculele s-au realizat cu programul DRACULA în MATLAB. 



0.5 ■ 

0 4 - 

0 3 - 

0 2 - 




-0.2 - 
-0.3 - 
-0.4 - 
-0.5 - 

I I I I I L 

0 0.2 0.4 0.6 0.8 1 



Fig.4.9 
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In Tabelul 4.3 se compară rezultatele obţinute la turaţia de 6000 
rot/min modelând întregul arbore folosind elemente conice CON 12 şi CON 8 de 
tip Timoshenko, ţinând cont de faptul că elementul cilindric de arbore este un 
caz particular al elementului conic de arbore. 



Tabelul 4.3. Frecvenţe proprii obţinute utilizând elemente de arbore Timoshenko 



Mod 


Element CON 12 


Element CON 8 


Nr. 


de tip Timoshenko 


de tip Timoshenko 




Frecvenţa proprie [Hz] la 6000 rot/min 


1B 


0.288 


0.288 t 


1F 


1.894 


1.899 


1B bis 


2.177 


2.188 


lFbis 


27.045 


27.116 


2B 


20.622 


20.622 


2F 


33.330 


33.375 


2B bis 


26.099 


26.282 


2F bis 


58.671 


58.974 [ 



Se observă că rezultatele obţinute au valori apropiate pentru cele 

două situaţii. 

In figura 4.10 se arată forma primelor 8 moduri de precesie. 

In figura 4.11 se prezintă orbitele de precesie ale nodului 9, ales 
arbitrar, de pe arbore, pentru cele 8 moduri de precesie analizate. 

Diagrama Campbell din figura 4.12 prezintă dependenţa de turaţie a 
primelor 12 frecvenţe proprii amortizate, calculate folosind elemente conice 
Timoshenko de tip CON 12 pe porţiunea tronconică. 

Turaţiile critice amortizate sunt date de abscisele punctelor de 
intersecţie dintre curbele frecvenţelor proprii şi linia de sincronism (desenată 
punctat). 

Intersectând liniile din diagramă cu linia verticală de la 6000 
rot/min rezultă frecvenţele proprii din Tabelul 4.3. 

Diagrama Campbell din figura 4.13 prezintă dependenţa de turaţie 
doar a primelor 8 frecvenţe proprii amortizate. Diagrama a fost arătată pentru o 
claritate mai mare a imaginii. 
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Mode # 1 B 
0.2887 Hz 




Spin speed 6000 rpm 



Mode # 1 F 
1.8941 Hz 




Spin speed 6000 rpm 



Mode # 1 B bis 
2.1774 Hz 




Spin speed 6000 rpm 



Mode # 1 F bis 
27.0455 Hz 




Spin speed 6000 rpm 



Mode # 2 B 
20.6276 Hz 




Spin speed 6000 rpm 



Mode # 2 F 
33.3304 Hz 




Spin speed 6000 rpm 



Mode #2 Bbis 
26.0991 Hz 




Spin speed 6000 rpm 



Mode #2 F bis 
58.6709 Hz 




Spin speed 6000 rpm 



Fig.4.10 
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Campbell Diagram 




0 1000 2000 3000 4000 5000 6000 



Spin speed, rpm 
Fig.4.12 

Campbell Diagram 




0 1000 2000 3000 4000 5000 6000 

Spin speed, rpm 



Fig.4.13 



135 



Dana-Codruţa Vişan Modelarea dinamică a sistemelor rotor-lagăre 



Secţiunea tronconică a fost modelată prin patru elemente conice 
CON 12 şi apoi CON 8, de tip Timoshenko, această variantă de discretizare 
fiind satisfăcătoare din punctul de vedere al preciziei rezultatelor. 

Procedura de modelare a secţiunii tronconice este foarte simplă, şi 
mai mult, rezultate la fel de precise se pot obţine şi în urma modelării acestei 
secţiuni în doar două elemente finite conice. 

Nu este recomandată modelarea secţiunii tronconice printr-un 
singur element conic, deoarece se produc erori mari. 
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4.2.3. Exemplul 3 

In acest exemplu este prezentată o consecinţă a conicizării unui 
arbore. Prin conicizare se înţelege o creştere a unghiului de înclinare suprafeţei 
exterioare a arborelui, în condiţiile păstrării aceluiaşi volum de material. 

Această „conicizare" este utilă la sculele aşchietoare pentru găurire 
şi frezare, pentru că astfel, frecvenţele vibraţiilor de încovoiere şi, deci, 
rigiditatea sculelor aşchietoare pot fi mărite considerabil faţă de cele ale 
arborilor de secţiune constantă având acelaşi volum şi fabricaţi din acelaşi 
material [207, 208]. 

Se consideră un arbore de formă tronconică tubulară, ca în figura 
4.14, din oţel, având volumul de 59 cm . 




Fig.4.14 



Materialul din care este confecţionat rotorul are modulul de 
elasticitate E = 207 GPa, densitatea p = 7700 kg/m , iar turaţia rotorului este de 
400 rad/sec, adică 63,66Hz. 
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Se fac următoarele notaţii: 

bi = raza exterioară mare a trunchiului de con; 
b 2 = raza exterioară mică a trunchiului de con; 
L = lungimea arborelui tronconic; 
t = grosimea peretelui arborelui tronconic. 

Se defineşte „raportul de înclinare" TR, sau „conicitatea arborelui", 
în modul următor: 

TR = 100(b 1 -b 2 ) 
L 

Pentru diverse valori ale raportului de înclinare, respectiv TR = 0, 
TR = 1, TR = 2, TR = 3 şi utilizând modelări cu diferite numere de elemente, 
adică 2 şi 4 elemente, s-au calculat prima pulsaţie proprie a arborelui, cu ajutorul 
programului DRACULA şi s-a reprezentat grafic variaţia acesteia cu conicitatea 
arborelui. 



Pentru ca arborele studiat să-şi păstreze volumul constant de 59 
cm , s-au menţinut constante următoarele mărimi: 

t = 5.4 mm; 
L = 240 mm; 
bi+b 2 = 20 mm. 

Mărimi variabile au fost doar razele exterioare ale arborelui 
tronconic, bi şi b 2 . 

In figura 4.15 s-a arătat valorile primei pulsaţii proprii a rotorului 
pentru diferite conicităţi şi discretizări ale acestuia. 
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b i = 1 Omm 
b 2 =10mm 
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Fig.4.15 



139 



Dana-Codruţa Vişan Modelarea dinamică a sistemelor rotor-lagăre 



Concluzia care se poate trage în urma acestui studiu este că, odată 
cu creşterea conicităţii arborelui, în cazul folosirii aceluiaşi volum de material, 
valoarea primei pulsaţii proprii a acestuia creşte, deci rigiditatea arborelui creşte, 
avantajul imediat fiind acela că o sculă aşchietoare, realizată în condiţiile 
amintite, poate lucra la turaţii mai înalte în condiţii de siguranţă. 

In figura 4.16 s-a reprezentat grafic dependenţa raportului 

pulsaţiilor proprii — cu i=0, 1, 2, 3, de conicitatea arborelui TR. 
Po 




0.995 1 ' ' ' ' ' 1 

0 0.5 1 1.5 2 2.5 3 

Conicitatea TR 



Fig.4.16 

Se observă tendinţa crescătoare a raportului p/po cu mărirea 
conicităţii arborelui, deci creşterea rigidităţii rotorului, respectiv a sculei 
aşchietoare. 
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4.3. Studiul modurilor mixte de precesie 

Exemplul de simulare numerică a fost prezentat pentru a ilustra 
modul în care un cuplaj de tip conservativ produce moduri de precesie mixtă. 
Studiul a fost efectuat cu ajutorul programului de calcul DRACULA, realizat în 
MATLAB, program bazat pe modelarea cu elemente finite a sistemelor rotor- 
lagăre, capabil să sorteze modurile cu precesie directă, inversă sau mixtă. 

In dinamica sistemelor rotor-lagăre, mişcarea de revoluţie a 
rotorilor se numeşte precesie directă (forward), dacă are acelaşi sens cu rotaţia 
rotorului în jurul propriei axe, şi precesie indirectă (backward) - dacă are sens 
contrar rotaţiei. 

Deşi sensul mişcării de precesie este o proprietate locală, ce 
caracterizează un punct dat al rotorului, şi nu una globală, modurile proprii de 
precesie ale unui rotor se numesc, pe scurt, moduri directe (cu precesie directă) 
şi moduri inverse (cu precesie inversă). 

Cu toate acestea este cunoscut faptul că, la o anumită turaţie, unele 
secţiuni ale rotorului pot avea o mişcare de precesie directă, în timp ce alte 
secţiuni au o precesie inversă, fără torsiunea arborelui [16]. Moduri de precesie 
mixtă (directă+inversă) au fost semnalate de Lund [111] şi menţionate recent de 
Jei şi Kim [206] şi Dias şi Allemang [217, 218]. Explicaţia constă în faptul că, la 
multe maşini cu rotor, modurile de precesie sunt într-adevăr moduri pure directe 
sau inverse, sensul precesiei fiind acelaşi în toate punctele în lungul rotorului. 
Aceasta permite o clasificare logică a modurilor proprii, care poate fi extinsă la 
moduri de precesie predominant (nu total) directe sau inverse. 

Sistemele rotor-lagăre nesimetrice şi cu lagăre anizotrope au 
moduri de precesie mixtă produse de cuplaje datorite momentelor giroscopice şi 
rigidităţilor de cuplaj ale lagărelor. Amortizarea din lagăre complică fenomenul, 
modificând alternanţa normală a modurilor inverse cu cele directe. 

In primul studiu detaliat în care se semnalează existenţa modurilor 
de precesie mixtă [111], Lund a analizat un compresor centrifugal cu 8 trepte, cu 
console la capetele rotorului şi cu centrul de greutate situat aproximativ la 
mijlocul distanţei între lagăre. Un astfel de sistem trebuie modelat cu cel puţin 
12 noduri, deci suficient de detaliat pentru a pune în evidenţă modurile mixte. 

La un rotor real, în lagăre cu alunecare, cea mai mică asimetrie 
produce valori diferite ale încărcării şi temperaturii uleiului în lagăre. Deci chiar 
în cazul unor lagăre constructiv identice, coeficienţii de rigiditate şi de 
amortizare sunt diferiţi pentru cele două lagăre ale unui rotor. 

Rezultă că pentru asimetrii uzual întâlnite în practică, curbele 
frecvenţelor proprii amortizate din diagrama Campbell nu se mai intersectează, 
cuplajele intermodale producând forme specifice de diagrame, în care curbele se 
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apropie una de alta şi se îndepărtează fără să se intersecteze („curve veering"). 
In zona de apropiere relativă a frecvenţelor proprii, variaţia bruscă (dar 
continuă) cu turaţia a formei modurilor proprii de precesie produce moduri 
mixte. 

Pentru a ilustra natura modurilor mixte şi evoluţia lor cu variaţia 
turaţiei rotorului, se consideră un rotor nesimetric pe lagăre anizotrope, similar 
cu cel studiat în [209]. Un disc rigid este montat la mijlocul unui arbore de 
secţiune constantă, de lungime 0,44 m, diametrul 9 mm, modul de elasticitate 
longitudinal 2 • 10 11 N/m 2 şi densitate 7800kg/m 3 , figura 4.17. 

Masa discului este 560 kg iar momentele de inerţie masice sunt 18 

2 2 

kgm cel diametral şi 36 kgm cel polar. 

Rotorul este rezemat la capete în lagăre ortotrope, ai căror 
coeficienţi de rigiditate şi amortizare sunt consideraţi independenţi de turaţie. 

S-au calculat modurile proprii de precesie pentru două grupuri de 
valori ale proprietăţilor dinamice ale lagărelor. 




Cazul a : Lagăre identice: 



k yy = 2,2 10 N/m, k„ = 1,1 10 N/m, 



c yy = 2,2 • IO 4 Ns/m, Czz = 1,1 • IO 4 Ns/m. 
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Cazul b : Lagăre puţin diferite: 

k' yy =2,15 - IO 8 N/m, k'^ =1,15 IO 8 N/m, 

c' yy =2,15 IO 4 Ns/m, c' a =1,15 IO 4 Ns/m, 

k^ y = 2,25 • IO 8 N/m, k^ = 1,05 • IO 8 N/m, 

c; y = 2,25 • IO 4 Ns/m , c" a = 1,05 • IO 4 Ns/m . 
Indicii prim şi secund denotă lagărul din stânga, respectiv lagărul din dreapta. 

Variaţia frecvenţelor proprii amortizate şi a rapoartelor de 
amortizare modale pentru primele două moduri, în funcţie de turaţie, este 
prezentată în figura 4.18 (Cazul a) şi figura 4.19 (Cazul b). Datorită ortotropiei, 
fiecare mod de precesie al rotorului constă de fapt din două moduri, unul cu 
precesie inversă, altul cu precesie directă, identificate pe curbe prin literele „B" 
(backward), respectiv „F" (forward). 

Când lagărele au aceleaşi caracteristici dinamice (Fig.4.18), curba 
2.B intersectează dreptele l.F şi l.B. Nu există cuplaj între modul conic 2.B şi 
modurile cilindrice l.F şi l.B. In figura 4.18.b, curbele rapoartelor de amortizare 
nu se intersectează. 

Când lagărele au caracteristici dinamice puţin diferite (Fig.4.19.a), 
curba 2B intersectează dreapta 1 .F la aproximativ 500 rot/min, dar se apropie de 
dreapta l.B la 1900 rot/min, apoi se depărtează. In figura 4.19.b, curbele 
modurilor 2.B şi l.F nu se intersectează, dar au un minim, respectiv un maxim, 
la 500 rot/min, dovedind o interacţiune între moduri, în timp ce curbele 2.B şi 
l.B se intersectează la 1900 rot/min, confirmând cuplajul modal. 

Evoluţia formelor proprii în zonele de turaţii cu interacţiune modală 
este prezentată în figurile 4.20 - 4.24. Deoarece rotorul a fost modelat cu 4 
elemente de arbore de lungimi egale, s-au trasat orbitele de precesie în cele 5 
noduri. De asemenea, s-a trasat deformata în spaţiu a rotorului la momentul t=0 
(linie continuă) şi forma modală la t=7i/2co, după un sfert de perioadă (linie 
întreruptă). Deci mişcarea de precesie în lungul elipselor se face de la punctul 
situat pe linia continuă, la punctul de pe linia întreruptă. în cazul modurilor 
mixte, în dreptul fiecărei elipse s-a notat sensul precesiei F-directă, B-inversă. 

In figura 4.20 se prezintă evoluţia modului l.M între 1700 şi 3000 
rot/min. Modul IM rezultă din cuplarea unui mod conic vertical, cu un mod 
cilindric orizontal care, odată cu creşterea turaţiei rotorului, se transformă în 
mod conic orizontal. 
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La 1700 rot/min, traiectoria discului de la mijlocul rotorului este 
aproape linie dreaptă, marcând trecerea de la precesie directă (în stânga) la 
precesie indirectă (în dreapta). La 1900 rot/min, axa orizontală a elipsei din 
stânga scade, modul cilindric devenind un mod conic cu nodul în afara rotorului. 

La 2000 rot/min, modul orizontal este un mod conic, cu nodul între 
capetele rotorului. Cele două moduri conice, vertical şi orizontal, intersectează 
axa rotorului în puncte diferite, în care orbitele degenerează în linii drepte. 
Modul mixt corespunzător are 3 zone alternative de precesie inversă, directă şi 
din nou inversă. Pe măsură ce turaţia creşte, nodul modului conic orizontal se 
apropie de mijlocul rotorului, modul mixt de precesie devenind predominant 
2.B. 

In figura 4.21 se prezintă evoluţia modului 2.M între 1700 şi 3000. 
La 1700 rot/min modul este încă 2.B. Creşterea turaţiei determină o modificare a 
modului conic orizontal. La 1880 rot/min, capătul din dreapta devine nod, deci 
elipsa este degenerată în linie dreaptă, apoi modul orizontal devine cilindric. 
Deoarece modul vertical rămâne conic, modul de precesie este mixt şi la 300 
rot/min. 

In figura 4.22 se arată evoluţia modului 2.M între 400 şi 600 
rot/min. Dacă în figurile 4.20 şi 4.21 modurile erau aproape plane, în acest caz 
se observă caracterul spaţial al modurilor de precesie ale rotorilor cu amortizare 
în lagăre. 

Deşi curbele frecvenţelor proprii se intersectează şi modul este 
predominant invers (2.B), la 490 şi respectiv 540 rot/min el devine mod mixt. 
Modul conic vertical evoluează prin deplasarea nodului spre capătul din dreapta 
al rotorului, în timp ce deformata la t=0 are o evoluţie spaţială mai complicată. 

In figura 4.23 se arată evoluţia modului 3.M între 480 şi 600 
rot/min, rezultând din combinaţia unui mod cilindric predominant vertical, cu un 
mod conic orizontal. Variaţia fazei modului conic între 490 şi 500 rot/min 
inversează brusc zonele de precesie directă şi inversă în lungul rotorului. 

După cum s-a mai arătat, caracterul mixt al unui mod de precesie 
poate fi omis dacă numărul punctelor în care se calculează sensul mişcării de 
precesie este insuficient. In figura 4.24 se arată al doilea mod propriu al 
sistemului analizat, la 1700 rot/min. Dacă modelul are numai 5 noduri, modul 
apare de precesie inversă. Dacă se consideră încă un nod, între disc şi mijlocul 
jumătăţii din dreapta rotorului, atunci orbita acestuia fiind parcursă în sens 
contrar celorlalte, se pune în evidenţă caracterul mixt al modului respectiv. 
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Concluzia principală a studiului este următoarea: ortotropia 
statorului produce rigidităţi dinamice diferite, deci deformate dinamice diferite, 
ale rotorului, în două direcţii ortogonale, la aceeaşi turaţie. Cuplarea a două 
forme proprii (chiar puţin) diferite ale rotorului din două plane ortogonale 
produce moduri de precesie mixtă. 

Acest lucru se poate explica mai uşor pentru sisteme giroscopice 
conservative, care au moduri de precesie plane. în acest caz, deformatele în două 
plane ortogonale corespund formelor modale la două momente având o diferenţă 
de timp de un sfert de perioadă de precesie. Dacă cel puţin una din cele două 
curbe intersectează axa rotorului pe lungimea acestuia, atunci, în punctul 
respectiv, orbita de precesie degenerează intr-o linie dreaptă, care separă zonele 
de precesie directă şi inversă în lungul rotorului. 

Probabil că modurile de precesie mixtă au fost neglijate în trecut 
deoarece sistemele rotor-lagăre considerate în exemplele teoretice erau fie 
simetrice, fie cu lagăre izotrope, iar mişcarea era analizată într-un număr relativ 
redus de puncte în lungul rotorului. în multe cazuri, modurile mixte sunt 
predominant directe sau inverse, cu zone reduse de precesie în sens contrar, 
astfel că, dacă numărul nodurilor într-un model cu elemente finite este redus, 
atunci nu se observă caracterul mixt al modului de precesie respectiv. 
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5. Rezultate experimentale 



In acest capitol, sunt prezentate rezultatele obţinute de autoare în 
urma unor încercări în laboratorul Catedrei de Rezistenţa materialelor din 
Universitatea "Politehnica" din Bucureşti, încercări efectuate pe un Rotor Kit, 
producţie a firmei Bently Nevada. 

In prima parte a capitolului este descris standul de încercări, cu 
posibilităţile şi caracteristicile funcţionale şi sunt puse în evidenţă 
particularităţile importante ale sistemului rotor-lagăre. 

In cea de a doua parte a acestui capitol sunt prezentate rezultatele 
experimentale obţinute în urma măsurărilor din laborator şi o comparaţie a 
acestora cu rezultatele obţinute, pentru acelaşi rotor, utilizând programul 
DRACUL A în MATLAB. 



5.1. Standul rotoric Bently Nevada 

Standul folosit este un ROTOR KIT BENTLY NEVADA (model 
nr. 24755). In continuare sunt prezentate caracteristicile generale ale standului şi 
o scurtă descriere a acestuia. 



Caracteristici generale ale standului: 



Caracteristici mecanice 

-Turaţia motorului 
-Dimensiuni 



-până la 
-lungime 
-lăţime 
-înălţime 



12 000rot/min; 
787 mm; 
152 mm; 
168 mm. 



Caracteristici electrice 

-Motorul electric 
-Traductorul fără contact 



-220 V, 50 Hz; 
-la ieşire: 18 V± 1 V 
la max. 100 mA. 
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Greutăţi 

-Rotor (fără fundaţie) 1 3,80 kg; 

-Discul mare 0,81 kg; 

-Discul mic 0,605 kg; 

-Lagărul de alunecare plus 
sistemul de ungere 1,87 kg; 

-Sistemul pentru modelarea 
fenomenelor aerodinamice 0,74 kg; 

-Placa cu amplificatoarele 
traductorilor 2,21 kg. 

Accesorii 

-Tahometru N 2603; 
-Osciloscop TEKTRONIX 2205. 



Descrierea standului: 

Standul este prezentat schematizat în figura 5.1. Se disting 
următoarele părţi mai importante: 

1. motorul de acţionare al arborelui; 

2. traductorul de control al turaţiei arborelui; 

3. traductorul pentru semnalul de fază; 

4. şurubul de simulare a frecărilor între rotor şi stator; 

5. discul mare şi cel mic; 

6. traductorii de deplasare, fără contact, montaţi pe direcţiile Y-Z; 

7. cutia de comandă a turaţiei motorului; 

8. cuplajul elastic motor-arbore; 

9. lagărul 

10. placa de bază. 

Standul a fost prins prin placa sa de bază cu şuruburi de o fundaţie 
realizată prin turnarea de ciment într-o cutie de lemn. In scopul manevrării 
uşoare a conexiunilor electrice, placa pe care se află montate amplificatoarele, 
care amplifică semnalele traductorilor, a fost fixată de o laterală a cutiei de 
lemn, figura 5.2. 

Turaţia motorului electric de acţionare poate varia continuu între 
limitele 1200 + 12000 rot/min, existând posibilităţi de modificare a ratei de 
creştere sau de scădere a turaţiei într-o plajă de turaţii prestabilită, sau de 
menţinere a ei la o valoare constantă. 
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Lagărele, notate cu 9 în figura 5.1, deşi aparent rigide, 
caracteristicile lor dinamice se apropie mai degrabă de cele ale lagărelor cu 
alunecare decât ale celor cu rulmenţi. Astfel, aşa cum se observă din figura 5.3, 
lagărul îi permite totuşi arborelui o rotire relativ mică - bucşa de bronz grafitat, 
fiind prinsă elastic de corpul lagărului, are posibilitatea de a se roti în jurul unei 
axe perpendiculare pe axa lagărului, atât cât îi permite jocul dat de cotele de 
execuţie ale celor două piese. 



corpul lagărului M de 



arborele 

\ 
\ 




bucşa c|e bronz 
/ grafitat 




Fig.5.3 

In plus, inelul de cauciuc duce la modificări importante ale 
caracteristicilor dinamice ale lagărului cu turaţia, atât ale rigidităţii, cât şi ale 
amortizării. Variaţia acestora cu turaţia arborelui a fost determinată experimental 
de către producător pentru acest tip de lagăr şi este reprodusă în figura 5.4, cu 
observaţia că valorile coeficienţilor dinamici se modifică odată cu accentuarea 
uzurii bucşei de bronz grafitat. 
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5.2. Modul de lucru 

Măsurătorile efectuate au avut ca scop determinarea practică a 
turaţiei critice a unui rotor din oţel având secţiuni conice şi compararea acesteia 
cu turaţia critică obţinută în urma simulării numerice cu programul DRACULA 
înMATLAB. 

Schema lanţului de măsură utilizat este prezentată în figura 5.5. 




Fig.5.5 

In figura 5.5 s-au reprezentat: 1 -Generator de frecvenţă GF22; 2- 
Amplificator LV103; 3-Vibrator electromagnetic VEB RFT; 4-Tijă; 5-Captor de 
forţă 8200 B&K; 6-Inel de bronz; 7-Amplificator de sarcină 2626 B&K; 8- 
Voltmetru ORV2-RFT; 9-Captor de acceleraţii 4370 B&K; 10-Vibrometru 2511 
B&K; 1 l-Filtru trece bandă 1621 B&K; 12-Arbore; 13-Osciloscop Tektronix 
2205 ; 1 4-Tahometru N 2603 . 
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Aparatura utilizată se poate observa şi în vederea de ansamblu din figura 5.6. 




Fig.5.6 



S-a procedat în modul următor: S-a măsurat deplasarea absolută pe 
lagărul rotorului cu ajutorul unui accelerometru piezoelectric tip 4370 B&K şi 
un vibrometru 251 1 B&K legat la un filtru trece bandă acordabil tip 1621 B&K 
(Fig.5.7). 




Fig.5.7 
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în figura 5.7 se poate observa şi cuplajul elastic ce există între 
motor şi arbore. 

Turaţia a fost măsurată cu tahometrul tip N 2603 şi captorii fără 
contact de tip Bently Nevada cu care este echipat standul. S-a putut vizualiza 
orbita mişcării de precesie a centrului arborelui într-o secţiune din vecinătatea 
lagărului, pe ecranul unui osciloscop Tektronix 2205 (Fig.5.8). 




Fig.5.8 



Orbita eliptică a mişcării de precesie a arborelui conic se observă 
mai clar în figura 5.9. 




Fig.5.9 
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Semnalul electric preluat de la vibrometru a fost înregistrat pe hard 
disk-ul unui laptop (Fig.5.10), prevăzut cu placă de achiziţie de date. 

— ^ — 

iî 



Fig.5.10 

Semnalul înregistrat a fost apoi prelucrat folosind Toolbox-ul 
Signal Processing din MATLAB. 

Arborele conic utilizat la efectuarea măsurătorilor, figurile 5.1 l.a,b, 
a fost confecţionat din oţel, în atelierul mecanic al Catedrei de Rezistenţa 
materialelor din Universitatea Politehnica Bucureşti. 




Fig.5.11.a 
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Fig.5.11.b 

Acest arbore are lungimea totală de 450 mm, diametrul maxim de 
21 mm, diametrul minim de 9 mm, iar porţiunile conice au o înclinare de 2°, 
figura 5.12. 




0 4 0.6 
Lungimea 



Fig.5.12 



Utilizând programul DRACULA în MATLAB şi discretizând 
fiecare porţiune tronconică în câte 4 elemente finite conice, prima turaţie critică 
calculată a rezultat de 7725 rot/min, adică 128,75 Hz. 
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5.3. Interpretarea rezultatelor obţinute 



Arborele cu porţiuni tronconice a fost montat pe standul Rotor Kit 
în lagăre de bronz grafitat (Fig.5.13). Folosind facilităţile oferite de către 
sistemul de comandă al motorului de antrenare al rotorului, turaţia acestuia a 
putut fi reglată şi menţinută constantă la valori discrete cuprinse în plaja 1200- 
12000 rot/min, cu un pas de 600 rot/min. Corespunzător fiecărei turaţii, s-au 
făcut înregistrări pe durata a trei secunde, cu o rată de eşantionare de o miime de 
secundă. Ulterior, datele achiziţionate au fost prelucrate utilizând tool-box-urile 
specifice (Signal Processing) din MATLAB, obţinându-se curbele de răspuns în 
frecvenţă. 




Fig.5.13 



Este cunoscut faptul că un proces aleator poate fi descris fie în 
domeniul timp, prin valorile unei funcţii de timp h(t), fie în domeniul frecvenţă, 
prin funcţia H(f), ca funcţie de frecvenţa f , în care -oo </ < +qo. Transformata 
Fourier a uneia dintre funcţii o dă pe cea de-a doua, ele putând fi scrise sub 
forma: 



H(f) = £h(t)e 2 " ift dt < 51 > 
h(t) = fH(f>- 2 *df < 52 > 
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In cazul de faţă, h(t) reprezintă variaţia deplasării în timp, măsurată 
pe lagărul rotorului şi înregistrată pe hard-disk-ul laptop-ului. O secvenţă a unei 
astfel de înregistrări este prezentată în figura 5.14 şi a fost făcută la turaţia de 
1800 rot/min. Privind graficul din figura 5.14, se observă că în semnalul 
înregistrat predomină componenta primei armonici (cea corespunzătoare 
frecvenţei de rotaţie a rotorului), peste care se suprapune o componentă de înaltă 
frecvenţă, aceasta ieşind însă din domeniul frecvenţelor de interes şi 
corespunzând, probabil, semnalelor induse de periile colectorului motorului 
electric. 



2 




TIMP [miimi de secunda] 



Fig.5.14 

Utilizând facilităţile oferite de către programele MATLAB, 
fişierele cu datele înregistrate au fost prelucrate numeric cu ajutorul funcţiilor de 
forma FFT - conform relaţiei (5.2) - obţinându-se astfel răspunsul în frecvenţă, 
răspuns reprezentat grafic în figura 5.15. Aşa cum se observă din această figură, 
spectrul de răspuns în frecvenţă are o singură componentă, corespunzătoare 
frecvenţei de 30 Hz, cauzată de dezechilibrul masic al rotorului (rotorul 
tronconic a fost realizat în atelierul mecanic al catedrei de Rezistenţa 
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materialelor, după care a fost montat direct pe stand, fără a i se face o 
echilibrare). 




40 60 80 

FRECVENTA [Hz] 



Fig.5.15 

Valoarea amplitudinii nu a prezentat în cazul de faţă vreun interes, 
urmărindu-se doar determinarea turaţiei critice a rotorului şi compararea ei cu 
cea obţinută teoretic, prin simulare numerică. Ca urmare, nu s-a făcut o 
etalonare a aparatelor de măsură în scopul obţinerii unei valori exacte a 
amplitudinii răspunsului la dezechilibru masic, interesând doar modificări 
cantitative ale acesteia, respectiv la trecerea de la o turaţie la alta a rotorului. 
Obţinerea unei valori maxime a amplitudinii măsurate va indica existenţa unei 
turaţii critice, aceasta fiind o metodă relativ simplă de determinare a turaţiilor 
critice. 

In cazul de faţă, prin suprapunerea curbelor de răspuns în frecvenţă, 
curbe obţinute prin prelucrarea datelor experimentale înregistrate, se poate 
vizualiza cu uşurinţă valoarea frecvenţei în dreptul căreia apare rezonanţa, deci 
valoarea acelei frecvenţe la care viteza de rotaţie a arborelui devine egală cu 
pulsaţia proprie a sistemului elastic rotor-lagăre (Fig.5.16). Astfel, folosind 
funcţii de interpolare (MATLAB) corespunzătoare coordonatelor punctelor de 
maxim, puncte aflate pe curbele de răspuns în frecvenţă specifice fiecărei turaţii 
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la care s-au făcut înregistrări, s-a obţinut o singură valoare critică în domeniul 
frecvenţelor studiate, f cr =126 Hz, care corespunde turaţiei critice n cr =7560 
rot/min (Fig.5.17). 
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Fig.5.16 
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Fig.5.17 

Comparând valoarea turaţiei critice obţinută experimental (7560 
rot/min) cu cea obţinută prin simulare numerică şi folosirea elementului conic în 
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metoda elementelor finite 7725 rot/min, a rezultat o diferenţă An cr = 165 rot/min, 
care corespunde unei erori de 2,18% raportată la valoarea măsurată 
experimental. 

O explicaţie, a faptului că turaţia critică obţinută prin simulare 
numerică are o valoare mai mare decât cea obţinută experimental, poate fi dată 
dacă se ia în considerare rigiditatea relativ mai mare a modelului cu elemente 
finite al sistemului, faţă de rigiditatea existentă în cazul real. Rezultate numerice 
apropiate de cele experimentale pot fi obţinute printr-o discretizare mai fină a 
rotorului, deci prin mărirea numărului de elemente finite, ceea ce va duce la 
mărirea elasticităţii sistemului rotor-lagăre. Este de altfel cunoscut faptul că cu 
cât se doresc rezultate ale modurilor proprii superioare, obţinute prin aplicarea 
metodei elementelor finite, mai apropiate de cele reale, cu atât mai fină trebuie 
efectuată discretizarea sistemului. 

De multe ori, reprezentările grafice ale curbelor de răspuns în 
frecvenţă se fac folosind scări logaritmice în scopul unei vizualizări optime a 
răspunsului şi în domeniul frecvenţelor joase şi înalte. Curbele de răspuns în 
frecvenţă din figurile 5.16 şi 5.17 au fost trasate folosind şi astfel de scări şi sunt 
prezentate în figurile 5.18 şi respectiv 5.19. Rezultă clar şi din aceste grafice că 
în plaja de turaţii 1800-12000 rot/min apare numai o singură turaţie critică, în 
dreptul valorii de 7560 rot/min. Cea de a doua turaţie critică nu a mai putut fi 
obţinută experimental deoarece valoarea ei depăşea cu mult turaţia maximă, 
capabilă, a motorului de antrenare a rotorului (~ 12230 rot/min). 




10° 10 1 10 2 10 3 
FRECVENTA [Hz] 



Fig.5.18 
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Fig.5.19 

Faptul că apar intersecţii ale curbelor de răspuns în frecvenţă, în 
domeniul frecvenţelor relativ înalte (în zonele descendente - cu pantă negativă - 
ale acestor curbe), poate fi explicat prin neliniarităţile existente în sistem şi 
cauzate în principal de prezenţa inelului de cauciuc aflat între bucşa de cupru 
grafitat al lagărului şi corpul acestuia ( a se vedea figura 5.3). 

Un alt mod de vizualizare al rezultatelor este cel în cascadă 
(waterfall) şi este prezentat în figura 5.20, în care curbele de răspuns în 
frecvenţă, corespunzătoare turaţiilor rotorului la care s-au făcut înregistrările, 
sunt reprezentate succesiv, una în spatele celeilalte. Pe abscisă sunt trecute 
frecvenţele care caracterizează răspunsul sistemului, iar pe ordonată turaţiile la 
care s-a făcut înregistrarea datelor experimentale. Dacă pe ambele axe s-ar fi 
ales aceeaşi mărime fizică (fie numai frecvenţa exprimată în Hz, fie numai 
turaţia - în rot/min), maximele de pe curbele de răspuns s-ar fi aliniat în lungul 
primei bisectoare din plan şi ar fi demonstrat astfel clar că ele sunt cauzate de 
existenţa dezechilibrului masic al rotorului. 

în concluzie, efectuarea de măsurări numai la o singură turaţie, sau 
la un număr relativ mic de turaţii ale rotorului, nu oferă informaţii referitoare la 
turaţiile critice ale sistemului. Pentru a obţine turaţiile critice sunt necesare un 
număr relativ mare de înregistrări, la diverse turaţii ale rotorului pentru a putea fi 
posibilă găsirea valorii maxime extreme a răspunsului, aceasta corespunzând 
turaţiei critice. 
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Fig.5.20 

Eroarea relativ mică între valorile obţinute prin simulare numerică 
şi cele experimentale oferă siguranţă în folosirea elementului conic implementat 
în programul de elemente finite, certificând astfel robusteţea acestuia. 
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6. Concluzii 



6.1. Modelarea dinamică a sistemelor rotor-lagăre 

Pentru modelarea dinamică a sistemelor rotor-lagăre, în lucrarea de 
faţă, din cele două metode posibile: metoda matricelor de transfer şi metoda 
elementelor finite, a fost abordată cea de a doua, datorită posibilităţilor mai largi 
de modelare a tuturor părţilor componente a maşinilor cu rotor şi a diverselor 
efecte ce apar în funcţionarea acestora. 

Obiectul tezei, însă, 1-a constituit numai modelarea sistemelor 
rotor-lagăre, adică modelarea arborelui, a discurilor şi a lagărelor, nu şi celelalte 
efecte privind interacţiunea rotorului cu etanşările, piedestalurile şi fundaţia 
maşinii. 

Din marea diversitate de elemente finite existente cu aplicabilitate 
în dinamica rotorilor, pentru modelarea arborilor s-a optat pentru elementul finit 
de tip bară cu două noduri la capete. Au fost stabilite matricea de rigiditate, 
matricea de inerţie coerentă, matricea giroscopică coerentă şi vectorul forţelor 
de dezechilibru 

Discurile rigide s-au modelat declarând existenţa şi caracteristicile 
acestora în nodurile elementelor. Au fost stabilite matricea de inerţie, matricea 
giroscopică şi vectorul forţelor de dezechilibru pentru discurile rigide axial 
simetrice. 

Lagărele au fost modelate, de asemenea, declarând în nodurile 
elementelor existenţa şi caracteristicile lor mecanice şi definind matricele celor 
opt coeficienţi dinamici elastici şi de amortizare. 

Problema modelării arborelui a fost mai pe larg dezbătută, fiind de 
altfel şi cea mai importantă. 

Ca urmare, la programul de Dinamica rotorilor în uz la Catedra de 
Rezistenţa materialelor din Universitatea Politehnica Bucureşti, care modela 
arborele prin elemente finite de bară de tip Bernoulli-Euler, s-a adăugat 
posibilitatea de modelare a arborilor cu elemente finite de bară de tip 
Timoshenko ce ţin cont atât de încovoiere, ca şi elementul de bară de tip 
Bernoulli-Euler, cât şi de deformaţiile datorită forfecării şi de inerţia la rotaţie. 

S-a implementat, în mediul MATLAB, un element de bară 
tronconic de tip Timoshenko, pentru a înlocui procedeul de modelare al arborilor 
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având porţiuni tronconice prin elemente cilindrice în trepte, cu modelarea 
directă prin elemente de bară tronconice. 

In literatura de specialitate există o multitudine de elemente finite 
pentru modelarea arborilor conici, dar experienţa căpătată în cursul studiilor 
efectuate, demonstrează că cele mai performante sunt: 1) elementul conic 
descris în lucrarea [58], preluat de autoare şi numit CON 8, care este un element 
finit de bară axial simetric (având secţiunea transversală circulară sau inelară), 
de tip Timoshenko, cu două noduri la capete, cu câte patru grade de libertate pe 
nod, adică deplasarea laterală şi rotirea secţiunii transversale, în câte două plane 
perpendiculare între ele şi 2) elementul conic descris în lucrarea [70] şi numit de 
autoare CON 12. Şi acest element se bazează pe elementul de bară de tip 
Timoshenko având două noduri la capete, secţiunea transversală circulară sau 
inelară, dar el are câte şase grade de libertate pe nod: deplasarea transversală, 
rotirea axei neutre şi unghiul de forfecare, în câte două plane perpendiculare 
între ele. 

Pentru fiecare din aceste două elemente finite conice s-au scris, în 
mediul MATLAB, subrutinele pentru calculul matricei de masă coerente, 
matricei de rigiditate şi matricei giroscopice coerente în vederea completării 
programului existent. 

Matricele elementului CON 8 au dimensiuni 8x8. Matricele 
elementului CON 12 au dimensiuni 12x12, ceea ce necesită o procedură de 
condensare statică pentru micşorarea ordinului matricelor la dimensiuni 8x8, cu 
scopul de a fi compatibile cu restul programului. 

Cu programul astfel realizat s-au făcut o serie de simulări numerice 
din care au rezultat câteva particularităţi: 

1 . Existenţa aşa numitului „spectru de frecvenţe de ordinul doi" al rotorilor 
modelaţi cu elemente de bară de tip Timoshenko. Spectrul de frecvenţe de 
ordinul doi apare atunci când rotirea secţiunii transversale datorită 
forfecării este de sens contrar rotirii produse de încovoiere. Când aceste 
rotiri au acelaşi sens există un mod de precesie clasic. 

2. Efectul de rigidizare introdus de conicizarea unui arbore. S-a studiat 
efectul variaţiei conicităţii unui arbore tubular de grosime constantă, 
realizat dinvolum de material constant, asupra rigidităţii acestuia. S-a 
observat că, rigiditatea arborelui creşte odată cu creşterea conicităţii, la un 
volum constant de material. 

3. Existenţa unor moduri de precesie mixtă, în care unele puncte ale 
rotorului au mişcare de precesie directă, iar alte puncte ale aceluiaşi rotor 
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au precesie inversă, fără a apărea torsiunea arborelui. Prin mod mixt de 
precesie se înţelege un mod compus rezultat din combinarea a două 
moduri diferita care au aceeaşi frecvenţă. Se arată că sensul mişcării de 
precesie nu este o proprietate globală ci una locală, deşi există moduri de 
precesie exclusiv directă sau inversă. Aceste efecte pot trece neobservate 
dacă numărul punctelor în care se calculează sensul mişcării de precesie 
este insuficient. Modurile mixte pot apare la sistemele cu lagăre 
anizotrope şi nesimetrice în lungul arborelui. Chiar când rotorul este axial 
simetric, ortotropia statorului face ca prin cuplarea unor moduri cu forme 
proprii diferite să rezulte moduri cu precesie mixtă. In cadrul studiului 
efectuat s-au rezolvat dificultăţile ce apar la trasarea automată a diagramei 
Campbell, utilizată la determinarea turaţiilor critice, reliefând faptul că 
pentru anumite tipuri de rotori apare alura convergent-divergentă a 
ramurilor diagramei (veering), despre care s-a vorbit în paragraful 4.3. 
Utilizând cele două elemente finite conice CON 8 şi CON 12, s-a realizat 
un studiu comparativ pornind de la modelarea unei porţiuni tronconice 
prin elemente conice sau cilindrice în trepte, comparând apoi utilizarea în 
calcule a elementului de bară tip Bernoulli-Euler sau tip Timoshenko. S-a 
ajuns la următoarele concluzii: a) Modelarea porţiunii tronconice a unui 
rotor cu elemente finite conice duce la obţinerea de valori mai mari ale 
pulsaţiilor proprii, în fiecare caz studiat, acestea convergând monoton de 
la valori superioare către valorile „exacte" (considerate a fi foarte 
apropiate de cele obţinute după discretizarea rotorului în 32 de elemente); 
b) Rezultatele obţinute în urma discretizării rotorului în elemente 
cilindrice în trepte sunt convergente, dar nu numai dinspre valori 
superioare, ci chiar şi dinspre valori inferioare către valori apropiate de 
cele „exacte". Convergenţa nemonotonă a acestui tip de modelare este 
mai evidentă pentru elementele Bernoulli-Euler; c) Elementul conic CON 
12, cu reducerea Guyan a gradelor de libertate pentru deformaţiile de 
forfecare, conduce la o precizie mai mare a valorilor primelor două 
moduri de precesie decât elementul conic CON 8. Elementul CON 8 
conduce la o precizie mai mare în calcul pentru valorile pulsaţiilor proprii 
corespunzătoare modurilor superioare de precesie; d) Modelarea secţiunii 
tronconice printr-un singur element conic produce erori mari, astfel încât 
nu este recomandată; e) Dacă inerţia la rotaţie şi forfecarea transversală se 
neglijează, atunci elementele CON 12 şi CON 8 conduc, în esenţă, la 
aceleaşi rezultate; f) Utilizarea elementelor finite conice pentru porţiunile 
de arbore tronconic conduce la creşterea preciziei calculelor şi la 
simplificarea procedurii de modelare. S-a măsurat experimental prima 
turaţie critică a unui rotor cu porţiuni tronconice, după efectuarea 
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calculului acesteia cu programul DRACULA în MATLAB în vederea 
validării programului. Totuşi turaţia critică obţinută prin simulare 
numerică, la o discretizare a unei porţiuni tronconice în patru elemente 
conice, a avut o valoare mai mare decât cea obţinută experimental, 
eroarea fiind de 2.18%, explicaţia fiind rigiditatea relativ mai mare a 
modelului cu elemente finite, faţă de rigiditatea existentă în cazul real, dar 
rezultatele numerice sunt cu atât mai apropiate de cele experimentale, cu 
cât discretizarea este mai fină. 

6.2. Contribuţii originale 

Lucrarea de faţă a apărut ca o necesitate a rezolvării unor probleme 
de modelare dinamică a sistemelor rotor-lagăre în cadrul preocupărilor Catedrei 
de Rezistenţa materialelor din Universitatea Politehnica Bucureşti. 

Contribuţiile principale ale autoarei sunt următoarele: 

1. Prezintă o sinteză a modelării barelor şi sistemelor de bare cu elemente 
finite de tip Timoshenko în analiza dinamică. 

2. Prezintă o clasificare şi o caracterizare detaliată a elementelor finite 
tronconice. 

3. Realizează şase subrutine de calcul în mediul MATLAB pentru modelarea 
elementelor conice de tip Timoshenko CON 8 şi CON 12, primul fiind un 
element de bară conic având două noduri la capete şi patru grade de 
libertate pe nod, cu matrice de ordin 8x8, iar al doilea un element cu două 
noduri la capete şi şase grade de libertate pe nod, care necesită şi o 
condensare statică pentru micşorarea ordinului matricelor de la 12x12 la 
8x8. 

4. Pune în evidenţă spectrul de frecvenţe de ordinul doi al arborilor modelaţi 
cu elemente de bară de tip Timoshenko. 

5. Studiază modurile mixte de precesie, în care unele puncte ale unui arbore 
au precesie directă, iar altele, simultan, au precesie inversă şi apoi face 
analiza diagramelor Campbell trasate pentru rotori cu precesie mixtă şi 
modelaţi prin elemente de tip Timoshenko. 
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6. Realizează un studiu comparativ între valorile măsurate experimental şi 
apoi calculate numeric, pentru o turaţie critică a unui arbore cu porţiuni 
tronconice. 

7. Completează programul DRACULA (Dinamica Rotarilor Asistată de 
CalcULAtor), sub forma unui Toolbox în mediu MATLAB, pentru a 
permite calculul răspunsului dinamic al retorilor modelând arborii cu 
elemente de bară Timoshenko tronconice sau chiar cilindrice, acestea din 
urmă fiind de fapt un caz particular al celor conice. 

8. Prezintă o amplă bibliografie (peste 200 de titluri) referitoare la domeniul 
studiat. 

Pe parcursul dezvoltării tezei, autoarea a participat la contracte de 
cercetare [211, 212, 213] cu tematică legată de preocupările din această lucrare, 
în care s-au publicat o serie de rezultate obţinute folosind programul de calcul D 
elaborat de autoare. 

Rezultate parţiale legate de tema abordată în teza de doctorat, au 
fost publicate în [179, 180, 181, 182, 183, 184, 186,210]. 

Studiile şi cercetările din ultimii ani au permis autoarei să-şi 
formeze o imagine de ansamblu asupra problematicii modelării dinamice a 
sistemelor rotor-lagăre şi au deschis perspectiva unor cercetări viitoare: 

1. Modelarea arborilor prin elemente finite diferite de cele de tip bară 
(modele de tip C°). 

2. Continuarea studiului modurilor mixte de precesie, având în vedere 
interesul manifestat pe plan internaţional faţă de această problemă [217, 
218]. 

3. Studiul influenţei flexibilităţii discului asupra comportării dinamice a 
maşinilor cu rotor. 

4. Studiul retorilor din materiale compozite. 
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Calculul funcţiilor de formă pentru elementul de arbore Bernoulli-Euler 



Calculul funcţiilor de formă 
în planul x-y 

Funcţiile de formă au următoarele proprietăţi importante: 

1 . Suma funcţiilor de formă are valoarea egală cu 1 în orice nod; 

2. Corespunzător coordonatei n, funcţia de formă N n =l, iar celelalte 
funcţii de formă au valori nule. 




Determinarea funcţiilor de formă N^) şi N x (^) 
în planul x-y 



relaţii: 



Pentru determinarea funcţiilor de formă se folosesc următoarele 



v(0 = A 3 ^ + A 2 4 2 + A^ + A 0 



De asemenea 



vfe)=N 1 fe)v i +N 2 feVi +N 3 fe)vj +N 4 (0i|/j 
[vfeMifefo +N 2 fe) Vi +N 3 (0vj +N 4 fe>|/j 
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Atunci 



v(0) = l ^A 0 =l 

\|/(<)) = () => ^ A, =0 

v(l) = 0 =>A 3 +A 2 +A 1 +A 0 =0 

v (i)=o =>^a 3+ ^a 2 + 1a 1 =o 



De unde A 0 = 1, A , = 0, 
apoi se rezolvă sistemul: 



A 2 + A 3 = -l 
2A 2 + 3A 3 = 0 



Ca urmare 



iar soluţiile finale sunt: 
A 0 = 1 

A, = 0 

A 2 = -3 
A 3 = 2 



N 1 (0=A3^ 3 +A 2 ^ 2 + A 1 ^ + A 0 



Deci funcţia de formă pentru translaţii este: 



Pe de altă parte 



Adică 



Deci funcţia de formă pentru rotaţii este: 
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Determinarea funcţiilor de formă N 2 (^) şi N 2 (^) 

în planul x-y 



v(0) = 0 
y(0)=l 
v(l) = 0 

v(i)=o 



A 0 =0 

A 3 + A 2 +A l +A 0 =0 

3 . 2.1. . 
— A, + — A 2 + — A, =0 
L L L 



De unde A 0 = 0, A , = -L, 

apoi se rezolvă sistemul: 



A 2 + A 3 
2A 2 + 3A 3 



Ca urmare 



Deci 



L 
L 



iar soluţiile finale sunt: 



A 0 = 


0 


A 1 = 


L 


A 2 = 


-2L 


A 3 = 


L 



N 2 fe) = A3^+A 2 ^ 2 + A^ + A 0 



Pe de altă parte 



N 2 fe) = Lfc-2^ 2 ^J 



Adică 



Deci 



N 2 fe)=^L^ 2 -^2L^ + iL 



N 2 fe) = 3^ 2 -4^ + l 
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Determinarea funcţiilor de formă N 3 (^) şi N 3 (^) 

în planul x-y 



v(0) = 0 
y(0) = 0 
v(l) = l 

v(i)=o 



A 0 =0 

^A, =0 
L 1 

A 3 + A 2 +A 1 + A 0 =1 

3 . 2.1. . 
— A 3 + — A 2 +— A, =0 
L L L 



De unde A 0 = 0, A , = 0, 
apoi se rezolvă sistemul: 



A 2 + A 3 



= 1 



2A 2 + 3A 3 = 0 



Ca urmare 



iar soluţiile finale sunt: 
A 0 = 0 

A, = 0 

A 2 = 3" 
A 3 = -2 



N 3 0;) = A3^+A^ 2 +A^ + A ( 



Deci 



N 3 fe) = 3^ 2 -2^ 



Pe de altă parte 



Adică 



N 3 ($)=-^ 2 +^3$ 



Deci 
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Determinarea funcţiilor de formă N 4 (^) şi N 4 (%) 

în planul x-y 

v(0) = 0 
¥ (0) = 0 
v(l) = 0 

v(i)=i 

De unde A 0 - 0, A , - 0, 
apoi se rezolvă sistemul: 

A 2 + A 3 = 0 
2A 2 + 3A 3 = L' 

Ca urmare 

N 4 (0=A 3 ^ + A 2 *; 2 +A^ + A 0 

Deci 

N 4 fe)=L(-^ 2 +^ 3 ) 

Pe de altă parte 
Adică 
Deci 

N 4 ©=3^ 2 -2^ 



A 0 =0 

^A, =0 
L 1 

A 3 +A 2 + A 1 + A 0 =0 

3 . 2 1 

— A 3 + — A 2 +— A, =1 

L 3 L 2 L 1 



iar soluţiile finale sunt: 



A 0 = 


0 


A 1 = 


0 


'a 2 = 


L 


A 3 = 


-L 
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Calculul funcţiilor de formă 
în planul x-z 

Folosind proprietăţile funcţiilor de formă, se calculează expresiile 
acestora şi în planul x-z. 



Determinarea funcţiilor de formă N x (^) şi Nx(^) 

în planul x-z 

Pentru determinarea funcţiilor de formă se folosesc următoarele 

relaţii: 

w(4)=A 3 4 3 +A 2 4 2 +A 1 4 + A 0 
'-<pfe)=f A 3 $ 2 +£a 2 $+±A 1 

De asemenea 

' wfe)=N 1 fe)w i -Nafefo +N 3 (0wj -N^fo 
><pfe) = NifeK -N 2 fe>Pi +N 3 fe)wj "N 4 fe>Pj 

Atunci 

w(0) = l =>A 0 =1 
-cp(()) = 0 =^> ^ A , = 0 
' w(l) = 0 =>A 3 +A 2 +A 1 +A 0 =0 
-cp(l) = 0 ^^A 3+ ^A 2+ iA 1= 0 

Deci A 0 = 1, A , = 0. 
Se formează acelaşi sistemul de ecuaţii: 

/ A 2+ A 3= "! 

*2A 2 + 3A 3 = 0' 



pentru care se obţin soluţiile: 



A 0 = 


1 


A 1 = 


0 


*A 2 = 


-3 


A 3 = 


2 
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Ca urmare, funcţia de formă pentru translaţii este: 

N 1 fe)=A3^+A 2 ^ 2 +A 1 ^ + A 0 

Deci 

N 1 fe) = l-3^ 2 +24 3 

Pe de altă parte, funcţia de formă pentru rotaţii este: 

-N 1 fe)=^A 3 ^ + ^A 2 ^ + iA 1 

Adică 
Deci 

Calculul pentru celelalte seturi de funcţii de formă se face 
asemănător. Se observă că funcţiile de formă corespunzătoare planului x-z sunt 
identice cu cele din planul x-y, doar că funcţiile de formă pentru rotirea 9 , 
notate în final cu indicele cp , au semnul „minus" spre deosebire de cele pentru 
rotirea \|/ : 

Funcţii de formă pentru deplasări, valabile atât în planul x-y cât şi în planul x-z: 



N 1 fe) = l-3§ 2 +2§ 3 
N 2 (0 = 4-2^ 3 ) 
N 3 (ţ)=3ţ 2 -2ţ 3 

n 4 ©=l(-^ 2 +4 3 ) 
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Funcţii de formă pentru rotiri în planul x-z. 



N 2(p ft) = 
N 4tp fe) = 



3^ 2 -4^ + l]=-N 2 fe), 
3^ 2 -2^]=-N 4 fe> 



Deci 



KJ=-LnJ 
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Calculul funcţiilor de formă pentru elementul de arbore Timoshenko 



Calculul funcţiilor de formă 
în planul x-y 

Funcţiile de formă au următoarele proprietăţi importante: 

1 . Suma funcţiilor de formă are valoarea egală cu 1 în orice nod; 

2. Corespunzător coordonatei n, funcţia de formă N n =l, iar celelalte 
funcţii de formă au valori nule. 



Pentru Vj = v(x = 0)= vfe = 0) ->N,(£) = l;N,(s)=l; 
Pentru Vl = v (x = 0)=V|/(^ = 0) ^N 2 fe) = l;N 2 fe) = l; 
Pentru Vj = v(x = L)= v(^ = l) ->N 3 ($)=1;N 3 ($)=1; 
Pentru Vj = V (x = L) = ¥ fe = l) ^N 4 fe)=l;N 4 fe) = l. 



relaţii: 



Determinarea funcţiilor de formă Nj(^) şi N x (^) 
în planul x-y 

Pentru determinarea funcţiilor de formă se folosesc următoarele 



De asemenea 



\i/(?)=— A.E 2 + — A 2 E + — f A, + — A, 



vfe)=N 1 fe)v i + N 2 (0v, +N 3 fe)vj +N 4 (0i|/j 
[vfe)=N 1 fe)v i +N 2 fe) M / i +N 3 fe)vj +N 4 feVj 
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Atunci 



v(0) = l 
v|/(0) = 0 
v(l) = 0 
y(l) = 0 



=>A 0 =1 
lf . K 
2 



A 1+^T A 3 



= 0 



A 3 + A 2 + A X + A 0 =0 



3 . 2 . 1 
— A, + — A 2 + — 
L L L 



A 1 + 



K 



v 



= 0 



Se rezolvă sistemul de ecuaţii: 



2A X + 
A 1 + 



KA 3 = 

A. = 



A 2 + 



2A, + 3A. = 



0 
-l 
0 



(A2.1) 



pentru care 



A = 



2 0 K 
111 

0 2 3 



= 2(l + K); 



A 1 = 



0 0 K 
111 

0 2 3 



= -2K; 





2 


0 


K 




2 


0 


0 


A 2 = 


1 


-1 


1 


= -6; A 3 = 


1 


1 


-1 




0 


0 


3 




0 


2 


0 



= 4 



Deci 



A 0 = 
A 1 = 

A 2 = 

A,= 



1 

A 

A2 
A 
Aş 
A 



K 



1 + K 

3 

1 + K 
2 

1 + K 
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Ca urmare 
Adică 



N 1 fe)=A3^ 3 +A 2 ^ 2 + A 1 ^ + A 0 
uv 1 + K^ 1 + K^ 1 + K^ 



Deci 



Ni ® = î7k[ 1+2 ^" 3 ^ 2+k(1 "^ 



Pe de altă parte 



2 3 J 



Adică 



1W L 1 + K^ L 1 + K^ L 



K K 2 > 
, ~1 + K + T'l + K, 



Deci 




Determinarea funcţiilor de formă N 2 (^) şi N 2 (^) 

în planul x-y 



A 0 =0 



= 1 



v(0) = 0 

y(o)=l => if A 1+ Îa 

v(l) = 0 =>A 3 +A 2 +A 1 +A 0 =0 

*Kl) = 0 => l A3 + L A2 + L K 



A i+y A 3 l = o 
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Sistemul de ecuaţii care se formează este următorul: 

A 1 + 



KA 3 = 



A, + 



A 3 = 



2L 
0 



2A 2 + 



3A. = -L 



Pentru acest sistem, determinantul coeficienţilor A = 2(l + K) este acelaşi ca 
pentru sistemul (A2.1), dar ceilalţi determinanţi A 1? A 2 , A 3 sunt diferiţi: 



A 1 = 



2L 0 K 
0 11 
-L 2 3 



= L(2 + K); 



A, = 



2 2L K 
10 1 
0 -L 3 



= -L(4 + K); 



A 3 = 



2 0 2L 
11 0 
0 2 -L 



= 2L 



Ca urmare 



A 0 = 
A 1 = 



A 2 = 



A 3 = 



0 

Aj_ 
A 
A2 
A 
A3 
A 



L(2 + K) 

2(1 +k) 

L(4 + K) 
2(1 +K) 

L 
1 + K 



N 2 © = A3^+A 2 ^ 2 +A^ + A ( 



Adică 



2V ' 1 + K 2(1 + K) 2(1 + K) 



Deci 
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Pe de altă parte 



Adică 



Ct f*\ 3 L 2 2 L(4 + K) £ l(L(2 + K) 



L 1 + K S L 2(l + K) s L^2(l + K) 2 1 + K 



K L 



Deci 



^ 2fe)= î7K^ 2 " 4 ^ +1+K(1 "^ 



Determinarea funcţiilor de formă N 3 (^) şi N 3 (^) 

în planul x-y 



v(0) = 0 
V (0) = 0 
v(l) = l 

v(i)=o 



A 0 =0 
1 f 



K 



Aj + — A 



= 0 



A 3 + A 2 +A 1 + A 0 =1 



-A, + — A, + 



1 



K 



Aj + — A 



= 0 



Sistemul de ecuaţii care se formează este următorul: 

r 2A x + 

Aj+ A 2 + 



KA 3 = 0 
A, = 1 



2A 2 + 3A 3 = 0 



Pentru acest sistem, determinantul coeficienţilor A = 2(l + K) este acelaşi ca 
pentru sistemul (A2.1), dar ceilalţi determinanţi A 1? A 2 , A 3 sunt diferiţi: 



A 1 = 



0 0 K 
111 

0 2 3 



= 2K; 



A 2 = 



2 0 K 
111 
0 0 3 



= 6; 
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A 3 = 



2 0 0 
111 
0 2 0 



= -4. 



Ca urmare 



Adică 



A 0 = 



A 1 = ^ = 



A 2 = 



A 3 = 



0 

Aj_ 
A 

A 

A^ 

A 



K 



1 + K 

3 

1 + K 
2 

1 + K 



N3fe) = A3^+A 2 ^ 2 +A^ + A ( 



N 3 fe) = — — ţ 2 + — £ + 0 

3VV 1 + K^ 1 + K^ 1 + K^ 



Deci 



Pe de altă parte 



N 3 fe)=| A 3 ^ + ^A 2 ^ + Ua i+ ^A 



Adică 



3V9/ L 1 + K^ L 1 + K^ L 



K K 2 



1 + K 2 1 + K 



Deci 
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Determinarea funcţiilor de formă N 4 (^) şi N 4 (%) 

în planul x-y 



v(0) = 0 
V (0) = 0 
v(l) = 0 

v(i)=i 



A 0 =0 

E [A 1+y A 3 J = 0 

A 3 +A 2 + A 1 + A 0 =0 
3 . 2 . lT K 



-A 3 + 



L J L 

Sistemul de ecuaţii care se formează este următorul: 

r 2A x + 

A 1 + 



■A, + 



— A, 



= 1 



KA 3 = 



A, + 



2A 2 + 



3A 3 = 



0 

= 0 
L 



Pentru acest sistem, determinantul coeficienţilor A = 2(l + K) este acelaşi ca 
pentru sistemul (A2.1), dar ceilalţi determinanţi A 1? A 2 , A 3 sunt diferiţi: 



A 1 = 



0 0 
0 1 

L 2 



K 
1 

3 



= -KL; 



A 2 = 



0 K 
0 1 
L 3 



= L(K-2); 



A, = 



0 0 



1 

2 



0 
L 



= 2L. 



A 0 = 
A 1 = 

A 2 = 



0 

Aj_ 
A 

A2 
A 

A3 
A 



KL 



2(1 +K) 
L(K-2) 

2(1 +K) 
L 

1 + K 



Ca urmare 



N 4 (0 = A 3 ^ + A 2 ^ 2 +A^ + A 0 
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Adică 



KL 



N (A= L £ 3 L(K-2) 2 ^ £ Q 

W4 ^~l + K^ 2(1 + K)^ 2(l + Kr + ° 



Deci 




Pe de altă parte 



N 4 (0 = f A 3 £ 2 +| A 2 ^ + f fA 1 +| A 3 



Adică 



~ 3 L K 2 2 L(K-2V 1 



L 1 + K L 2(1 + K) 2(l + K) 2 1 + K 



KL K L 

■ + - 



Deci 



S-au determinat funcţiile de formă pentru deplasări 
Nj(^), N 2 (^), N 3 (^), N 4 (^) care ţin cont de forfecare, deci pentru elementul 

Timoshenko. In expresiile acestora, pentru K=0 se obţin polinoamele Hermite de 
gradul trei, deci pentru elementul Bernoulli-Euler. 

De asemenea, s-au determinat funcţiile de formă pentru rotiri 
Nj©, N 2 (^)> N 3 (^)> N 4 (^) pentru elementul Timoshenko. 

Când K=0 între aceste două seturi de funcţii există relaţia: 
5^)=!^-, (i = l,...,4> 

Este utilă introducerea unui al treilea set de funcţii de formă care se 
folosesc la derivarea matricei de rigiditate şi se definesc astfel: 
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Acestea sunt: 



H=H-r 



dN 



N 2 fe) = 



1 K 



L 1 + K 
1 K 



2 1 + K 



Calculul funcţiilor de formă 
în planul x-z 

Folosind proprietăţile funcţiilor de formă, se calculează expresiile 
acestora şi în planul x-z. 



Determinarea funcţiilor de formă Nj(^) şi N x (^) 

în planul x-z 



relaţii: 



Pentru determinarea funcţiilor de formă se folosesc următoarele 



w(0 = A 3 ^ + A 2 ^ 2 + A^ + A ( 



De asemenea 



wfe) = N 1 fe)w i + N 3 feVj -N 4 fe)pj 

l-q>fe) = N 1 fe)w i -N 2 fe>Pi +N 3 feVj "N 4 fe>Pj 
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Atunci 



w 



(0) = 1 
-cp(0) = 0 
w(l) = 0 
- 9 (1)= 0 



^A 0 =l 

4(a 1+ |a 3 ) = o 

=> A 3 +A 2 + A x + A 0 =0 



3 . 2 . \( 

— A, +— A, + — 

L 3 L 2 L V 



K 



A i + — -3 



= 0 



Se formează acelaşi sistem de ecuaţii (A2.1) ca şi pentru determinarea funcţiilor 
N t fe>N 2 fe) din planul x-y: 

r 2A x + KA 3 = 

Ax+ A 2 + A 3 = 



2A 2 + 3A 3 



0 
1 
0 



pentru care se obţin aceleaşi soluţii: 



A 0 = 
A 1 = 



A 2 = 



1 

A 
A 

A 3 = ± 
3 A 



K 



1 + K 

3 

1 + K 
2 

1 + K 



Ca urmare, funcţia de formă pentru translaţii este: 



N 1 fe)=A3^ 3 +A 2 ^ 2 +A 1 ^ + A 0 



Adică 



Deci 



K 



1 + K 1 + K 1 + K 



N ife) = î^[ 1 + 2 ^- 3 ^ 2+K ( 1 -^] 



Pe de altă parte, funcţia de formă pentru rotaţii este: 



190 



Anexa 2 



-Nxfe) = 



L 1 + K 



7 2 



3 e 1' 

£ + — 

L 1 + K L 



K K 



1 + K 2 1 + K 



Deci 




Calculul pentru celelalte seturi de funcţii de formă se face 
asemănător. Se observă că funcţiile de formă corespunzătoare planului x-z sunt 
identice cu cele din planul x-y, doar că funcţiile de formă pentru rotirea cp au 
semnul „minus" spre deosebire de cele pentru rotirea \|/ : 
Funcţii de formă pentru deplasări, valabile atât în planul x-y cât şi în planul x-z. 



N 



1 ® = î+ 1 k[ 1 " 3 ^ 2+2 ^ + k(1 "^ 



N 2 fe) = 



1 + K 



Lfe-2S 2 +S 3 )+K^-$ 2 ) 



N3ft)= î+ 1 K^ 2 " 2 ^ +K ^ 



1 



1 + K 

Funcţii de formă pentru rotiri în planul x-z 

1 



l(-^ 3 )+k^(-^) 



N 2<P fe)=- î ^N 2 -^ + l + K(l-0]=-N 2 fe), 
N 3, fe) = " [f (" 6^ 2 + 6$)] = -N 3 (4), 



N 4cp fe) = " ^ N 2 " ^ + K J = "N 4 fe> 



Deci 



KJ=-LnJ 

In urma calculelor efectuate pentru determinarea funcţiilor de 
formă pentru cele două tipuri de elemente de arbore, se observă că funcţiile de 
formă pentru elementul Bernoulli-Euler se pot obţine dacă în expresiile 
funcţiilor de formă corespunzătoare elementului Timoshenko factorul K devine 
egal cu zero. Aceeaşi legătură va exista şi între matricele de masă, giroscopică şi 
de rigiditate ce se vor determina pentru cele două tipuri de elemente de arbore. 
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Calculul submatricelor de masă şi de rigiditate 
pentru elementul de arbore Bernoulli-Euler 



Pentru stabilirea formelor explicite ale submatricelor de masă, 
giroscopică şi de rigiditate, este nevoie de calcularea integralelor următoare, 
pentru fiecare tip de element de arbore în parte: 

JLNfe)J T LNfe)jd ? 

o 

o 



dN 


T 


dN 









o 

In cele ce urmează se vor rezolva aceste integrale atât pentru 
elementul Bernoulli-Euler, cât şi pentru elementul Timoshenko (Anexa 4). Se va 
ţine cont de faptul că funcţiile de formă sunt diferite pentru cele două tipuri de 
elemente, şi că pentru elementul Bernoulli-Euler, în formulele generale de mai 
sus, termenii care ţin de inerţia la rotaţie şi de forfecare sunt nuli. 

Submatricea coerentă de masă pentru elementul de arbore Bernoulli-Euler 

Funcţiile de formă pentru elementul Bernoulli-Euler sunt 

următoarele: 
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Funcţii de formă pentru 
rotaţii 



N 2 fe)=3^ 2 -4^ + l 
N 3 fe)=j;(-^ 2 +^) 



Funcţii de formă pentru 
translaţii 

N 2 (0 = L(^-2^ 3 ) 
N 3 (0=3^ 2 -2^ 3 

n 4 ©=l(-4 2 + ^ 3 ) 

[N 4 £)=3^ 2 -2^ 

Din relaţia (3.138), submatricea de masă pentru un element de 
arbore, în cazul general, are forma: 

t» s ]=4LNfe)J T LN(|)Jd| + Ml.(LNfe)J r LNft)k 

0 o 
1 

Se poate calcula integrala: J |_N(£)_| T |_ N fe)J d £ = 



-f 



l-3£ 2 +2^ 3 

4-2^ 2 + ^ 3 ) 

3£ 2 -2£ 3 

_l(-^ 2 H 3 ) J 



•[l-3^ 2 +2£ 3 l(^-2^ 2 +^ 3 ) 3^ 2 -2£ 3 l(-^ 2 +£ 3 )J^ = 



}(l-3^+2^ 3 ) 2 d^ 

0 


jL(l-3^+2^ 3 ). 

0 


j(l_3^ + 2^ 3 ). 

0 

•fe 2 -2^ 


jL(l-3^ 2 + 2^ 3 ). 

0 




0 


Jl(^-2^ + ^ 3 ). 

0 

■(3^-2^ 


jL 2 (|-2^ + | 3 ). 

0 






0 


Jl(3^ 2 -2^ 3 ). 

0 


simetrică 






0 
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420 



156 22L 54 -13L 
4L 2 13L -3L 2 



156 -22L 
sim 4L 2 

Al doilea termen din expresia submatricei de masă este nul pentru elementul 
Bernoulli-Euler, pentru că acesta nu ţine cont de inerţia la rotaţie, deci ji , care 
reprezintă masa de rotaţie pe unitatea de lungime, este egal cu zero. 



Submatricea coerentă de masă pentru elementul de arbore 
Bernoulli-Euler este: 

"156 22L 54 -13L 

î 



sim 



4L 2 13L -3L 2 
156 -22L 
4L 2 



Submatricea de rigiditate pentru elementul Bernoulli-Euler 

Din formula generală a subatricei de rigiditate, pentru elementul 
Bernoulli-Euler, al doilea termen este nul, deoarece acest element nu ţine cont 
de forfecare. 



dN 


T 


dN 









dS + GA s L|LNj[N_k 



Atunci, expresia submatricei de rigiditate pentru acest element este: 

I n im n im 

Derivatele în raport cu \ ale funcţiilor de formă pentru rotaţii sunt: 

hm-*) 



dN 


T 


dN 






_ d ^_ 



dN 



L 

6^-4 

i(-12^ + 6) 
6^-2 
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Se poate calcula integrala: 



J 



dN 


T 


dN 









1 

-J 



6^-4 

±(-12$ + 6) 
6^-2 



^(12^-6) (6^-4) I(- 12^ + 6) (6^-2) d^ = 



jA(l^-6) 2 d$ 

0 ^ 


jf(uţ-«). 

0 ^ 

fe-4)d^ 


o ,v 

(-12^ + 6^ 


ji(!2|-6). 
o ,v 

fe-2)d^ 




J(6^-4) 2 d^ 

0 


J o L 

•(-12^ + 6)d^ 


J(6§-4). 

0 

fe-2)d^ 






K(-12^ + 6) 2 d^ 
J o L 


)f(-124 + 6). 
fe-2)d^ 


simetrică 


"12 6L 


-12 6L " 


J(6^-2) 2 d^ 

0 



4L 2 



sim 



-6L 2L 2 
12 -6L 
4L 2 



Submatricea de rigiditate pentru elementul de arbore Bernoulli- 



Euler este: 



dN 


T 


dN 






_ d ^_ 



d^= 



EI 



12 6L -12 6L 
4L 2 -6L 2L 2 



sim 



12 -6L 
4L 2 
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Calculul submatricelor de masă, giroscopică şi de rigiditate 



pentru elementul Timoshenko 



Submatricea de masă coerentă pentru elementul Timoshenko 

Din relaţia (3.138), submatricea de masă pentru un element de 
arbore, în cazul general, are forma: 

t» s ]=HLjLNfe)J T LNfe)Jd^ALjLNfe)J r LNft)k 



Funcţiile de formă pentru translaţii pentru elementul Timoshenko 
sunt următoarele: 



N 



1 ^ = î7k[ 1 " 3 ^ 2+2 ^ +k(1 "^ 



1 + K 



1 + K 



L (_^3) +K L(_^ 2 ) 



± 

Se poate calcula integrala: J |_N(£)_| T |_ N fe)J d £ = 



[l-3^ 2 +2^+K(l-^)] 
[^-2? + ?) + K±(ţ-?)] 

[l(-^h 3 Kk|(-^ + ^)] 



[ 1 _3^ + 2^ + K(l-^)] [l(^-2^ + ^)+k|(^-^)] [3^ 2 -2^ + K^] + ^)+k(-^ + ^)] 
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-+- * 

r 


^ jc 

i * jjj^ 

•JLX" -+- 

+ ~zr^ 
+ 


1 

sa » 

-+- +^ 
•» «»» 


1 

M 

1 

M 

+- 










*f -+• 

ST V 


*J | *-* 

M 

■+■ -4- 
■ 1 


SC 




"«_#» 

■+■ 

^ 1 


1 

-J 1 

M 

■+- 

1 






1 

+- 

-+- 






1 
■I 

i 
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Al doilea termen din expresia submatricei de masă ţine cont de 
funcţiile de formă pentru rotaţii. 

Funcţiile de formă pentru rotaţii pentru elementul Timoshenko sunt 

următoarele: 



N 2 ft) = î ^W-4Ş + l + K(l-|)] 



Se poate calcula, apoi, integrala: 



(1 + K) 2 I 



[(l-4^ + 30+K(l-0] 
I(-6^ + 60 
_^ 2 -2^ + K^) 



[i(6^ 2 -60] [(l-4$ + 3^ 2 ) + K(l-0] [±-(-6ţ 2 + 6ţ] + 
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Atunci, submatricea de masă pentru elementul de arbore Timoshenko are 
următoarea expresie: 



k]= 



(1 + K) 2 



420 



156 22L 54 
4L 2 13L 
156 



-13L 
-3L 2 
-22L 
4L 2 



+ K 



(a,L 



120 



84 11L 

2L 2 



36 
9L 
84 



-9L 
-2L 2 
-11L 

2L 2 



+ K 



120 



40 5L 20 
L 2 5L 2 
40 



-5L 
-L 2 
-5L 

L 2 



+ 



30L 2 



36 3L 

i 4L 2 



-36 3L 
-3L -L 2 
36 -3L 
4L 2 



+ 



+ K 



pL 



6L 2 



0 - 



3L 0 
L 2 3L 
0 



-3L 
-L 2 
3L 

L 2 



K 



jlL 



6L 2 



0 0 

; 2l 2 



o o 

0 L 2 

o o 

; 2l 2 



Sau altfel scrisă, 



unde: 



m 



ii 



sim 



m 
m 



s 

12 

S 
22 



m 
m 
m 



s 




13 




S 




23 


m s 24 


S 


s 


33 


m 34 








m 44 



= (l56 + 294K + 140K 2 )x T + 36a R 
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mj 2 =(22 + 38.5K + 17.5K 2 )La T +(3-15K)La R 

= (54 + 126K + 70K 2 )x T - 36<x R 
m s 14 =-(l3 + 31.5K + 17.5K 2 )La T +(3-15K)La R 
m L = (4 + 7K + 3.5K 2 )L 2 a T + (4 + 5K + 10K 2 )L 2 a R 

m 24 = -(^ + 7K + 3.5K 2 )L 2 a T - (l + 5K - 5K 2 )L 2 a R 

m 33= m il 

m 34 = ~ m n 
m s 44 = m s 22 

unde am făcut notaţiile: 
(0.L 

T ~ 420(l + k) 2 

Şi 

_ Al 

aR ~30L 2 (l + K) 2 



Submatricea giroscopică pentru elementul de arbore Timoshenko 

Din relaţia (3.138), submatricea giroscopică pentru un element de 
arbore, are forma: 
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[g s ]=ApLjLNfe)J I LNfe)k = 2tn R S ], 

o 

deci preluând din relaţia matricei de masă (3.160) determinată anterior, termenii 
care conţin factorul p.L şi înmulţindu-i cu doi, se obţine submatricea 
giroscopică pentru elementul de arbore Timoshenko: 



b']= 



(1 + Kf 



30L 2 



36 3L 

i 4L 2 



-36 3L 
-3L -L 2 
36 -3L 
4L 2 



+ 



+ K 



6L 2 



0 -3L 0 -3L 
L 2 3L -L 2 
0 3L 



+ K 



6L 2 



0 0 0 0 

1 2L 2 0 L 2 

0 0 

; 2l 2 



sau 



'ii 



>12 
>22 



§33 



sim 



§14 
§24 
§34 
§44 



în care, cu aceleaşi notaţii făcute anterior a x şi a R , elementele submatricei 

au următoarele valori: 

§ii=72a R 

g^ 2 =2(3-15K)La R 
§13 = "72a R 



g J 4 =2(3-15K)La R =gJ 2 
g^ 2 =2(4 + 5K + 10K 2 )L 2 a R 
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§23 — §12 



g^ 4 =-2(l + 5K-5K 2 )L 2 a R 



§33 — §11 



§34 — §23 



§44 — §22 



Submatricea de inerţie pentru elementul de arbore Timoshenko 



Din relaţia (3.155), submatricea de inerţie pentru un element de 
arbore, are forma: 



dN 


T 


dN 









dS + GA s LjLNfLNk. 



Pentru calculul acestei submatrici este nevoie de : 

Derivatele în raport cu ţ, ale funcţiilor de formă pentru rotaţii ale elementului 
Timoshenko: 





1 




1 + K 


dN 2 






1 + K 


dN 3 


1 




1 + K 


dN 4 






1 + K 



hm-*) 



Se poate rezolva integrala: | 



o L 



dN 


T 


dN 


_d£_ 







d^= 
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^-6) 

(l + K) 2 J I(_i24+6) 
6^-2 + K 

1(12^-6) (64-4-K^) ^(-12^ + 6) (6^-2 + K) 



(1 + Kf 



iA( i2 ^-6) 2 ^ 

0 L 


ji(l2Ş-6). 

0 

•(6^-4-K^ 


0 

• (-12^ + 6^ 


ji(l2ţ-6). 

0 

.te-2 + K)N; 




J(6^-4-K^) 2 d^ 

0 


}l(6^-4-K0- 

0 

•(-12^ + 6)d^ 


}(6^-4-K4). 

0 






}A(- 12 ^ + 6) 2 d4 

0 L 


jf(-12ţ + «). 

0 

•(6£-2 + K)d£ 


simetric 






j(6^-2 + K) 2 d^ 

0 



Este nevoie, de asemenea, de diferenţele dintre funcţiile de formă pentru rotaţii 
şi derivatele în raport cu \ ale funcţiilor de formă pentru translaţii ale 
elementului Timoshenko: 

K 1 

= 

1 1 + K L 

1 K 

N 2 = 

2 1 + K 

K 1 

N, = — - 

3 1 + K L 

1 K 

. 4 2 1 + K 



pentru a rezolva integrala: 



206 



Anexa 4 






iu* 

0 








K 2 




0 






(1 + K) 2 






0 






simetrică 






0 



Atunci, forma finală a submatricei de rigiditate pentru elementul de arbore 
Timoshenko este: 



i i:i 

l + K V 



12 6L 

; 4l 2 



-12 6L 
-6L 2L 2 
12 -6L 
4L 2 



+ K 



0 0 0 

; l 2 o 
o 



o 

-L 2 
O 
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